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1.  INTRODUCTION  AND  STATEMENT  OF  THE  PROBLEM 


The  evaluation  of  the  potential  performance  of  gas-driven  gun  propulsion  concepts  requires  an  accurate 
determination  of  the  thennochemical  properties  of  the  propelling  gases.  The  principal  quantities  of  interest 
are  as  follows:  (1)  the  effective  gas  temperature  (flame  temperature  in  the  absence  of  electrical  energy); 
(2)  the  ratio  of  the  frozen  specific  heats  (y);  (3)  the  covolume;  and  (4)  the  ballistic  energy,  which  is 
defined  by 


Ballistic  Energy  =  .  (1) 

Thennochemical  codes  designed  for  ballistic  applications  (i.e.,  utilizing  nonideal  equations  of  state) 
such  as  BLAKE  (Freedman  1982).  F.  Volk’s  ICT  code  (Volk  and  Bathelt  1978).  or  SNFE’s  BAGHEERA 
(Kotlar  1991)  have  proven  very  effective  in  calculating  these  quantities  for  standard  solid  propellant  gun 
systems.  Unfortunately,  calculating  thermochemical  values  for  gases  produced  by  the  plasma/propellant 
mixture  in  the  electrothermal-chemical  (ETC)  ballistic  process  is  complicated  by  die  likelihood  of 
achieving  gas  temperatures  well  in  excess  of  5,000  K. 

Two  potential  difficulties,  which  can  result  in  incorrect  thennochemical  values  being  calculated,  arise 
when  die  gas  temperature  greatly  exceeds  5,000  K.  First,  a  significant  concentration  of  ions  may  be 
formed  in  die  gas  at  these  temperatures;  in  their  present  forms,  neither  BLAKE  nor  die  ICT  code  takes 
into  account  the  presence  of  ions  and  thur  concomitant  effect  on  the  thermodynamics. 

Second,  die  representations  of  the  temperature  dependence  of  the  thermodynamic  data  used  by  the 
codes  in  performing  the  calculations  may  not  be  valid  much  beyond  this  temperature.  This  point  requires 
further  discussion.  General-purpose  thermodynamics  codes,  such  as  the  ones  mentioned,  must  have 
available  as  functions  of  temperature  the  thermodynamic  quantities,  heat  capacity  at  constant  pressure  (Cp), 
enthalpy  (H).  and  entropy  (S),  for  500-1,000  species.  In  order  to  minimize  computer  memory 
requirements,  these  quantities  are  not  themselves  stored  directly.  Since  H  and  S  are  directly  obtained  by 
integration  of  C^(T)<fr  and  [Cyrjr/T,  respectively,  it  is  sufficient  to  store  only  C^,  plus  two  integration 
constants.  Furthermore,  even  Cp  itself  is  not  stored  directly,  instead,  die  coefficients  of  an  empirical 
representation  are  used. 
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There  is  no  objective  way  of  deciding  which  representation  is  the  best  Gordon  and  McBride  (1971) 
used 


CpOV*  *  c  +  <f*T  +  e*'!'2  +  f*T®  +  g*T*  . 


(2) 


More  recently,  McBride  (1992)  has  been  using 


CpCry*  =  a*T2  +  b*T*  +  c  +  d*T  +  e^T2  +  frT3  +  g*!4 .  (3) 

The  authors  of  the  TIGER  code  (Cowpeithwaite  and  Zwisler  1973),  and  subsequently,  BLAKE,  used 

CpCW  =  B(l)  +  B(2)*0  +  BO)*©2  +  B(4)*03 

+  B(5)*(l/0)  +  B(6)*(l/0)2  +  BCD*(l/9)3  .  (4) 

In  Eqs.  2-4,  T  is  the  absolute  temperature  and  St  is  the  universal  gas  constant  in  appropriate  units 
(originally,  calories/mole-K;  more  recently,  joules/mole-K).  In  Eq.  (4),  6  is  T/ljOOO. 

In  all  cases,  the  constants  a-g  or  B(l)-B(7)  are  determined  by  linear  least  squares  fitting  to  the  heat 
capacity  data  in  the  JANAF  Tables  (Chase  et  al.  1986).  Earlier  editions  of  die  JANAF  Tables  terminated 
at  5,000  K  and.  therefore,  so  did  die  fittings  in  BLAKE  and  the  ICT  programs.  Extrapolating  linear  least 
squares  fittings  is  known  to  be  a  risky  procedure.  (This  difficulty  can  often  be  avoided  [Chase  1985]  by 
including  an  estimated  heat  capacity  at  a  temperature  well  outside  die  upper  range,  but  this  was  not  done 
in  the  BLAKE  fittings.)  Any  inaccuracies  arising  from  such  extrapolation  will  be  propagated  through  die 
calculations,  possibly  resulting  in  significant  errors  in  the  values  determined  for  die  other  thermochemical 
properties. 

Even  when  die  total  amount  of  added  electrical  energy  is  relatively  low  so  drat  die  overall  average 
gas  temperature  remains  below  5,000  K,  transient  high  gas  temperatures  (over  5,000  K)  with  their 
(possibly)  resulting  inaccurate  thermochemical  properties  may  occur.  (The  time  scale  of  these  transients, 
while  short  compared  to  the  gun  firing  cycle,  is  nevertheless  long  enough  so  that  thermochemical 
equilibrium  is  maintained  in  the  chamber.)  This  is  illustrated  in  Figures  1  and  2.  Figure  1  shows  the 
calculated  instantaneous  electrical  energy  density  (kilojoules  of  input  electrical  energy  per  gram  of 
propellant  consumed)  for  an  experimental  30-mm  ETC  gun  firing.  One  sees  electrical  energy  densities 
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greater  than  five  occurring  for  about  1.2  ms  during  the  first  1.4  ms  of  the  ballistic  cycle.  Figure  2  shows 
that  electrical  energy  densities  greater  than  five  (for  the  M30  propellant  used  in  die  experiment)  result  in 
calculate  gas  temperatures  in  excess  of  5,000  K.  Thus,  even  though  the  average  plasm  a/propellant  gas 
temperature  is  approximately  3,800  K  (using  an  average  electrical  energy  density  of  1.16  kJ/g)  during  a 
significant  portion  of  the  ballistic  cycle  (1 2  ms  of  die  total  6.4  ms,  Figure  1),  the  effective  gas 
temperature  would  exceed  5,000  K.  Wren  and  Oberie  (1993)  have  shown  that  ignoring  the  effects 
with  the  high  transient  gas  tempera. ires  and  varying  thennochemical  properties  can  result  in 
significant  differences  in  interim'  ballistic  simulations  and  in  inverse  analyses  of  experimental  data  for 
ETC  systems. 

The  objective  of  the  present  work  was  to  produce  a  unified  thermodynamics  code  and  its 
accompanying  library  that  can  prod  *.  all  of  the  customary  ballistic  thennochemical  data  for 
ETC  propellant  systems  at  tempei.  ires  up  to  10,000  K  and  preferably  15,000  K.  Specifically  discussed 
are  the  following:  (1)  the  overall  approach  and  rationale  to  producing  the  desired  thennochemical  code; 
(2)  evaluation  of  the  effect  of  including  ions  on  calculated  thennochemical  values;  (3)  modifications  to 
the  BLAKE  thermodynamics  code;  and  (4)  results  of  several  test  cases  utilizing  the  unified 
thennochemical  code  with  libraries  over  different  temperature  ranges. 

2.  APPROACH 

In  order  to  achieve  the  desired  objectives,  the  following  steps  were  adopted:  (1)  determine  if  the 
presence  of  ions  at  high  temperatures  significantly  affects  values  for  computed  thennochemical  data; 
(2)  based  upon  the  results  of  (1),  decide  which  thermodynamic  code  to  modify;  and  (3)  modify  the 
selected  code. 

2.1  nisriwsinn  of  Step  1.  To  determine  the  potential  effect  of  ions,  the  National  Aeronautics  and 
Space  Administration  Lewis  Research  Center’s  (NAS  A-Lewis)  family  of  codes  (CET86,  TRAN72,  CET89, 
and  variants)  (Gordon  and  McBride  1971;  Svehla  and  McBride  1973)  was  selected  for  use  since  they  can 
easily  work  with  ions.  Before  evaluating  the  effect  that  ions  may  have  on  the  computed  thennochemical 
properties,  however,  it  was  necessary  to  obtain  thermodynamic  data  for  species  produced  in  typical  ETC 
gun  propellants  up  to  10,000  K.  Fortunately,  Bonnie  J.  McBride  of  NASA-Lewis  in  Cleveland,  OH,  was 
able  to  provide  to  the  Army  Research  Laboratory  a  set  of  coefficients  that  would  permit  computations  up 
to  15,000  K  (and  20,000  K  in  some  cases)  for  133  neutral  and  ionic  species  containing  C-H-O-N. 
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Since  ballistic  energy  (defined  earlier)  is  the  quantity  actually  utilized  in  most  ballistic  simulations, 
it  was  selected  as  the  figure  of  merit  for  ascertaining  the  effect  of  ions.  A  difference  of  more  than  2% 
in  computed  ballistic  energy  (with  and  without  the  formation  of  ions  permitted)  was  selected  as  the 
decision  criterion.  Results  are  presented  in  the  next  section. 

2.2  Discussion  of  Step  2.  We  consider  only  two  codes,  BLAKE  and  CET86,  the  next-to-latest 
version  of  the  NASA's  Lewis  Research  Center’s  chemical  equilibrium  and  transport  codes.  Each  code 
has  its  own  advantages  and  disadvantages.  CET86  is  an  outstanding  program  that  is  easy  to  nm  and 
readily  works  with  ions.  Unfortunately,  it  utilizes  only  die  ideal  gas  equation  and,  therefore,  cannot 
estimate  covolumes,  a  necessity  for  gun  calculations.  If  inclusion  of  ions  proved  to  be  necessary,  then 
the  CET86  code  would  be  adopted,  and  a  modification  or  separate  program  would  be  developed  to 
compute  covolumes.  On  die  other  hand,  if  die  formation  of  ions  is  not  significant  at  temperatures  up  to 
10,000  K,  then  the  BLAKE  code  would  be  used,  and  an  extended  library  (valid  up  to  at  least  10,000  K) 
would  be  prepared  for  it 

2.3  Discussion  of  Step  3.  As  will  be  shown,  the  formation  of  ions  was  found  not  to  be  important 
(at  least  for  die  propellant  systems  considered),  and  consequently  the  BLAKE  code  was  selected.  Only 
a  summary  of  the  modifications  to  BLAKE  will  be  provided  in  this  report;  more  detailed  information  will 
be  contained  in  a  forthcoming  report 

3.  EFFECT  OF  IONS 

As  mentioned  previously,  die  NASA-Lewis  code  CET86  was  utilized  to  determine  die  effect  of  ions 
on  the  calculated  thermochemical  data.  A  variant  of  CET86  called  MuCET*  was  chosen  for  this  task. 
Unfortunately,  the  format  of  McBride's  high-temperature  thermochemical  data  differed  hum  that  used  in 
all  NASA-Lewis  thermochemistry  programs  through  CET89  and,  hence,  MuCET.  It  was  therefore 
necessary  to  produce  a  modified  version  of  MuCET  to  accommodate  the  new  format  for  the 
high-temperature  thermochemical  data;  this  modified  version  of  MuCET  was  named  EXMuCET. 
Comparisons  between  MuCET  and  EXMuCET  at  temperatures  below  5,000  K,  using  the  same 
thermochemical  library,  were  performed  to  ensure  the  accuracy  of  EXMuCET.  Overall,  die  comparisons 
were  excellent,  but  small  differences  in  computed  values  were  observed  at  temperatures  below  about 


*  MnCET  is  •  version  of  NASA’s  CET86  adapted  for  easy  me  with  IBM-cofnpetible  desktop  computers.  It  was  ptodneed  by 
Eli  Freedman  it  Associates  md  is  licensed  to  the  U.S.  Army. 
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3,000  K;  these  were  caused  by  the  occurrence  or  nonoccurrence  of  condensed  phases.  Since  EXMuCET 
is  intended  for  use  primarily  for  temperatures  above  3,000  K.  the  differences  associated  with  inclusion  or 
noninduskm  of  condensed  phases  are  not  considered  important;  their  presence  is  ignored  in  all  of  the 
subsequent  comparisons  discussed  here.  Some  standard  propellants  (e.g.,  M30)  contain  inorganic  additives 
that  result  in  the  formation  of  condensed  phases  at  chamber  conditions.  When  necessary,  these 
formulations  have  had  such  additives  deleted  from  their  compositions  when  entered  into  the  codes;  they 
are  referred  to  as  pseudopropellants.  Additional  details  can  be  found  in  the  appendix. 


To  determine  the  importance  of  ions  cm  the  thermodynamic  calculations,  eight  propellant  systems 
(listed  in  Table  1)  were  investigated  with  the  formation  of  ions  both  permitted  and  forbidden  over  a  wide 
range  of  electrical  energy  inputs,  which  produced  computed  temperatures  in  excess  of  13,500  K. 


Table  1.  Propellant  Systems  Investigated  With  and  Without  Ions  Permitted 


1.  Pseudo  M30 

2.  WC890 

3.  JA2 

4.  M9 
3.  Ml 

6.  Decalin  +  white  turning  nitric  add 

7.  Decalin  +  70% 

8.  Red  fuming  nitric  add  +  unsymmetrical  dimethyl  hydrazine 


The  differences  in  ballistic  energies  for  foe  two  conditions  (ions  included  vs.  no  ions  induded)  for 
pseudo-M30  and  JA2  over  a  wide  range  of  dectrical  energy  densities  (kilojoules  of  electrical  energy  added 
per  gram  of  propellant)  are  shown  in  Figure  3.  This  figure  shows  foal  there  is  virtually  no  difference  up 
to  10,000  K  (solid  square),  corresponding  to  dectrical  energy  densities  of  22-23  kJ/g.  The  results  are 
similar  for  the  other  six  propellant  systems  investigated.  The  maximum  percent  differences  in  ballistic 
energy  for  the  eight  propellant  systems  at  each  electrical  energy  density  (foe  propellant  system  producing 
foe  maximum  percent  difference  may  differ  at  each  dectrical  energy  density)  are  displayed  in  Figure  4. 
(See  Appendix  G  of  foe  overall  appendix  for  details  on  all  eight  propellant  systems.)  As  would  be 
expected,  this  curve  is  similar  in  appearance  to  foe  ballistic  energy  difference  curve  of  Figure  3. 

In  fact,  for  temperatures  approaching  12,000  K  (electrical  energy  density  of  34  kJ/g),  foe  maximum 
percent  difference  is  less  than  1.0%.  Thus,  based  upon  foe  criteria  stated  in  foe  previous  section,  the 
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inclusion  of  ions  in  tee  thermodynamic  calculations  is  unimportant.  Therefore,  the  program  of  choice  to 
modify  for  ETC  propellants  is  BLAKE  with  a  suitably  prepared  library  valid  over  the  range 
3.000-10000  K  (or  higher). 

The  contents  of  the  new  library  (which  will  be  called  die  extended  library  here,  in  distinction  to  the 
standard  library)  were  limited  by  die  following  two  constraints:  BLAKE  cannot  deal  with  more  than 
29  gaseous  species  in  any  one  run,  and  there  were  no  data  for  compounds  of  potassium  or  magnesium  (the 
two  elements  most  commonly  occurring  in  inorganic  additives)  in  McBride’s  files.  Therefore,  the  new 
library  was  limited  to  31  gases,  and  no  condensed  phases  were  included.  These  31  gases  comprise 
29  species  containing  C,  H,  N,  or  O,  or  combinations  of  them,  plus  helium  and  argon.  The  limit  of 
29  gases  means  that,  in  the  few  cases  in  which  either  helium  or  argon  is  included  in  the  chamber,  a  REJect 
instruction  would  be  needed.  To  avoid  this  requirement  (which  is  unpopular  with  some  users),  die 
CONstituent  instructions  for  the  species  O3  and  HjOC  have  been  displaced  by  six  columns  (resulting  in 
their  omission),  but  their  STR  instructions  have  been  left  intact  These  species  are  almost  never  significant 
in  ETC  computations.  If  a  need  for  either  of  them  should  arise,  it  is  a  simple  matter  to  restore  diem. 

4.  NBLAKE  AND  ITS  LIBRARY 

The  revised  program  has  been  temporarily  named  NBLAKE  to  make  it  easier  for  testers  to  keep  both 
the  old  and  new  programs  in  their  computers.  When  testing  is  completed,  the  name  will  revert  to  BLAKE, 
version  218.0. 

The  only  major  change  to  the  program  concerns  the  way  it  handles  the  library.  The  primary  concern 
was  to  attempt  to  minimize  the  chance  of  using  the  standard  library  in  place  of  the  extended  library  and 
vice  versa.  The  selection  of  the  library  is  made  by  the  user  prior  to  running  die  program,  either  by 
explicitly  copying  die  appropriate  library  from  storage  to  die  library  file,  or  (preferably)  by  choosing  one 
of  two  batch  files  (each  specific  to  a  particular  library)  that  accomplishes  the  same  result  The  only  way 
that  users  will  be  aware  that  they  have  made  the  wrong  choice  will  be  if  they  attempt  to  use  the  extended 
library  on  a  formulation  containing  elements  other  than  C,  H,  O,  or  N  (the  program  will  immediately 
abort).  It  seemed  desirable,  at  the  least,  to  print  information  identifying  the  library  that  was  used.  This 
information  is  now  stored  in  the  first  line  of  die  alphsnumeric  file  of  each  library.  As  s  result,  die  new 
library  cannot  be  used  with  any  preceding  version  of  BLAKE,  and  older  versions  of  the  library  will  not 
work  with  die  new  program. 
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A  few  revisions  were  made  in  the  standard  library  as  follows:  the  data  for  the  condensed  phases  of 
KjS,  KjCOj,  and  KjSC^  were  revised:  the  species  HNO,  Aia2,  and  AICI3  were  added;  and  the  species 
CjH  and  C^N  were  removed. 

Some  lesser  changes  were  made  to  the  program  as  follows.  (1)  The  output  from  the  COMpodtkm 
instruction  has  been  reformatted  to  fit  into  80  columns.  (2)  Many,  but  not  aH.  of  the  messages  printed 
by  die  program  have  been  changed  from  all  UPPER  CASE  to  the  conventional  mixture  of  Upper  rod 
Lower  Cue.  (3)  The  DATe  instruction  is  no  longer  active.  The  program  will  ignore  it  if  used  but  will 
not  issue  an  error  message.  (4)  The  CMT  (CoMmenT)  instruction  has  been  deleted.  The  program  will 
ignore  it  if  used  but  will  not  issue  an  error  message.  (5)  The  FORmula  instruction  has  been  changed  to 
permit  a  formula  name  to  be  redefined  in  foe  same  run.  (6)  The  default  of  foe  ECHo  instruction  has  been 
changed  from  "ofT  to  "on."  Also,  the  appearance  of  foe  echoed  instructions  has  been  changed.  (7)  The 
program  can  now  optionally  read  a  file  containing  ingredient  dma.  (8)  Some  changes  have  been  made  in 
the  list  of  prestored  ingredients;  some  other  changes  had  been  made  prim  to  the  present  work.  (9)  The 
FORmula  instruction  was  modified  to  permit  foe  optional  entry  of  energies  in  joulesAnole  instead  of  foe 
default  kilocalories/mole.  (10)  The  GUN  instruction  was  modified  to  permit  foe  direct  entry  of  electrical 
energy  into  foe  system. 

Further  details  on  all  of  these  changes  will  be  found  in  the  program  documentation  being  prepared. 

5.  TEST  CASES  AND  RESULTS 

Once  the  program  was  modified,  two  important  questions  remained:  (1)  Did  foe  extended  temperature 
library  produce  equivalent  results  with  the  standard  library  for  overlapping  temperature  ranges 
(3,000-5,000  K)?  and  (2)  To  what  extent  do  the  calculated  thermochemical  values  change  in  using  the 
extended  temperature  library  vs.  the  standard  library?  To  address  these  questions,  calculations  were 
performed  using  the  two  thermochemical  data  libraries  with  varying  electrical  energy  densities  (equivalent 
to  varying  the  temperature)  for  five  propellant  systems.  Three  of  the  systems,  pseudo-M30,  WC890,  and 
Decalin/70%  HjC^.  were  exothermic,  while  the  other  two,  methanol  rod  water,  were  endothermic.  Results 
for  three  of  foe  propellant  systems  will  be  presented  in  this  paper  to  illustrate  the  effects  of  calculations 
using  the  different  libraries. 
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The  percent  difference  in  ballistic  energy  is  computed  using  the  following  formula: 

^  ivtf  _  B.E.  (Standard  Lib)  -  B.E.  (Extended  Lib)  A  1fift  m 

B.E.  (Standard  Ub) 

Temperature  and  ballistic  energy  results  at  selected  electrical  energy  densities  for  paeudo-M30  are 
presented  in  Table  2  with  the  percent  difference  in  ballistic  energy  drown  in  Figure  5. 

Table  2  and  Figure  5  rfrow  that  for  pseudo-M30  in  the  overlapping  temperature  range  (3,000-5,000  K) 
foe  differences  in  computed  temperatures  and  ballistic  energies  are  less  than  1%.  Both  the  standard  and 
extended  BLAKE  libraries  use  foe  same  fitting  formula.  The  standard  library  is  essentially  the  one  that 
has  been  used  with  BLAKE  for  years;  it  was  based  on  JANAF  data  that  were  current  at  foe  time.  The 
extended  library,  on  foe  other  hand,  is  new,  it  is  based  on  McBride’s  independent  results.  Considering 
foe  difference  in  the  sources  of  foe  data  for  foe  two  libraries,  we  feel  that  foe  computed  thermodynamic 
properties  are  virtually  foe  same  in  the  overlapping  temperature  ranges.  Thus,  the  libraries  produce 
equivalent  results,  as  one  would  hope,  in  foe  overlapping  temperature  range.  The  percent  differences  in 
ballistic  energies  above  5,000  K,  where  the  libraries  do  not  have  overlapping  ranges,  are  not  significantly 
larger  in  magnitude  for  temperatures  up  to  10,000  K  (2.4%  absolute  maximum)  but  increase  in  magnitude 
to  exceed  +6%  as  the  temperature  increases.  Therefore,  for  pseudo-M30,  the  use  of  the  extended 
temperature  library  does  not  appear  to  have  a  significant  impact  on  calculated  thennochemical  values  until 
gas  temperatures  exceed  10,000  K.  However,  these  types  of  results  are  very  dependent  on  the  propellant 
system,  as  illustrated  by  foe  final  two  propellant  systems  discussed  in  this  paper. 

Results  from  similar  calculations  for  Decalin  plus  70%  hydrogen  peroxide  are  given  in  Table  3  and 
Figure  6.  As  can  be  seen  from  the  table  and  figure,  foe  differences  in  foe  computed  values  are  quite 
significant  outside  of  foe  overlapping  temperature  range  (below  3,000  K  and  above  SjOOO  IQ  with  percent 
differences  in  ballistic  energy  up  to  25%  (less  than  2%  in  foe  overlapping  temperature  range). 
Temperature  differences  are  not  as  severe  as  for  ballistic  energy.  This  propellant  system  animates  the 
potential  errors  which  can  arise  when  foermochexnical  data  are  extrapolated. 

Both  previous  propellant  systems  discussed  were  for  exothermic  propellants.  To  determine  foe  effect 
of  using  the  different  libraries  for  an  endothermic  propellant  system,  the  final  system  discussed  win  be 
methanol-  Results  are  provided  in  Table  4  and  Figure  7.  The  differences  in  this  case  are  even  more 
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Table  3.  Computed  Temperature  and  Ballistic  Energy  Results  for  Decalin  Plus  70% 
Hydrogen  Peroxide  Using  the  Standard  and  Extended  Temperature  Libraries 
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Table  4.  Computed  Temperature  and  Ralliuric  Energy  Results  for  Methanol  Using  the  Standard 
and  Extended  Temperature  Libraries 
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pronounced  than  in  the  previous  two  ctm.  Such  differences  are  apparently  more  pronounced  with 
carbon-deficient  synems.  The  reason  for  this  may  be  explored  at  a  later  dale. 

6.  SUMMARY 

Results  presented  hi  this  paper  have  shown  that  the  formation  of  ions  can  he  safely  ignored  in 
thennodynamic  computations  for  ballistic  applications  at  temperatures  up  to  at  least  10,000  K.  Using  data 
from  NASA’s  Lewis  Research  Center,  a  new  library  for  BLAKE  for  use  at  temperatures  between  3,000  K 
4  at  least  10,000  K  has  been  created  Some  minor  modifications  have  been  made  in  the  program  to 
ittaie  its  use  with  ETC  systems.  Comparison  of  results  below  3,000  K  and  above  5,000  K  show  some 
significant  differences  in  computed  thennochemkal  data  between  the  two  libraries,  especially  with 
carixm-defident  propellant  systems. 
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APPENDIX; 

PREPARATION  OF  A  THERMODYNAMICS  PROGRAM 
AND  LIBRARY  FOR  USE  WITH  ETC  SYSTEMS 


Aside  from  minor  format  changes,  this  appendix  is  presented 
in  its  original  form  without  editorial  changes  or  comments. 
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ABSTRACT 


The  goal  of  this  raft  was  to  produce  a  thermodynamics  program  and  its  accompnying  library  for  use 
with  electrothermal-combustion  propellant  (ETC)  systems  at  temperatures  up  to  10,000  K.  It  was  first 
neoeasary  to  decide  whether  ions  had  to  be  taken  into  account  in  such  a  program.  For  this  propose,  a  new 
program  called  EXMuCET  was  produced,  h  is  a  modification  for  use  at  temperatures  above  WXX)  K  of 
an  earlier  program  called  MuCET,  a  microcomputer  version  off  the  NASA-Lewis  thermochemistry  program 
CET86.  This  program  was  applied  to  9  different  propellant  mixes;  the  results  showed  that  the  formation 
of  ions  could  be  safely  neglected  at  temperatures  up  to  10,000  K  at  the  gas  densities  to  be  expected  in 
ETC  formulations.  This  meant  that  ‘Blake’  would  be  the  program  used. 

A  ‘Blake*  library  for  the  temperature  range  4,000  K-l 0,000  K  was  created  and  tested.  The  testing 
consisted  of  comparing  results  from  the  new  library  with  those  obtained  using  the  accepted  standard 
library  in  the  overlap  temperature  range,  4,000  K-5,000  K.  The  results  of  die  testing  depend  to  a  small 
extent  on  whether  the  REJect  instructions  are  used.  Overall  the  differences  in  impetus  between  the  two 
sets  of  results  are  less  than  03%. 
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I.  INTRODUCTION 


A.  Background. 

1.  An  clcctrothcnnal -combustion  (ETC)  gun  derives  put  (bin  not  usually  all)  of  its  energy  from  a 
high-power  electrical  discharge.  One  of  the  consequences  of  this  discharge  is  that  initially  the  chamber 
temperatures  are  very  high  (10,000-20,000  K),  which  may  cause  ionization. 

2.  The  computation  of  die  performance  of  any  gun  requires  for  its  input  thermodynamic  data  about 
its  propellant.  The  ‘Blake’  code  is  one  of  the  most  widely  used  programs  for  such  calculations,  but  it  has 
drawbacks  for  application  to  ETC  problems: 

a.  The  heat  capacities  of  the  products  in  its  library  are  fitted  as  a  function  of  temperature  only 
up  to  5,000  K. 

b.  The  empirical  formula  that  it  uses  for  fitting  is  not  well-suited  fin  covering  the  entire 
temperature  range  from,  say,  500  K  to  10,000  K. 

c.  This  code  cannot  deal  with  systems  containing  ions. 

B.  Goals.  The  objective  of  this  task  was  to  provide  a  unified  thermodynamics  code  that  can  produce  all 
of  the  customary  thennochemical  data  for  candidate  electrothermal  propellants.  This  code  is  to  have  the 
following  features: 

1.  It  will  be  a  unified  thermodynamics  code  that  produces  all  of  die  customary  thennochemical  data 
for  candidate  electrothermal  propellants. 

2.  It  will  cover  the  temperature  range  from  1,000  K  to  10,000  K  or  higher. 

3.  It  will  furnish  die  frozen  gamma  of  the  equilibrium  composition. 

4.  It  will  estimate  die  covolume  of  the  equilibrium  system. 

H.  APPROACH 

A.  Possible  Codes.  There  are  essentially  only  two  choices  for  the  code,  either  ‘Blake'  or  the 
NASA-Lewis  code,  CET89  (or  variants  of  it  or  of  its  predecessor,  CET86). 

1.  CET89  is  an  outstanding  program  with  major  advantages  and  major  disadvantages. 

a.  It  can  easily  work  with  ions. 

b.  BUT:  It  works  only  for  ideal  gases,  and  therefore  cannot  estimate  covolumes.  If  it  were  to 
be  adopted,  a  modification  or  a  separate  program  would  have  to  be  developed  to  supply  this  deficiency. 

2.  In  most  respects  ‘Blake’  does  not  nm  as  easily  as  CET89,  but  it  was  designed  for  use  primarily 
with  non-ideal  gases  and  the  propellant  covolume  is  one  of  its  primary  outputs.  If  it  can  be  shown  that 
the  formation  of  ions  is  not  of  importance  at  temperatures  up  to  10,000  K,  then  ‘Blake’  will  be  chosen. 
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There  will  remain  the  lesser  task  of  providing  an  extended  temperature-range  library  for  use  with  the 
program. 

B.  Subtasks.  The  task  was  broken  down  into  the  following  subtasks: 

1.  Search  for  thermodynamic  data  at  temperatures  up  to  10,000  K  for  significant  species. 

2.  Devise  a  way  of  deciding  foe  extent  to  which  ions  are  formed  at  temperatures  below  15,000  K. 

3.  Based  on  the  these  results,  choose  either  CET89  or  ‘Blake’. 

a.  If  ions  are  significant,  then  choose  CET89  and  devise  a  means  of  estimating  covolumes. 

b.  If  ions  are  not  significant,  then  choose  ‘Blake*  and  fit  foe  thermodynamic  data  accordingly. 

C.  Thermodynamic  Data.  The  quest  for  thermodynamic  data  was  ended  even  before  it  started. 

1.  Bonnie  McBride,  National  Aeronautics  and  Space  Administration,  Lewis  Research  Center 
(NASA  LRC),  Cleveland,  Ohio,  kindly  furnished  to  the  U.S.  Army  Research  Laboratory  (ARL)  a  am  of 
coefficients  that  permit  computations  up  to  15,000  K  (and  20.000  K  in  some  cases)  for  about  a  hundred 
species. 

2.  The  format  of  these  tables  was  different  firm  that  used  in  all  of  foe  NASA  thermodynamic 
programs.  Hence  it  was  necessary  to  modify  one  of  them  to  accept  foe  new  format.  MuCET,  a 
modification  of  CET86  for  microcomputers,  was  selected  for  this  purpose.  The  new  program  was  named 
EXMuCET. 


IB.  PREPARING  EXMuCET 

A.  Description  of  the  Formats. 

1.  In  all  of  foe  NASA-Lewis  thermochemistry  programs  through  CET89,  the  heat  capacity,  Cp, 


enthalpy,  H,  and  entropy,  S,  were  represented  by  the  empirical  formulas 

Cp/R  *  a  +  b*T  +  c*TA2  +  d*TA3  +  e*TM  {1} 

H/R*T  *  a  +  b*T/2  +  c*T*2/3  +  d*T*3/4  +  e*TM/5  +  f/T  [2] 

S/R  «  a*lnCD  +  b *T  +  c*TA2/2  +  d*TA3/3  +  c*TM/4  +  g  {3} 

R  is  the  gas  constant,  so  Cp/R  and  S/R  are  dimensionless;  likewise,  H/R*T  is  dimensionless.  It  should 
be  noted  thin  Cp,  H,  and  S  are  functionally  related: 

H  =  JeprfT  +  constant;  and  (41 

S  *  j[Cp/T)dT  +  constant  (5) 


The  integration  constants  in  Eqns.  (4)  and  {5}  are  denoted  by  f  and  g  in  Eqns.  {2}  and  (3). 
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2.  The  fitting  range  is  split  from  300  to  1,000.  and  1X00  to  5000;  the  values  of  the  three  functions 
are  forced  to  match  at  the  boundary  (1,000  K).  This  makes  for  2  sett  of  a.  b.  ....  g.  (For  ballistic 
applications,  it  might  be  more  useful  to  make  the  fitting  ranges  1,000-2^00, 2.500-6.000.  but  this  point 
will  now  be  punued  at  this  time.) 

3.  In  theory  fids  same  format  could  be  used  for  higher  temperatures  (meaning  up  to  at  least  10X100  K, 
and  perhaps  tip  to  20,000  K).  Instead.  McBride  improved  the  accuracy  by  splitting  the  fitting  range  into 
3  parts,  with  boundaries  of  1XXX)  and  6,000.  This  format  will  be  called  here  foe  "intermediate  format." 

4.  Later,  she  adopted  a  different  format,  called  here  the  "new  format": 

Cp/R  •  afl*2  +  b/TAl  +  c  ♦  d*T  +  e*TA2  +  f*TA3  +  g*TM  {6} 

H/R*T  ■  -aXTA2  -Kb/T)*ln(T)  +  c  +  d*T/2  +  e*TA2/3  +  f*TA3/4 

♦  g*TN/5  +  h/T  |7) 

S/R  *  -a/(2*TA2)  -  b/T  +  c*lnCO  +  d*T  +  e*TA2/2  +  f*TA3/3 

+  g*TM/4  +  i  (8) 

H  and  S  are  related  to  Cp  just  as  they  were  before.  The  constants  a.  b, ....  g  are  of  course  different 
from  a.  b, ....  g  in  Eqns.  { 1 }— (3);  h  and  i  are  integration  constants.  Throughout  this  report,  it  is  this 
format  (Eqns.  {6}-{8})  that  will  be  referred  to  as  die  "new"  format 

5.  For  die  record  it  should  be  noted  that  McBride  included  provision  for  an  8th  term  (j*TA5)  in 
Eqn.  (6)  with  corresponding  terms  in  (7)  and  (8).  To  date  this  term  has  not  been  used. 

B.  Details  of  the  Changes. 

1.  Examination  of  CET86  showed  that  only  3  subroutines  had  to  be  modified:  UTHERM,  SEARCH, 
and  CPHS.  Technically.  Subroutine  HCALC  should  also  have  been  modified,  but  it  was  decided  not  to 
do  so.  This  subroutine  serves  to  compute  die  enthalpy  of  formation  of  some  (not  all)  input  species.  This 
feature  is  somewhat  useful,  but  its  function  is  largely  superseded  by  Program  PREP,  which  was  produced 
last  year  by  EFftA.  (See  the  documentation  of  PREP  for  further  details.) 

2.  The  programming  was  complicated  by  die  need  to  allow  for  die  future  addition  of  die  term  j*TA5 
in  Cp. 

3.  In  addition,  two  potential  difficulties  were  uncovered. 

a.  The  ranges  of  the  original  format  were  usually  298.15  (or  300)  K  to  1 XJOO,  and  1XXX)  to 
5,000  (or  6,000)  K,  at  least  for  gases.  For  liquids,  either  there  was  only  one  range,  300  K  to  boiling  point, 
or  else  the  two  ranges  were  300  to  1,000, 1 XXX)  to  boiling  point,  depending  on  whether  die  boiling  point 
was  below  or  above  1,000  K.  The  case  with  solids  was  similar,  either  300-melting  point,  or  300-1 XXX), 
1,000-melting  point  In  all  cases,  however,  the  actual  ranges  were  pan  of  the  input  for  each  species  and 
thus  were  available  to  the  program  throughout  its  execution.  Therefore,  when  computations  were  being 
ran  at  temperature  Tx,  the  program  could  decide  if  Tx  was  above  or  below  the  melting  (boiling)  point  of 
a  solid  (liquid),  and  take  appropriate  action. 
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b.  In  the  new  format,  the  ranges  are  fixed:  300-1,000,  1,000-6,000.  aid  6.000-10.000.  (Note 
that  all  three  ranges  are  not  necessarily  used  for  every  species;  the  actual  number  of  ranges  used  is  pat 
of  the  input.)  No  provision  has  been  made  for  inputting  melting  or  boiling  points.  In  theory  this  is  not 
necessary,  because  only  the  phase  with  the  most  negative  free  energy  car  form.  In  practice,  however, 
owing  to  artifacts  of  the  fitting  process,  it  is  posable  for  the  program  to  find 
thermodynamically-impossiNe  phases.  This  problem  has  been  side-stepped. 

c.  The  second  difficulty  was  similar  in  some  ways.  The  upper  ends  of  the  fitting  ranges  are  not 
the  same;  while  some  of  them  are  20,000  K,  most  of  them  are  only  5,000  or  6,000  K.  The  actual  value 
depends  on  the  available  thermochemical  data,  which  in  turn  reflea  physical  realities. 

d.  The  new  program  originally  included  a  tea  that  ensured  that  the  thermochemical  functions  for 
a  given  species  were  not  evaluated  outside  the  species’  fitting  range.  The  reason  for  including  this  tea 
was  indisputable'.  Owing  to  die  unknown  behavior  of  foe  fitted  thermodynamic  functions  outside  their 
range  of  validity,  it  is  possible  for  the  computation  to  come  up  with  a  species  under  conditions  where  it 
is  thermodynamically  impossible  for  it  to  exist.  Eventually  it  was  discovered  that,  in  foe  course  of 
converging,  the  program  occasionally  produced  intermediate  temperatures  that  were  outside  foe  fitting 
ranges  of  some  species.  (NASA’s  CET86  exhibits  similar  behavior.)  The  range  check  was  deleted  and 
the  program  worked  well. 

C.  Initial  Testing  Using  a  Modified  "OUT  Library. 

1.  Testing  was  accomplished  by  comparing  results  of  EXMuCET,  the  new  program,  with  those  from 
MuCET,  a  version  of  CET86  produced  by  EF&A  specifically  for  use  with  microcomputers.  It  was 
essential  that  this  testing  procedure  compare  only  programs,  not  thermodynamic  libraries.  For  this 
purpose,  2  additional  programs  were  written. 

a.  Program  OLDtNT  converted  the  standard  library  furnished  by  NASA’s  Lewis  Research  Center 
for  use  with  CET86  to  the  intermediate  format  described  in  ID  A  3  above. 

b.  Program  INTNEW  converted  the  library  in  foe  intermediate  format  to  foe  new  format  Note 
that  these  conversions  were  accomplished  strictly  by  reformatting  foe  libraries,  and  adding  zeros  for  foe 
coefficients  of  powers  of  T  that  were  not  present  in  foe  original  format.  No  refitting  was  done  by  either 
program. 


c.  These  two  programs  were  themselves  tested  by  inspection  of  the  results.  As  an  example  of 
this  transformation,  here  are  foe  entries  for  CH3: 


Original  Format: 


CH3 

0.28400327E  01 
0.16449813E  OS 
-0.18859236E-08 


J  6/69C  l.H  3.0  0.0  0.G 
0.60869086E-02-0.2 1 740338E-05 
0.5505675  IE  01  0.34666350E  01 
0.66803 182E- 12  0.16313104E  05 


300.000  5000.000  15.03470R 

0.36042576E-09 -0.22725 300E-13 
0.38301 845E-02  0.101 16802E-05 
0.24172192E  01  0.< 


I  M.I.M.M 


1 

2 

3 

4 
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New  Formic 


CH3 

2  J  6/69  C 
298.130 


1.00H  3.000  .000 

1000.000  7  -2.0  -1.0 

.346663500+01 
.668031820-12  .000000000+00 

1000.000  3000.000  7  -2.0  -1.0 


M  I  U  U  l)»ITU 


i  I  1 1  i  i  1 1  i 11 1  iWh  iXtYu;<i»TTU 


360425760-09  -22725300013 


.00  .00 
.0  1.0  2.0 
.38301843002 
.163131040+03 
JO  10  10 
284003270+01 
10+00 


III.IIJ  M 


0  13.03470  .OOOOOOOOE+OO 


MjTiTrriTTtT 


.10116802003  -.18839236008 

241721920+01 

3.0  4.0  .0 

.60869086002  -21740338003 
.164498130+05  23QS67S10+O1 


The  placement  of  the  high  temperature  data  before  the  low  temperature  ones  in  the  original  format  was 
one  of  the  little  details  that  complicated  an  otherwise  straightforward  programming  chore. 

d.  Using  these  new  programs,  die  standard  NASA  library,  FQJE40AT,  was  transformed  to  the 
new  format;  the  resulting  library  was  named  NEWFILE.DAT.  A  number  of  test  cases  were  run  with 
MuCET  and  EXMuCET  using  F1LE4DAT  or  NEWFILE.DAT,  as  appropriate. 

e.  Overall  the  comparison  was  excellent,  but  a  problem  caused  by  the  occurrence  (or  non¬ 
occurrence)  of  condensed  phases  was  revealed.  The  magnitude  of  this  problem  in  a  typical  case  can  be 
seen  by  comparing  the  output  from  3  runs. 

Run  A:  MuCET  with  FDLE4.DAT,  no  species  omitted. 

Run  B:  MuCET  with  FILE40AT,  all  condensed  species  omitted. 

Run  C:  EXMuCET  with  NEWFILE.DAT. 


RUN:  ABC 

THERMODYNAMIC  PROPERTIES 


P,  MPA 

21426 

214.49 

214.49 

T,  DEG  K 

2999.8 

2998.8 

2998.8 

RHO,  KG/CU  M 

200.0 

200.0 

200.00 

H.KJ/KG 

-50827 

-308.88 

-308.88 

U.KJ/KG 

•1381.33 

-1381.35 

-1581.35 

S,  KJ/(KGXK) 

9.4105 

9.4104 

9.4104 

CP,  KJ/(KGXK) 

1.9638 

1.9644 

1.9644 

CP  FROZN,  J/G-K 

1.8672 

1.8668 

1.8668 

EQUILIB  GAMMA 

12275 

12275 

12275 

FROZEN  GAMMA 

12369 

12370 

12370 

IMPETUS.  J/G 

1072.8 

10722 

10722 

B  ALLERGY*,  J/G 

45282 

4523.6 

4523.6 

*  B  allergy  *  Impetus/(gamma  -  1) 
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A  more  sensitive  comparison  is  ftimished  by  the  comparison  of  the  computed  mole  fractions: 


Run  A  B  C 


N2 

0.27696 

1 

N2 

0.27698 

1 

N2 

0.27698 

CO 

0.27467 

1 

CO 

0.27469 

1 

CO 

0.27469 

H20 

0.24136 

1 

H20 

0.24163 

1 

H20 

0.24163 

H2 

0.13166 

1 

H2 

0.13173 

1 

H2 

0.13173 

C02 

0.06932 

1 

C02 

0.06933 

1 

C02 

0.06933 

HF 

0.00179 

1 

HF 

0.00149 

1 

HF 

0.00149 

H 

0.00124 

1 

H 

0.00123 

1 

H 

0.00123 

OH 

0.00077 

1 

OH 

0.00077 

1 

OH 

0.00077 

NAOH 

0.00076 

1 

NAOH 

0.00077 

1 

NAOH 

0.00077 

NH3 

0.00038 

1 

NH3 

0.00038 

1 

NH3 

0.00038 

NA 

0.00016 

1 

ALF20 

0.00023 

1 

ALF20 

0.00023 

AL203(L) 

0.00011 

1 

NA 

0.00017 

1 

NA 

0.00017 

NO 

0.00011 

1 

NO 

0.00011 

1 

NO 

0.00011 

HCN 

0.00011 

1 

HCN 

0.00011 

1 

HCN 

0.00011 

ALF20 

0.00008 

1 

AL02H 

0.00007 

1 

AL02H 

0.00007 

HCORAD 

0.00006 

1 

HCORAD 

0.00006 

1 

HCORAD 

0.00006 

(HCOOH)2 

0.00005 

1 

(HCOOH)2 

0.000Q5 

1 

(HCOOH)2 

0.00003 

HCHO 

0.00004 

1 

HCHO 

0.00004 

1 

HCHO 

0.00004 

NAH 

0.00003 

1 

NAH 

0.00003 

1 

NAH 

0.00003 

NAF 

0.00003 

1 

NAF 

0.00003 

1 

NAF 

0.00003 

HNCO 

0.00003 

1 

HNCO 

0.00003 

1 

HNCO 

0.00003 

AL02H 

0.00002 

1 

NACN 

0.00001 

1 

NACN 

0.00001 

NACN 

0.00001 

1 

NH2 

0.00001 

1 

NH2 

0.00001 

NH2 

0.00001 

1 

CH4 

0.00001 

1 

CH4 

0.00001 

CH4 

0.00001 

1 

ALOH 

0.00001 

1 

ALOH 

0.00001 

1 

ALOF 

0.00001 

1 

ALOF 

0.00001 

f.  Evidently  the  effect  is  real  but  quite  small.  Since  ETC  propellants  do  not  fonn  any  condensed 
phases  at  chamber  conditions,  the  problem  is  not  serious. 


g.  Another  way  of  looking  at  dns  same  problem  is  to  consider  the  effect  of  adding  electrical 
energy  to  a  propellant  mix  that  forms  graphite  'jiolid  carbon)  in  the  absence  of  added  energy.  Such  a 
composition  is  a  stoichiometric  mixture  of  decalin  with  70%  hydrogen  peroxide. 
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Using  FILE4.DAT,  MuCET  gives 


DEC  ALIN  +  709b  HYDROGEN  PEROXIDE  +  VARYING  ENERGY 


CHEMICAL  FORMULA 

WTFRAC 

ENERGY 

KJ/KG-MOL 

S 

T 

(K) 

FUEL  C  10.000  H  18.000 

1.00000 

-230538. 

L 

298.1 

OXIDH  2.000  0  2.000 

.70000 

-187778. 

L 

298.1 

OXIDH  2.000  0  1.000 

.30000 

-285830. 

L 

298.1 

THERMODYNAMIC  PROPERTIES 

Elec  Nigy,  J/G 

0 

1000. 

2000. 

U.J/G 

-5145.85 

-4145.85 

-3145.85 

P,  MPA 

99.797 

123.77 

150.32 

T,  DEG  K 

1271.7 

1431.0 

1586.9 

CP  EQ,  J/G-K 

7.1134 

8.2267 

7.9425 

CP  FROZ.  J/G-K 

2.7930 

2.8810 

2.9433 

EQUHJB  GAMMA 

1.1326 

1.1427 

1.1589 

FROZEN  GAMMA 

1.1634 

1.1766 

1.1918 

IMPETUS,  J/G 

499.0 

618.9 

751.6 

BALLERGY,  J/G 

3052.9 

3503.7 

3919.1 

MOLE  FRACTIONS 


HCHO 

0.00002 

1 

HCHO 

0.00005 

1 

HCHO 

0.00010 

(HCOOH)2 

0.00002 

1 

9 

1 

T5 

0.00002 

1 

(HCOOH)2 

0.00003 

CH4 

0.26007 

1 

CH4 

0.24128 

1 

CH3 

0.00001 

CH30H 

0.00001 

1 

CH30H 

0.00001 

1 

CH4 

0.22032 

CO 

0.10092 

1 

CO 

0.19191 

1 

CH30H 

0.00002 

C02 

0.09098 

1 

C02 

0.06541 

1 

CO 

0.27833 

C2H4 

0.00004 

1 

C2H4 

0.00015 

1 

C02 

0.04009 

ETHANE 

0.00052 

1 

CH3CHO 

0.00001 

1 

C2H4 

0.00001 

H2 

0.14913 

I 

ETHANE 

0.00068 

1 

KETENE 

0.00001 

ICO 

0.22097 

1 

n-C3H6 

0.00001 

1 

C2H4 

0.00039 

C(GR) 

0.17731 

1 

PROPANE 

0.00001 

1 

CH3CHO 

0.00001 

1 

H2 

0.21653 

1 

ETHANE 

0.00081 

I 

ICO 

0.16898 

1 

n-C3H6 

0.00002 

1 

C(GR) 

0.11495 

1 

PROPANE 

0.00001 

1 

H2 

0.28297 

1 

H2 

0.34901 

1 

ICO 

0.11972 

1 

ICO 

0.07599 

1 

C(GR) 

0.05715 

1 

C(GR) 

0.01526 
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THERMODYNAMIC  PROPERTIES 


Elec  Nigy,  J/G 

3000.  1 

4000. 

1 

5000. 

U.J/G 

-2145.85  1 

-1145.85 

1 

-145.85 

P,  MPA 

180.97  1 

217.04 

1 

257.55 

T.DEGK 

17632  1 

1978.7 

1 

2234.5 

CP  EQ.  J/G-K 

6.7622  1 

5.5394 

1 

4.7839 

CPFROZ.  J/G-K 

2.9984  1 

3.0536 

1 

3.1056 

EQUILIB  GAMMA  1.1789  1 

12074 

1 

12124 

FROZEN  GAMMA  1.2065  1 

12189 

1 

12278 

IMPETUS.  J/G 

904.8  1 

10852 

1 

1287.8 

BALLERGY,  J/G 

4381.9  1 

4957.0 

1 

5651.9 

MOLE  FRACTIONS 


HCHO 

0.00015 

1 

HCHO 

0.00022 

1 

HCHO 

(HCOOH)2 

0.00002 

1 

(HCOOH)2 

0.00001 

I 

(HOOOH)2 

CH3 

0.00005 

1 

CH3 

0.00017 

1 

CH3 

CH4 

0.19312 

1 

CH4 

0.15654 

1 

CH4 

CH30H 

0.00002 

1 

CH30H 

0.00002 

1 

CH30H 

CO 

0.34392 

1 

CO 

0.37579 

1 

CO 

C02 

0.02050 

1 

C02 

0.00900 

1 

C02 

C2H4 

0.00005 

1 

C2H4 

0.00030 

1 

C2H4 

KETENE 

0.00003 

1 

KETENE 

0.00007 

1 

KETENE 

C2H4 

0.00089 

1 

C2H4 

0.00177 

1 

C2H3RAD 

CH3CHO 

0.00002 

1 

CH3CHO 

0.00002 

1 

C2H4 

ETHANE 

0.00090 

1 

C2H5 

0.00001 

1 

CH3CHO 

n-C3H6 

0.00004 

1 

ETHANE 

0.00068 

1 

C2H5 

PROPANE 

0.00001 

1 

PROPYNE 

0.00001 

1 

ETHANE 

ALLENE 

0.00001 

1 

C3H3RAD 

0.00001 

1 

C3H5RAD 

0.00001 

1 

PROPYNE 

0.00006 

1 

n-C3H6 

0.00009 

1 

ALLENE 

0.00005 

1 

Using  NEWFDLE.DAT,  EXMUCET  gives 

DECAUN  +  70%  HYDROGEN  PEROXIDE  +  VARYING  ENERGY 


CHEMICAL  FORMULA 

FUEL  C  10.000  H  18.000 
OXIDH  2.000  0  2.000 
OXIDH  2.000  0  1.000 


WTFRAC 

ENERGY 

J/KG-MOL 

S 

T 

(K) 

1.00000 

•230538. 

L 

298.1 

.70000 

187778. 

L 

298.1 

.30000 

285830. 

L 

298.1 

0.00028 

0.00001 

0.00053 

0.12068 

0.00002 

0.38422 

0.00385 

0.00144 

0.00013 

0.00002 

0.00313 

0.00003 

0.00004 

0.00078 
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THERMODYNAMIC  PROPERTIES 


ELEC  NRGY,  J/G 

0.0  I 

1000. 

1 

2000. 

U.J/G 

-5145.85  1 

-4145.85 

1 

-3145.85 

P,  MPA 

94.913  1 

119.50 

1 

147.73 

T.DEGK 

1220.6  1 

1392.1 

1 

1567.4 

CP  EQ,  J/G-K 

6.8719  1 

7.3472 

1 

7.1432 

CP  FROZ,  J/G-K 

2.7409  1 

2.8463 

1 

2.9285 

EQUILIB  GAMMA 

1.1364  1 

1.1534 

1 

1.1730 

FROZEN  GAMMA 

1.1653  I 

1.1776 

1 

1.1918 

IMPETUS,  J/G 

474.6  1 

597.5 

1 

738.6 

BALLERGY.  J/G 

2870.8  1 

3364.9 

1 

3851.5 

MOLE  FRACTIONS 


HCHO 

0.00002 

1 

HCHO 

0.00005 

1 

HCHO 

0.00010 

(HCOOH)2 

0.00002 

1 

(HCOOH)2 

0.00002 

1 

(HCOOH)2  0.00002 

CH4 

0.43043 

1 

CH4 

0.34248 

1 

CH3 

0.00001 

CH30H 

0.00001 

1 

CH30H 

0.00001 

1 

CH4 

0.26689 

CO 

0.18553 

1 

CO 

0.26984 

1 

CH30H 

0.00002 

C02 

0.15400 

1 

C02 

0.08467 

1 

CO 

0.32442 

C2H4 

0.00013 

1 

KETENE 

0.00001 

1 

C02 

0.04286 

CH3CHO 

0.00001 

1 

C2H4 

0.00030 

1 

C2H2 

0.00001 

ETHANE 

0.00151 

1 

CH3CHO 

0.00001 

1 

KETENE 

0.00001 

n-C3H6 

0.00001 

1 

ETHANE 

0.00138 

1 

C2H4 

0.00058 

PROPANE 

0.00002 

1 

n-C3H6 

0.00002 

1 

CH3CHO 

0.00002 

H2 

0.10379 

1 

PROPANE 

0.00002 

1 

ETHANE 

0.00119 

H20 

0.12453 

1 

H2 

0.18055 

1 

n-C3H6 

0.00003 

• 

1 

H20 

0.12064 

1 

PROPANE 

0.00001 

1 

H2 

0.26417 

1 

1120 

0.09966 

THERMODYNAMIC  PROPERTIES 

ELEC  NRGY,  J/G 

3000.  1 

4000. 

1 

5000. 

U.J/G 

-2145.85  1 

•1145.85 

1 

-145.85 

P,  MPA 

180.16  I 

217.04 

1 

257.55 

T.DEGK 

1758.8  1 

1978.7 

1 

2234.5 

CP  EQ,  J/G-K 

6.4774  1 

5.5394 

1 

4.7839 

CP  FROZ,  J/G-K 

2.9959  1 

3.0536 

1 

3.1056 

EQUILIB  GAMMA 

1.1922  1 

1.2074 

1 

1.2124 

FROZEN  GAMMA 

1.2062  1 

1.2189 

1 

1.2278 

IMPETUS,  J/G 

900.8  1 

10852 

1 

1287.8 

BALLERGY,  J/G 

4368.4  1 

4957.0 

1 

5651.9 
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MOLE  FRACTIONS 


HCHO 

0.00016 

1 

HCHO 

0.00022 

1 

HCHO 

0.00028 

(HCOOH)2 

0.00002 

1 

(HCOOH)2 

0.00001 

1 

(HCOOH)2 

0.00001 

CH3 

0.00005 

1 

CH3 

0.00017 

1 

CH3 

0.00053 

CH4 

0.20485 

1 

CH4 

0.15654 

1 

CH4 

0.12068 

CH30H 

0.00002 

1 

CH30H 

0.00002 

1 

CH30H 

0.00002 

CO 

0.35719 

1 

CO 

0.37579 

1 

CO 

0.38422 

C02 

0.02032 

1 

C02 

0.00900 

1 

C02 

0.00385 

C2H2 

0.00006 

1 

C2H2 

0.00030 

1 

C2H2 

0.00144 

KETENE 

0.00003 

1 

KETENE 

0.00007 

1 

KETENE 

0.00013 

C2H4 

0.00101 

1 

C2H4 

0.00177 

1 

C2H3RAD 

0.00002 

CH3CHO 

0.00002 

1 

CH3CH0 

0.00002 

1 

C2H4 

0.00313 

ETHANE 

0.00101 

1 

C2H5 

0.00001 

1 

CH3CHO 

0.00003 

n-C3H6 

0.00005 

1 

ETHANE 

0.00088 

1 

C2H5 

0.00004 

PROPANE 

0.00001 

1 

PROPYNE 

0.00001 

1 

ETHANE 

0.00078 

H2 

0.34404 

1 

ALLENE 

0.00001 

1 

C3H3RAD 

0.00001 

H20 

0.07115 

1 

C3H5RAD 

0.00001 

1 

PROPYNE 

0.00006 

1 

n-C3H6 

0.00009 

1 

ALLENE 

0.00005 

1 

PROPANE 

0.00001 

1 

C3H5RAD 

0.00003 

1 

H 

0.00002 

1 

n-C3H6 

0.00015 

1 

HCORAD 

0.00001 

1 

PROPANE 

0.00001 

1 

H2 

0.41098 

1 

1.3-C4H6 

0.00001 

1 

H20 

0.04406 

1 

H 

0.00009 

1 

HCORAD 

0.00002 

1 

H2 

0.45978 

1 

H20 

0.02462 

At  temperatures  of  1,587  K  and  below,  MuCET  shows  C(gr)  in  the  products;  there  are  moderate 
differences  in  the  thermodynamic  properties  computed  by  the  two  programs  in  the  same  region.  This 
difference  decreases  at  1,759  K  and  vanishes  at  1,979  K.  This  point  could  have  been  made  even  stronger 
by  adding  energy  in  increments  smaller  than  1,000  J/g.  Since  EXMuCET  will  be  used  primarily  for 
temperatures  above  5,000  K,  this  problem  with  the  omission  of  condensed  phases  is  not  important  and  will 
be  ignored. 


D.  Initial  Testing  Using  the  Extended  "New"  Library. 


1.  Early  tests  with  EXMuCET  using  McBride’s  extended  library  gave  terrible  results. 

a.  The  most  serious  problem  was  that  foe  new  (extended  range)  library  gave  flame  temperatures 
around  9,300  K! 
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b.  Furthermore,  many  of  the  test  casts  would  not  even  run  with  the  extended  library.  The  output 
gave  the  uninformative  message,  SINGULAR  MATRIX.  This  was  eventually  traced  to  the  lack  of  data 
for  compounds  of  potassium,  sodium,  fluorine,  aluminum,  and  magnesium  in  the  new  library.  These 
compounds  are  used  as  flash  suppressants  or  stabilizers  in  many  propellants.  Accordingly  all  of  them  were 
removed  from  all  of  the  compositions;  the  new  compositions  are  labelled  "pseudo-. . . "  to  emphasize  this 
change. 

2.  The  previous  testing  of  EXMuCET  using  the  transformed  standard  library  had  given  no  indications 
of  trouble  (see  '  jn  C,  above).  Based  on  those  results,  we  looked  elsewhere  for  die  problem.  The  first 
hypothesis  w&>  .at  there  was  a  problem  with  die  database.  In  view  of  the  well-deserved  high  reputation 
of  NASA’s  Lewis  Research  Center,  this  hypothesis  seemed  very  unlikely,  but  it  could  not  be  ruled  out 
a  priori.  As  will  be  seen  below,  this  hypothesis  was  false. 

3.  Testing  the  data  base  required  writing  a  new  program,  CALCTH.FOR,  whose  purpose  was  to 
generate  tables  of  the  thermodynamic  functions.  The  results  for  N  atoms  were  obviously  incorrect  This 
problem  was  readily  corrected.  It  turned  out  that  there  was  a  similar  problem  in  Subroutine  EXUTHERM 
in  EXMuCET.  which  was  also  corrected. 

E.  Testing  EXMuCET  with  Different  Libraries.  It  was  shown  above  that  EXMuCET  gives  the  same 
results  with  McBride’s  standard  library  transformed  to  the  new  format  as  MuCET  does  with  die  standard 
untransformed  library.  The  more  interesting  question  remained.  Would  EXMuCET  give  die  same  results 
at  temperatures  below  5,000  K  with  the  extended  library  as  MuCET  does  with  die  standard  library? 

1.  To  this  end,  a  DOS  batch  file,  TESTLIBS.BAT,  was  written  (sec  listing  in  Appendix  A).  This  file 
first  creates  two  binary  libraries  for  use  with  EXMuCET.  The  first  library  was  taken  from  the  NASA 
library  distributed  with  CET89  (suitably  transformed  for  use  with  EXMuCET — see  above  for  details).  The 
second  library  was  McBride’s  extended  library.  The  two  binary  libraries  were  saved  for  later  use.  The 
same  test  cases  were  run  with  each  library  and  the  results  compared. 

2.  It  was  deemed  necessary  not  only  to  run  a  number  of  different  test  cases,  but  also  to  run  each  of 
them  with  a  series  of  added  electrical  energies  so  as  to  be  sure  of  covering  die  temperature  range  up  to 
5,000  K  and  even  a  little  beyond.  Generating  the  input  for  such  tests  was  too  tedious  to  be  undertaken 
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by  hand,  so  a  QuickBasic  program,  MAKEINP.BAS,  was  written  that  automttcd  most  of  the  work.  bis 
lined  in  Appendix  B. 

a.  The  user  inputs  the  desired  starting  energy,  the  energy  increment,  tod  the  final  energy  (all  in 
Joules/g).  The  program  reads  a  template  containing  the  input  for  8  propellants,  including  the  NAMELIST, 
and  then  produces  a  properly-formatted  input  file  for  EXMuCET  for  all  the  compositions  at  each  energy; 
next,  it  prodr  "es  a  second,  almoa-identical  file,  except  that  in  this  case,  the  NAMELIST  contains  the 
option,  IONS=T. 

b.  The  eight  propellants  or  propellant  systems  used  were: 

1.  Pseudo  M30  4.  M9  7.  Decalin  ♦  70%  HjQ* 

2.  Pseudo  WC890  5.  Ml  8.  RFNA*  ♦  UDMH* 

3.  JA-2  6.  Decalin +WFNA* 

*  RFNA  «=  red  fuming  nitric  add.  UDMH  *  unsymmetrical  dimethyl  hydrazine.  WFNA  »  white 
fuming  nitric  acid. 

Details  of  these  compositions  are  given  in  Appendix  C,  where  the  template  is  primed. 

c.  The  title  of  each  case  lists  the  added  energy  as  "#  Joules/g”;  tbe  program  changes  the  "#"  to 
amount  of  added  electricity  chosen  by  the  user.  Tbe  NAMELIST  contains  die  option,  here,  die 
program  properly  changes  the  "#"10  the  correct  total  energy  of  the  system,  which  includes  both  the  initial 
energy  of  die  composition  as  well  as  tbe  added  electrical  energy.  For  simplicity  all  of  the  energies  are 
in  Joules. 

3.  The  SUMMARY  files  for  different  outputs  produced  by  EXMuCET  are  then  compared  by  a 
QuickBasic  program,  COMPARED  AS,  which  is  listed  in  Appendix  D. 

The  various  files  differ  from  each  other  according  to  which  thermodynamics  library  was  used  to 
produce  them,  and  whether  or  not  ions  were  permitted  to  form. 
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a.  COMPARE  identifies  the  thermodynamics  libraries  being  used  and  the  names  of  the  input  files, 
which  arc  printed  for  positive  identification  of  the  output 

b.  It  then  proceeds  to  compare  temperatures,  impetus,  gammas,  and  ballistic  energies  (labelled 
"ballergy"  in  the  output)  for  any  two  specified  files;  it  also  computes  various  measures  of  agreement 

c.  The  results  obtained  after  the  final  corrections  to  the  program  were  better  overall,  but  sdll 
showed  problems,  at  least  at  the  lower  temperatures.  Examination  of  the  file  OUTPUT  revealed  the 
occurrence  of  the  message  "ADD  H20  (L)*  even  though  the  approximate  flame  temperature  was  3,200  K! 
The  upper  range  for  die  thennodynamics  data  for  HjO  (Hq)  in  the  library  is  400-600  K;  owing  to  an 
apparent  defect  in  the  programming,  this  fact  is  not  being  used  by  the  program  to  eliminate  HgO  (lag)  from 
consideration.  This  defect  remains  elusive  in  spite  of  additional  effort  Pending  its  removal,  all  condensed 
phases  (including  graphite)  have  been  eliminated  from  the  database.  This  is  admittedly  inelegant;  on  the 
other  hand,  EXMuCET  is  intended  for  use  at  temperatures  above  5.000  K,  where  no  condensed  phases 
exist 


4.  The  final  results  of  die  various  comparisons  are  given  in  Appendixes  E,  F,  and  G.  The 
nomenclature  for  the  various  files  is  as  follows.  The  root  of  all  names  is  TSTLB.  To  this  is  appended 
in  order 


S  or  O  for  SUMMARY  or  OUTPUT,  respectively,  and 
NO  or  10  for  NO  [ions]  or  10  [ions  included],  respectively. 

5.  There  are  only  small  differences  in  die  results  produced  using  the  standard  (CET89)  library 
compared  to  those  produced  using  the  extended  library. 

a.  Further  evidence  that  the  two  libraries  are  in  virtually  complete  agreement  with  each  other  is 
shown  by  a  comparison  of  die  amounts  of  products  produced  by  each.  Appendix  H  gives  such  a 
comparison  for  one  case,  composition  1  (pseudo  M30)  at  loading  density  *  0.2  gfcc  and  with  5,000 1 fg 
of  added  energy.  This  is  an  extremely  sensitive  test,  so  the  agreement  is  quite  satisfying. 


39 


library  is  used,  die  omtafon  or  inclusion  of  km  has  no  effect  Therefore  the  program  of  choice  for  use 
with  ETC  propellants  it  ‘Blake'  with  a  suitably  prepared  library. 

F.  Intercomparisom  of  Various  Thermodynamic  Tobies  and  Coefficients. 

1.  The  principal  databaac  on  which  the  present  work  is  baaed  is  the  set  of  codBcfcnta  for 
approximately  133  C-H-O-N  specks  obtained  from  McBride  (NASA  LRQ.  The  vabie  of  this  set  k>  the 
present  work  cannot  be  overestimated;  nothing  could  have  been  done  without  iL  Nevertheless,  when 
tables  of  thermodynamic  data  for  28  C-H-O-N  species  were  belatedly  seat  to  ARL  by  die  National 
Standard  Data  Reference  System  of  the  National  Institute  for  Standards  and  Technology  (NIST)  (formerly 
the  National  Bureau  of  Standards),  it  appeared  prudent  to  compare  them  with  tables  produced  using  the 
LRC  coefficients.  These  species  are  listed  in  Appendix  I. 

To  this  end  two  new  programs  were  written. 

a.  XCALCTH4  is  a  variant  of  the  program  CALCTH  previously  mentioned  in  Sect  IILD.3  above. 

i.  The  input  to  XCALCTH4  was  an  edited  subset  of  the  LRC  coefficients.  This  editing 
reduced  it  to  the  28  species  in  die  NIST  tables  and  converted  their  names  so  that  they  were  identical  in 
both  sets.  This  was  necessary  to  ensure  that  identical  species  were  being  compared. 

ii.  The  output  from  XCALCTH4  was  a  table  of  the  temperature.  T.  and  die  three  principal 
thermodynamic  functions:  the  heat  capacity.  Cp(T);  the  entropy,  S(T);  and  the  enthalpy  of  formation  at 
298  K  plus  the  increment  in  enthalpy  from  298  K  to  T,  Delta  Hf(298)  +  H(t)  -  H(298).  Each  of  these 
entries  was  set  off  by  commas  to  make  them  suitable  for  input  to  die  next  program.  COMPARE3. 

b.  COMPARE3  is  a  BASIC  program  that  compares  the  two  sets  of  tables  at  every  temperature 
for  each  species.  It  was  programmed  to  take  into  account  the  fact  that  the  increments  in  temperatures  in 
the  two  tables  might  be  different  (in  which  case  only  die  values  at  the  larger  increments  were  used);  and 
the  fact  that  one  table  might  end  before  the  other.  This  was  because  some  of  die  species  in  the  LRC  set 
were  fitted  only  to  6,000  K,  while  the  NIST  tables  went  uniformly  to  10400  K. 
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c.  The  results  of  the  comparison  are  given  in  Appendix  J.  which  gives  the  mean  ahaotot 
deviation  and  the  rooc-mean-square  percent  deviations  for  the  three  principal  thermodynamic  functions. 
The  table  has  been  sorted  in  order  of  Increasing  RMS  percentage  error  in  the  enthalpy. 

i.  Overall  the  comparison  is  satisfactory:  Only  three  species  have  RMS  errors  in  their 
enthalpies  exceeding  5%,  and  they  are  all  ions,  which  will  be  of  little  or  no  importance  in  determining 
die  range  of  validity  of  EXMuCET. 

ii.  There  is  no  way  of  deciding  from  this  comparison  alone  which  set  of  data  is  die  belter. 
The  important  point  is  that  they  are  consistent 

IV.  ESTABLISHING  THE  NEW  LIBRARY 

A.  The  next  step  was  the  determination  of  die  29  species  that  would  constitute  the  new  ‘Blake’  library. 
This  limit  is  imposed  by  die  program  itself:  it  could  be  raised  to  about  35  without  too  much  extra  work, 
but  there  does  not  appear  to  be  any  compelling  reason  for  this  change. 

1.  To  this  end.  8  propellant  systems,  listed  in  Table  1,  were  selected. 

Table  1.  Systems  Studied  for  Determining  Important  Species 

1.  Pseudo  M30  5.  Ml 

2.  Pseudo  WC890  6.  Decalin  +  WFNA 

3.  JA-2  7.  Decalin  +  70%  HjOj 

4.  M9  8.  RFNA  +  UDMH 

These  systems  were  then  evaluated  with  EXMuCET  for  a  wide  range  of  added  energies  (horn  0  to  as 
much  as  40jOOOJ/g);  the  highest  resulting  temperature  was  12,924  K.  The  output  fiom  these  calculations 
was  then  sorted  in  order  of  decreasing  mole  fraction  for  each  species,  and  then  all  entries  with  mole 
fractions  less  than  1  x  10'3  were  deleted.  The  results  are  shown  in  Table  2. 

The  column  labelled  "Maximum"  is  self-explanatory.  The  column  labelled  "Minimum"  lists,  for  a  given 
species,  the  smallest  computed  mole  fraction  greater  than  1  x  10*3  for  that  species,  and  the  corresponding 
temperature. 
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MINIMUM* 


Table  2.  Oonosiinifcni  of  Various  Spades  ftoduce 
with  Varying  Added  Eim^kf 


MAXIMUM 


Mole 

Mole 

£io& 

£oc 

Temp  (K) 

Ex 

TempflC) 

1. 

Ha 

03485 

3675. 

0.0173 

12924. 

2. 

CO 

03062 

2491. 

0.0764 

11793. 

3. 

CH+ 

0.4312 

1217. 

0.0014 

8017. 

4. 

HjO 

0.4278 

3739. 

0.0022 

12924. 

5. 

H 

0.3543 

11793. 

0.0011 

2950. 

6. 

0 

02857 

12924. 

aooio 

3812. 

7. 

Ha 

02774 

3010. 

0.0129 

12398. 

8. 

COj 

0.2068 

3812. 

0.0010 

8017. 

9. 

OH 

01446 

6281. 

0.0031 

6145. 

10. 

N 

0.0926 

12383. 

0.0017 

6281. 

11. 

C 

0.0801 

12416. 

0.0011 

7990. 

12. 

NO 

0.0519 

7517. 

0.0016 

3575. 

13. 

Oa 

0.0400 

6281. 

0.0013 

5091. 

14. 

Acetylene 

0.0331 

5041. 

0.0012 

2233. 

15. 

CN 

0.0305 

12383. 

aooi5 

8465. 

16. 

NH 

00223 

11793. 

0.0013 

6281. 

17. 

CH 

0.0189 

8882. 

0.0016 

8162. 

18. 

qn 

0.0182 

8882. 

0.0010 

11424. 

19. 

HCO 

0.0109 

8882. 

0.0010 

5388. 

20. 

cq 

0.0104 

8882. 

0.0011 

9532. 

21. 

CHj 

0.0064 

5041. 

0.0033 

8882. 

22. 

Ca 

0.0064 

12416. 

0.0012 

10977. 

23. 

NHj 

0.0044 

11793. 

0.0010 

6776. 

24. 

HCN 

0.0037 

9532. 

aooio 

6599. 

25. 

Ethylene 

0.0034 

3675. 

0.0010 

5041. 

26. 

CCO  (radical) 

0.0029 

12416. 

0.0012 

10977. 

27. 

Os 

O.OQ23 

12924. 

0.0013 

12647. 

28. 

Hq 

0.0022 

6281. 

0.0010 

9604. 

29. 

Vinylidene# 

0.0022 

6145. 

0.0014 

8882. 

niTJWB  fnr  •BUfr.' 

W  W  M 

~wa  m  ivi 

UNO 

0.0019 

8465. 

0.0010 

10409. 

[ELECTRONS 

0.0018 

12924. 

0.0010 

12035.] 

NCN 

0.0018 

11456. 

0.0012 

9280. 

ETHANE 

0.0015 

1217. 

0.0015 

1217. 

NH3 

0.0014 

11793. 

0.0011 

12454. 

•  "Minimum”  here  reft  is  to  the  smallest  mole  fraction  greater  dun  1  x  10"3. 

•  Hus  is  the  radical  HjOC 
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The  significance  of  Table  2  can  be  made  dear  by  examining  some  of  the  eariy  entries  in  more  detail. 
Hydrogen  (Hj)  is  an  impoitam  product  whose  concentration  readies  a  maximum  (0.55)  at  3,675  K  and 
then  declines  to  0.017  at  12,294  K,  which  is  the  highest  temperature  studied.  Carbon  monoxide  (CO) 
behaves  similarly.  Methane  (CH^  reaches  its  maximum  concentration  at  a  lower  temperature,  and  then 
declines  until  it  becomes  0.001  at  8,017  K;  its  mole  fraction  at  all  higher  temperatures  is  less  than 
0.00100.  Water  (HjO)  behaves  similarly  to  Hj.  Hydrogen  atom  (H)  is  different:  Its  mole  fraction  is  less 
than  0.00100  at  temperatures  below  2,950  K;  its  concentration  increases  with  increasing  temperature  until 
it  reaches  a  maximum  at  11,793  K.  Oxygen  atom  (O)  behaves  similarly. 

In  its  present  form,  'Blake'  is  limited  to  working  with  29  gaseous  species  at  any  one  time.  Therefore, 
based  on  the  results  in  Table  2,  the  29  spedes  that  had  maximum  concentrations  equal  to  or  greater  than 
0.0022  anywhere  in  the  temperature  range  1,217  <  T  <  12,924  K  were  chosen  to  be  species  in  the 
extended  'Blake*  library. 

2.  Thermodynamic  data  for  the  five  permanent  gases,  N2.  CO.  CO2.  Hj,  and  HjO  were  taken  from 
McBride’s  tables. 

3.  Thermodynamic  data  for  all  of  the  other  species  were  obtained  by  expansion  of  die  coefficients 
in  the  extended  library  that  McBride  had  earlier  supplied  to  the  ARL. 

4.  The  fitting  of  these  data  to  the  form  used  by  'Blake*  was  accomplished  using  variants  of 
W.  H.  Zwisler’s  program,  STARFIT,  to  which  various  measures  of  the  goodness  of  fit  had  been  added. 

5.  The  thermodynamic  data  required  for  the  fittings  are  the  heat  capacity  as  a  function  of  temperature 
over  the  desired  range,  and  the  enthalpy  of  formation  and  the  entropy  at  some  reference  temperature. 

a.  In  theory  the  "desired  temperature  range"  could  be  chosen  to  extend  from  300  K  to  10,000  K, 
but  this  choice  would  not  be  suitable  for  'Blake*  in  general  because  the  available  thermodynamic  data  are 
limited  to  compounds  of  C,  H,  O,  and  N.  It  was  deemed  more  useful  to  retain  the  original  (old)  ‘Blake* 
library  for  almost  all  calculations,  and  to  obtain  an  extended  (new)  library  only  for  use  with 
ETC  formulations  (which  usually  do  not  contain  inorganic  additives). 
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b.  In  the  initial  attempts  at  fitting.  298.15  K  was  laid  wired"  into  the  program  as  the  reference 
temperature;  the  staffing  temperature  was  varied  from  300  K  to  5.000  K.  Some  fittings  were  made  with 
a  higher  reference  temperature.  They  showed  that  the  RMS  error  of  the  fit  decreased  as  the  reference 
temperature  increased,  so  the  fitting  programs  were  modified  to  make  both  the  reference  temperature  and 
the  starting  temperature  inputs  to  the  program.  Survey  runs  showed  that  die  RMS  error  was  least  when 
these  two  temperatures  were  the  same. 

c.  The  fust  of  these  modified  programs,  TIGFIT5,  which  is  listed  in  Appendix  K,  was  used  to 
fit  the  thermodynamic  data  for  the  five  permanent  gases.  The  second  program,  TIGF1T5X.  listed  in 
Appenlix  L,  was  used  for  fitting  all  of  the  other  species.  Two  programs  were  used  to  accommodate  die 
different  input  formats.  It  would  not  have  been  difficult  to  adapt  one  of  the  programs  to  accept  different 
input  formats,  but  it  was  simpler  and  faster  to  use  different  programs. 

6.  Certain  adjustments  had  to  be  made  in  the  input  files  in  order  to  make  the  comparisons  meaningful. 
Even  when  limited  to  C-H-N-O  species,  there  is  not  a  one-to-one  correspondence  between  die  two 
libraries.  All  in  all.  there  are  17  species  that  appear  only  in  one  or  the  other  of  the  libraries  but  not  in 
both.  REJect  instructions  for  these  17  species  were  inserted  into  all  runs  in  order  to  ensure  that  the  same 
species  appeared  regardless  of  which  library  was  used. 

7.  Seven  compositions  were  selected  for  the  comparison: 

1)  Pseudo  M30*  5)  Ml 

2)  JA-2  6)  Decalin  +  WFNA 

3)  Pseudo  M5  7)  Decalin  +  70%  H202 

4)  Pseudo  M9 

*  Pseudo  here  and  elsewhere  means  that  all  inorganic  additives  were  omitted  from  die 

formulation. 

8.  The  next  matter  to  be  settled  was  the  starting  temperature  (the  final  temperature  would  always  be 
10,000  K). 

a.  The  standard  (old)  ‘Blake’  fits  cover  the  range  500  K-5,000  K,  and  are  probably  valid 
somewhat  beyond  that  range,  say  to  5,500  K.  Thus  it  would  be  possible  to  have  the  extended  (new) 
library  start  at  5,000  K.  Nevertheless  it  seemed  preferable  to  have  it  start  at  a  lower  temperature. 
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b.  A  further  requirement  was  established,  that  there  be  little  or  no  difference  in  results  obtained 
with  either  library  in  their  overlapping  temperature  regions.  The  computed  quantities  concerned  were 
adiabatic  flame  temperature,  impetus,  and  ballergy,  which  is  the  ratio  Impetus/(gamma  •  1).  Ideally  it 
would  have  been  desirable  to  make  such  comparisons  for  extended  libraries  with  starting  temperatures 
from  2.000  K  to  5,000  K,  but  this  was  not  feasible.  Two  starting  temperatures  were  selected,  3  XXX)  and 
4,000  K.  The  corresponding  libraries  were  named  30BLKLYBUB  and  40BLKLYBUB,  which  were 
used  to  form  the  binary  libraries  30BLKLYB.DAT  and  40BLKLYB.DAT 

9.  The  comparisons  were  made  by  plotting  the  percent  differences  in  the  three  quantities  against  the 
temperature  computed  by  die  extended  (new)  library.  The  resulting  plots  for  the  two  libraries  were  then 

compared  by  eye.  These  plots  are  shown  in  Fig-l-A,  Fig-l-B . Rg-l-G.  It  is  readily  seen  that  the 

better  match  is  found  between  the  standard  (old)  library  and  the  extended  library  that  starts  at  3,000  K_ 
This  library  was  therefore  chosen  as  the  final  outcome  of  the  present  task.  The  names  of  die  chosen 
alphanumeric  and  binary  libraries  were  changed  to  XBLKLYB.LIB  and  XBYNLYB.DAT,  respectively. 

10.  As  a  matter  of  interest  the  effect  of  the  REJect  instructions  was  examined  by  comparing  die 
results  with  and  without  them.  The  resulting  (dots  are  shown  in  Fig-2-A,  Fig-2-B,  ....  Rg-2-G.  A 
noticeable  difference  is  found  rally  in  the  last  case.  It  is  caused  by  the  appearance  of  solid  carbon  when 
it  is  permitted. 


V.  NBLAKE 


A.  Although  it  was  not  required  in  the  Statement  of  Work  for  this  task,  it  seemed  appropriate  to  make 
some  changes  and  improvements  in  ’Blake’,  but  the  Took  and  feeT  of  the  program  have  not  been 
changed. 

B.  A  detailed  description  of  the  changes  is  given  in  the  documentation  for  the  revised  program.  The 
principal  changes  concern  the  library;  die  net  result  is  dial  the  new  library  is  incompatible  with  all 
previous  versions  of  die  program,  and  previous  versions  of  the  program  are  incompatible  with  the  new 
library.  This  change  was  made  deliberately  in  order  to  improve  the  paper  trail  leading  from  output  bade 
to  the  library  on  which  it  was  based. 
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Intentionally  left  blank. 
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7.  Difference  in  Quantity  (old  lib  vs.  oxtondod  lib)  !  7.  Difference  in  Quantity  (old  lib  vs.  oxtondod  lib) 


Standard  (Old)  Library  vs.  Extended  library  SO 


■is 

2S 


- - - 1 - 

i  i 

_ - - 

..**i  i  i 

.*i  i  l  l 


3j0  3S 


4S  SjO 


Temperature/ 1000  (K) 

Fig.  1  -30;  M30  (no  sddod  Inorganics)  -  17  Rejects 


Standard  (Old)  Library  vs.  Extended  Library  40 


SJO  3S  4JD  4S  SS  SS 

Temperature/ 1 000  (K) 

Fig.  t  -40:  M30  (no  inorganic  sddKvos)-  1 7  Reject* 


Figure  1  -A  M30  (no  inorganic  additives)  -17  Rejects 


47 


X  Difference  in  Quantity  (old  lib  vs.  oxtoodtd  lib)  !  X  Difference  in  Quantity  (old  lib  vs.  oxtsndsdW 


Standard  (Old)  Library  vs.  Extandad  Library  SO 


SjO  3£  40  4£  SJO  SJ5  60 

Temperature/ 1000  (K) 

F)g.S-30:  JA-2  - 17  W|»ct» 


Standard  (Old)  Library  vs.  Extandad  Library  40 


SO  3£  40  4£  SO  S£  60 

Temperature/ 1000  (K) 

Fig.  240:  JA-2  - 17  Bejectt 


Figure  1  -B:  JA-2  —  17  Rejects 
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Added  Energy  (J/g)  J  Added  Energy  (J/g) 


Z  Difference  in  Quantity  (old  lib  vs.  oxtondod  lib)  !  Z  Difference  in  Quantity  (old  lib  vs.  oxtondod  lib) 


Standard  (Old)  Library  m*  ffilandad  Library  SO 


SJD  3S  4J0  4JS  SJO  SJS  SJO 

Temperature/ 1000  (K) 

Fig l  MO;  MS  fee  kwgonlr  addUNop)  - 17  Wsjscte 


Standard  (Old)  Library  w  Extended  Library  40 


Temperature/ 1000  (K) 

WgMO.MS  (mo  Inorganic  addMvoo)  -17 


Figure  1  - C :  MS  (no  inorganic  additives)  - 1 7  Rejects 


49 


Added  Energy  (J/g)  '  Added  Energy  (J/g) 


7.  Difference  in  Quantity  (old  lib  vs.  extended  lib)  ;  7.  Difference  in  Quantity  (old  lib  vs.  extended  lib) 


Standard  (Old)  Library  ml  Extended  Library  30 


Temperature/ 1000  (K) 

Fig.  4-30:  119  (no  added  Inorganic*)  - 17  Refects 


Standard  (Old)  Library  vs.  Extended  Library  40 


Figure  1  -Or  M9  (no  inorganic  additives)  - 1 7  Rejects 
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Added  Energy  (J/g)  !  Added  Energy  (J/g) 
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X  Difference  in  Quantity  (ota lib  vs.  txfndtd  lib)  \  X  Difference  in  Quantity  (old lib  vs. 


Z  Difference  in  Quantity  (old  lid  vs.  oxtondod  lib)  i  Z  Difference  in  Quantity  (old  lib  vs. 
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X  Difference  in  Quantity  (old  lib  vs.  extended  lib)  I  X  Difference  in  Quantity  (oMM»  vs. 


Standard  (Old)  Ubrary  i «,  Extended  Library  SO 


Temperature/ 1000  (K) 

/fe  M80  (no  added Inorgeniet)  -  No  Rejects 


Figure  2-A:  M30  (no  inorg  additives)  --  Rejects  vs.  No  Rejects 
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Added  Energy  (J/g)  J  Added  Energy  (J/g) 


X  Difference  in  Quantity  (old  Hb  vs.  oaondtd  lib)  !  X  Difference  in  Quantity  (otd  ttb  vs. 
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Z  Difference  in  Quantity  (old  lib  vs.  oxtondod  lib)  i  Z  Difference  in  Quantity  (old  lib  vs.  o 


Library  SO 


Temperature/1000  (K) 

Hg.1000:  MS  (no Inorganic adiMh’ma)  -No  Rejects 


Figure  2-C:  M5  (no  inorg  additives)  ~  1 7  Rejects  vs  No  Rejects 
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Added  Energy  (J/g)  |  Added  Energy  (J/g) 


7.  Difference  in  Quantity  (old  lib  vs.  oxtondod  W  !  Z  Difference  in  Quantity  (S old  Mb  vs. 


Standard  (Old)  Library  w  Extandad  Library  SO 


33  40  43  SJO  53  BO 

Temperature/ 1000  (K) 

Fig.  1 1  -30:  M«  (mo  Inofffoale  oddbbmo)  -  No  Ha)aeta 


Figure  2-D:  M9  (no  inorg  additives)  - 1 7  Rejects  vs  No  Rejects 
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Added  Energy  ( J/g )  |  Added  Energy  (J/g) 


7.  Difference  in  Quantity  (eld  lib  vs.  extended  lib)  !  Z  Difference  in  Quantity  (old  lib  vs.  extended  lib) 


Standard  (Old)  Library  re.  Extended  Library  SO 


Temperature/1000  (K) 

mg.  Ml  -  No  ftojacta 


Figure  2-E:  Ml  -17  Rejects  vs  No  Rejects 
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X  Difference  in  Quantity  (old  lib  vs.  oxtondod  lib)  !  X  Difference  in  Quantity  (old  lib  vs. 


Figure  2-F:  Decalin  +  WFNA  -  1 7  Rejects  vs  No  Rejects 
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7  Difference  in  Quantity  (old lib  vs.  oxtondod  lib)  !  7  Difference  in  Quantity  (old  lib  vs. 


Standard  (Old)  Library  v*  Extended  Ubrary  SO 


Temperature/ 1000  (K) 

Fig.  14-00:  Decade  ♦  HiOt  -  We  Hejecte 


Figure  2-G :  Decalin  +  H2O2  -  1 7  Rejects  vs  No  Rejects 
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Appendix  A. 

Listing  of  TESTLIBS.BAT 

REM  A  batch  file  to  compare  old  &  new  McBride  libraries  up  to  6,000  K 
@ECHO  Off 

ton  11  September  1992 

tern  This  batdi  file  tuns  EXMUGET  for  4  cases: 

tern  o  without  or  with  ions  permitted,  and 

rent  o  using  old  (standard)  and  new  (extended)  thermo  libraries. 

rem  The  following  8  files  are  produced: 

rem  o  tstlbsno.old  and  tstlbsio.old  (summary  w/  old  lib ) 

rem  *no  ions  'ions 

rem  o  tstlbono.old  and  tstlboio.old  (output  "  "  "  ) 

rem  o  tstlbsno.new  and  tstibsio.new  (summary  w/  new  lib ) 

rem  o  tstlbono.new  and  tstlboio.new  (output  "  "  "  ) 

rem  Then  the  summary  files  are  operated  on  by  COMP ARE2 .BAS,  producing: 

rem  o  summary  1. tab:  compares  no  ions/ions  &  old  library; 

rem  o  summary2.tab:  compares  no  ions/ions  &  new  library; 

rem  o  summary3.tab:  compares  old  &  new  libraries,  no  ions; 

rem  o  summaiy4.tab:  compares  old  &  new  libraries,  with  ions; 

rem  Start  timer 
c^'nortonXtm  start  start  /log  fl 
echo. 

rem  Is  only  a  comparison  wanted? 
if  «%r  =  "Q"  goto  quikk 

rem  Verify  whether  new  binary  libraries  are  to  be  created 

if  "%i"  =  ""  goto  RUN1 

if  "%1"  =  "C  goto  makelib 

if  not  "»1"  =  "c"  goto  RUN1 

:MAKEUB 

els 

echo. 

echo. 

echo  Confirm  that  you  warn  a  CREATION  run,  or  press  *C  to  abort: 

echo. 

pause 

rem  Start  with  standard  FILE4.DAT;  convert  it  to  new  format 
echo. 

echo  Copying  standard  FILE4.DAT  to  EXFILE4.DAT 

copy  dtefisource\libraryvfile4.dat  file4.dat  >  nul: 

echo. 

echo  Converting  FILE4.DAT  to  ENTFILE.DAT 

@echo  on 
oldint  file4.dat 
@echo  off 
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echo. 

echo  Converting  INTFILE.DAT  to  NEWFILE.DAT 

<g>echo  on 

intnew 

<§>echo  off 

echo. 

echo  Add  identifying  header  to  NEWFILEDAT 

copy  old_lib.hed  +  newfile.dat  tempo  >  nul: 

del  newfile.dat 

ran  tempo  newfile.dat 

echo. 

echo  Now  form  binary  library 

copy  libjnmin  input  >  nul: 

@echoon 
exmucet  newfile.dat 
del  newfile.dat 
del  intfile.dat 
del  file4.dat 
(©echo  off 
echo. 

echo  Save  binary  library  for  future  use. 

copy  exthnnlb.dat  newstdJib  >  nul: 

rem  Now  repeat  with  McBride’s  extended  library 
echo. 

echo  Copying  McBride’s  extended  FILE4.DAT  to  EXFILE4.DAT 

copy  d^f\sourceMibrary'exfile4.dat  exfile4.dat  >  nul: 
echo. 

echo  Now  form  binary  library 

copy  lib_run.in  input  >  nul: 

@echo  on 
exmucet  exfilc4.dat 
(©echo  off 
echo. 

echo  Save  Unary  library  for  future  use. 

copy  exthnnlb.dat  exfile4  Jib  >  mil: 
goto  nin2 

:RUN1 

els 

echo. 

REM  echo  Confirm  NO  need  to  (re)create  libraries,  or  press  X?  to  abort 
echo. 

REM  pause 

rem  Start  here  if  libraries  have  previously  been  formed 
:RUN2 
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echo  Copy  standard  binary  libary  to  EXTHRMLB.DAT 

copy  newstdJib  exthnnlb.dat  >  nul: 

echo  Copy  first  set  of  test  cases  to  INPUT 

copy  NcefexibimVno JonsJn  input  >  nul: 

echo. 

echo  ...and  now  nm  EXMUCET  using  OLD  library 

echo. 

@echo  on 
exmucet 
@echo  off 
echo. 

echo  Save  output 

copy  summary  tstlbeno.old  >  rail: 
copy  output  tstlbooo.old  >  nul: 

echo. 

echo  Now  copy  new  set  of  test  cases  (thermo  library  unchanged) 
echo. 

echo.  Copying  ions  in  to  input 
copy  ionsin  input  >  nul: 
echo. 

echo  Now  re-iun  EXMUCET 

@echo  on 

exmucet 

@echo  off 

echo. 

echo  Save  output 

copy  summary  tstlbsio.old  >  nul: 
copy  output  tstlboio.old  >  mil: 

rem  Now  repeat  these  calculations  using  McBride's  extended  library 
echo. 

echo  Copy  previously-formed  new  binary  library 

copy  exfile4  Jib  exthnnltxdat  >  nul: 

echo. 

echo  Run  the  same  test  file  again 

echo. 

exmucet 

@echo  off 

echo. 

echo  Save  output  for  comparison 

copy  output  tstlboiojrew  >  mil: 
copy  summary  tstibskuiew  >  nul: 

echo. 

echo  Copy  FIRST  test  of  test  cases... 
echo. 


copy  no  Joskin  input  >  nul: 
echo. 

echo  ...  and  ran  EXMUGET 

exmucet 

Gecfao  off 

echo. 

echo  Save  output  for  comparison 

copy  output  tsfibonojiew  >  nul: 
copy  aummary  tadbnojiew  >  nul: 
ran  Stop  timer 

ctaortoriStm  stop  done  /I  >  dummy 

cfNnoitortan  stop  done  /I 

copy  dummy  df>ce^exfotm\iime_nm  >  nul: 

del  dummy 

echo. 

echo  Delete  scratch  files  (if  any) 
del  a*.* 

:quikk 

echo. 

echo.  Run  COMPARED  AS  4  times 
echo. 

copy  testlibl  Jn  testjopt  >  nul: 
qb/run  'cetexfonnNaourcetaxnpare 
copy  summaiy.tab  summary  1.  tab  >  nul: 
copy  testlib2.in  test.opt  >  nul: 
qb/run  \cetexform\souiceVompare 
copy  summary.tab  summary2.tab  >  nul: 
copy  tcsUib3Jn  testjopt  >  nul: 
qb/nm  VcebexfohnNsourceNcompare 
copy  summary.tab  summary3.tab  >  mil: 
copy  tesdib4Jn  testjopt  >nut 
qb/run  VcetexfonnNsource^oinpare 
copy  summary  .tab  summary4.tab  >  nul: 

if  "%1"  «  "Q"  goto  qq 
echo. 

etho  For  safety’s  sake,  ensure  that  EXTHRMLB.DAT  is  (new)  extended  library 

copy  exfilc4Jib  exthrmlb.dat  >  nul: 

echo. 


•qq 

if  "%r  «  "Q"  caU  wdstar  dfcetexfonn\summary3.tab 
if  "%i”  «os  "Q"  goto  endd 

call  wdstar  d^cetexform\time._nin 

:endd 

d: 

ctfceflexfbnn 
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Appendix  B. 

Listing  of  PROGRAM  MAKE JNP.BAS 


NOTE:  In  order  to  Improve  the  legibility  of  the  listings  of  die  QUICK  BASIC  prograns  given  in  this 
report,  the  underscore,  has  been  used  as  a  continuation  symbol  in  some  lines  that  would  otherwise 
be  too  long. 

EXAMPLE:  A  line  such  as 

PRINT  #2,  num.kasc2(kisc);  temp2(kase);  gmptus2(kase);  V 
will  be  printed  here  as 

PRINT  #2,  num.kase2(ka«c);  V;  temp2(lcase);  V;  _ 
gmptus2(kase);  V 


*  A  program  to  make  huge  run  streams  for  checking  EXMUCET 

*  Written  by  EF&A  in  Aug/92 

’  Saved  as  MAKEJNP.BAS  (S  September  1992] 

DIM  e.zeiT>(15) 

DEFINT I-N:  DEFSNG  A-H,  O-Z 


DATA  -1538.79  ,  -1587.72  ,  -2488.30  ,  -2302.44,  -2103.22 
DATA  -2316.85  ,  -2625.75  .  -5145.85  ,  -1573.22 
DATA  -99999.,  0.  0.  0.  0 . 0 


FOR  i  =  1  TO  15:  READ  e.zeroO):  IF  (ejeroCO  *  -99999!)  THEN  num.e  =  i  -  I 
NEXT 

GLS  :  FOR  i  =  1  TO  7:  PRINT  :  NEXT 

INPUT  "Enter  in  order.  E  start.  Delta  E,  and  E  final:  ”,  _ 

energyz,  delta.e,  e .final 


FOR  Linp  =  1  TO  2 
CLOSE  #1:  CLOSE  #2 
IF  (tinp  =  1)  THEN 

file.inputtS  =  "dfscefiexfonn\template.r 
file.outputS  *  "d^cef«xibnn\noJons.in" 

ELSE 

file.inputt$  as  "dtetexformNtemplate.2H 
file.outputS  «  "d.VetNexformNions.in" 
END  IF 

OPEN  file.outputS  FOR  OUTPUT  AS  #2: 
energy  =  0! 
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recycle  :  OPEN  He  Ji^iS  FOR  INPUT  AS  #1 
IF  (energy  o  0!)  THEN  GOTO  tfara 
PRINT"  "; 

FOR  i  *  1  TO  INT(1. 05  *  eiinal  /  deltas  +  1):  _ 

PRINT  enerzy  delta*  *  0  - 1);  V; :  NEXT: 

PRINT 

PRINT  #2."  "; 

FOR  i  *  1  TO  INT(1.QS  *  (e.final  /  delta*)  +  1):  _ 

PRINT  #2,  enerzy  +  delta*  •0*1): 

NEXT:  PRINT  #2,  "  " 

PRINT  #2."  -9999." 

energy  *  enerzy 
start!'  kk  =  0 

rcadl:  LINE  INPUT  #1.  lyneS:  IF  (EOF(l))  THEN  GOTO  new.energy 
FOR  i  *  1  TO  LEN(lyne$) 

IF  (MIDSOyneS,  i.  1)«  *#")THEN 
num  « i  - 1:  GOSUB  makeline  1:  GOTO  read2 
END  IF 
NEXT 
print  lyneS 

PRINT  #2.  lyneS:  GOTO  rcadl 
rcad2:  LINE  INPUT  #1.  lyneS 
FOR  i  «  1  TO  LENOyneS) 

IF  (MIDSOyneS.  i,  2)«  "U«"  AND  MIDSOyneS.  i  +  2. 1)  -  "#")  THEN 
num  ■ i  - 1:  GOSUB  makeline2:  GOTO  readl 
END  IF 
NEXT 

’  print  lyneS 

PRINT  #2.  lyneS:  GOTO  rcar*2 

new.energy:  energy  *  energy  +  delta*:  PRINT  "Energy  = ";  energy 
DF  (energy  >  e.final)  THEN 
GOTO  nexxt2 
ELSE 

CLOSE  #1:  GOTO  recycle 
END  IF 
nexxt2:  NEXT 
SYSTEM 

makelinel: 

’  a  subroutine  to  insert  die  added  energy  into  the  tide  line 
dumS  =  LEFTSOyneS,  num) 

dum$  *  dumS  +  STR$(energy)  +  "  J/g  Elec":  PRINT  #2,  dumS: 

PRINT  dumS:  RETURN 

makeline2: 

*  a  subroutine  to  insert  die  revised  energy  into  die  namelist  line 
kk*kk+  1 

dumS  *  LEFTSOyneS,  num) 

dumS  *  dumS  +  "U»"  +  STRS(energy  +  e*ero(kk))  +  T 
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PRINT  #2,  dumS: 
PRINT  dumS 
RETURN 


Intentionally  left  blank. 


Appendix  C. 

Template  of  the  Teat  Cases  for  Compering 
Different  Thermodynamics  Lfonriec 


TITLE 

1.  Pseudo M30- NO  IONS  ALLOWED  +  #  J/g  Added 


REACTANTS 


C  6. 

H  7349 

0  9.901 

N  2.451 

27.90 

•169.16E3S298.15 

F 

C  3. 

H  5. 

0  9. 

N  3. 

22.42 

-8.86E4  S298.15 

F 

C  1. 

H  4. 

N  4. 

01 

47.54 

-221E4  S298.15 

F 

C  17. 

H  20. 

N  2. 

0  1. 

1.49 

-2J1E4  S298.15 

F 

Cl. 

.10 

aO  S298.15 

F 

C  2. 

H  6. 

O  1. 

25 

-66.42E3L298.15 

F 

NAMELISTS 

&INPT2  UV«TJRH0-02,  SIUNIT-TJONS-F.U-#/ 
TITLE 

2.  M30A1  -  NO  IONS  ALLOWED  +  «  J/g  Added 


REACTANTS 


C  6. 

H  7.739 

0  9.521 

N  2.261 

27.9 

- 173.77E3S298.15 

F 

C  3. 

H  5. 

0  9. 

N  3. 

22.42 

-88.6E3  L298.15 

F 

Cl. 

H  4. 

02. 

N  4. 

46.84 

-22.1E3  S298.15 

F 

C  17. 

H  20. 

0  1. 

N  2. 

1.49 

-25.1E3  S298.15 

F 

C  2. 

H  6. 

0  1. 

.25 

-66.42E3  L298.15 

F 

Cl. 

.1 

0.  S298.15 

F 

NAMELISTS 

&INPT2  UV=T.RHCWL2.  SIUNnr«T40NS«F,U-#/ 


TITLE 

3.  Pseudo  WC890  —  NO  IONS  ALLOWED  ♦  #  J/g  Added 


REACTANTS 

C  6. 

H  7.549 

0  9.901 

C  3. 

H  5. 

0  9. 

C  12. 

N  1. 

H  11. 

C  16. 

H  22. 

0  4. 

C  7. 
Cl. 

H  6. 

0  4. 

N  2.451 

79.98 

N  3. 

10.2 

1.11 

7.65 

N  2. 

.08 

.14 

-169.16E3S298.15  F 
-8.86E4  S298.15  F 

3.107E4  S298.15  F 

-201. E3  S298.15  F 

.1.71E4  S298.15  F 

0.0  S298.15  F 


NAMELISTS 

&INPT2  UV*T,RHO><X2,  SIUNIT*T,I0NS-F,U*4/ 
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TITLE 

4.  JA-2  (FRG)  ~  NO  IONS  ALLOWED  +  #  J/g  Added 


REACTANTS 


C  6. 

H  7.395 

0  10.209 

N  2.605 

59.5 

-165.41E3S298.15 

F 

C  3. 

H  5. 

0  9. 

N  3. 

14.9 

•88.6E3  L298.15 

F 

C  4. 

H  8. 

0  7. 

N  2. 

24.8 

-103.56E3S298.15 

F 

C  2. 

H  6. 

0  1. 

.05 

-66.42E3  L298.15 

F 

C  14. 

H  14. 

N  2. 

0  1. 

.7 

-255X3  L298.15 

F 

NAMELISTS 

&INPT2 

!  UV=TJRH(WL2,  SIUNIT*TJONS»F,U»#  / 

TITLE 

5.  M9- 

NO  IONS  ALLOWED  +  «  J/g  Added 

REACTANTS 

C  6. 

H  7.329 

0  10.341 

N  2.671 

57.55 

-163.81E3S298.15 

F 

C  3. 

H  5. 

0  9. 

N  3. 

39.86 

-88.6E3  L298.15 

F 

C  17. 

H  20. 

O  1. 

N  2. 

.75 

-25.1E3  S298.15 

F 

C  2. 

H  6. 

0  1. 

.35 

-66.42E3  L298.15 

F 

NAMELISTS 

&INPT2 

1 

m 

l 

1 

i 

TITLE 

6.  Ml- 

NO  IONS  ALLOWED  +  *  J/g  Added 

REACTANTS 

C  6. 

H  7.364 

0  10.271 

N  2.636 

83.11 

-164.65E3S298.15 

F 

C  7. 

H  6. 

0  4. 

N  2. 

9.77 

-17.1E3  S298.15 

F 

C  16. 

H  22. 

0  4. 

4.89 

-201.4E3  L298.15 

F 

C  2. 

H  6. 

0  1. 

.05 

-66.42E3  L298.15 

F 

C  12. 

H  11. 

N  1. 

.98 

31.07E3  L298.15 

F 

H  2. 

O  1. 

JO 

-285831.  L298.15 

JF 

NAMELISTS 

&INPT2  UV=T,RHO=0.2,  SIUNIT=T,IONS=FIU=#/ 


TITLE 

7.  DECAUN  +  WFNA  -  NO  IONS  ALLOWED  +  #  J/g  Added 


REACTANTS 

HI.  N 1.  0  3. 

C  10.  H  18. 

N  1.  0  2. 

H  2.  0  1. 


.8088  -41.61E3  L298.15 

.178  -55.1E3  L298.15 

.02466  31711.  G298.15 

.010686  -285831X298.15 
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si  si  *«  *« 


NAMELISTS 

&INPT2  UV-T.RHO-O.2,  SIUN1T*TJ0NS=F,U»#/ 


TITLE 

8.  DECALIN  +  70%  HYDROGEN  PEROXIDE  -  NO  IONS  ALLOWED  +  •  J/g  Added 


REACTANTS 
C  10.  H  18. 

H  2.  0  2. 

HI  0  1. 


1.  -55.1E3  L298.15  F 

.7  -44.88E3  L298.1S  O 

.3  -68.315E3L298.15  O 


NAMELISTS 

&INPT2  UY*TJtH0-O2.  SIUNIT*TJONS*=F,U=#/ 


TITLE 

9.  RFNA  +  UDMH  -  NO  IONS  ALLOWED  +  #  J/g  Added 
REACTANTS 

N  1.6294  H  13722  O  4.6950  .7  -64860.  L298.15  F 

C  2.  H  8.  N  2.  .3  11.9E3  L298.15  F 


NAMELISTS 

&INPT2  UV=T,RHO=0.2,  SIUNIT=TJONS*=F,U=#/ 
STOP 
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Intentionally  left  blank. 
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Appendix  D. 

Listing  of  Program  COMPARE.BAS 

‘  A  program  to  compere  output  files  from  EXMUCET 
’  Written  by  EF AA  in  Oct/92 

1  Saved  as  DAEXMUCETCOMPARE.BAS 

*  This  version:  01-MAR-1994 

-  Error  in  file  names  in  DATA  statement  corrected 

-  Changes  begun  on  13-FEB-94  made  internally  consistent  A  completed 
~  Some  cosmetic  changes  made. 

’  13  February  1994: 

*  On-screen  directions  added. 

*  Use  of  external  file  for  options  discarded. 

’  Hie  names  changed. 

’  9  February  1994: 

*  File  names  changed 

*  Modified  [27-OCT-93]  to  remove  bugs  and  to  simplify  the  programming; 

’  The  output  was  improved. 

’  Modified  [11 -SEP-92]  to  indude  identification  of  foe  thermo 
’  libraries  used  to  generate  foe  outputs  being  compared. 

CLS  :  DEFINT I-N:  DEFSNG  A-H,  O-Z 
DIM  templ(999),  gmptus 1(999),  balleigy 1(999) 

DIM  num.kase 1(999),  num.kase2(999),  e.zero(99) 

DIM  temp2(999),  gmptus2(999),  ballergy2(999) 

DIM  Delta.g(999),  Delta.pct.t(999),  Deha.pcLb<999) 

DIM  Delta.t(999).  Ddta.b(999) 

energy  .save  =  -1!:  no.ener  «  Ch  no.cmpsns  =  8 

*  Determine  type  of  ran.  There  are  3  types:  Compare  results  from 

*  McBride's  standard  library  without  ions  or  with  ions  permitted,  and 
’  compare  foe  standard  library  with  foe  extended  one  without  or  with 

*  ions  permitted. 

get.  option: 

COLOR  1. 11:  CLS  :  LOCATE  3, 7 

PRINT  "  :  COLOR  1,  7:  PRINT  "TYPE  1"; :  COLOR  1, 11 

PRINT  "  Compare  results  from  the  standard  library  with 

PRINT  SPACER  17);  "results  from  foe  extended  library,  no  ions  permitted." 
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LOCATE  d.  7 

PRINT  "  "; :  COLOR  1.  7:  PRINT  "TYPE  2"; :  COLOR  1.  11 

PRINT  "  Compare  results  from  die  standard  library  with 

PRINT  SPACE$(17);  "results  from  the  extended  library,  ions  permitted." 

LOCATE  9.  7 

PRINT  "  "; :  COLOR  1.  7:  PRINT  "TYPE  3"; :  COLOR  1.  11 
PRINT "  Compare  results  from  the  extended  library 
PRINT  SPACE$(17);  "without  or  with  ions  permitted.” 

LOCATE  14. 17:  COLOR  13. 9:  PRINT  "Enter  your  option  1. 2.  or  3:  " 
wain: 

a$  *  INKEY$:  IF  (LEN(a$)  >  0)  THEN  GOTO  wain 
CLS  :  opt  *  VAL(a$) 

IF  (opt  <  1  OR  opt  >  3)  THEN  CLS  :  PRINT  CHR$(7):  GOTO  getoption 

COLOR  13. 5:  LOCATE  9. 23:  PRINT  "  Option  chosen  opt 
LOCATE  11. 19:  PRINT  "WAIT.  Data  are  being  processed.  ” 

COLOR  15. 9 
IF  (opt  =  1)  THEN 
options  *  _ 

"  »>  EXTENDED  LIBRARY  vs.  STANDARD  LIBRARY.  NO  IONS  PERMITTED  «<" 

READ  filel  jRimS.  file2.sum$ 

READ  filel.outS.  file2.outS 
outextS  *  M" 

ELSEIF  (opt  =  2)  THEN 
options  =  _ 

"  »>  EXTENDED  LIBRARY  vs.  STANDARD  LIBRARY.  IONS  PERMITTED  «<" 

FOR  i  =  1  TO  4:  READ  JUNKS:  NEXT 
READ  filel  jum$,  file2.sum$ 

READ  filel.outS.  file2.out$ 
outextS  *  n2K 
ELSE 

options  *  _ 

"  »>  EXTENDED  LIBRARY:  IONS  EXCLUDED  vs.  IONS  INCLUDED  <«" 

FOR  i  s  1  TO  8:  READ  JUNKS:  NEXT 
READ  filel  .sumS,  fik2.sum$ 

READ  filel.outS.  file2.out$ 
outextS  -  ".3" 

END  IF 
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file$  *  "OUTPUT"  +  ouLextS 
OPEN  fileS  FOR  OUTPUT  AS  #9 
fileS  *  TNTERMED"  +  ouLextS 
OPEN  fileS  FOR  OUTPUT  AS  #33 
fileS  *  "BRIEF  +  ouLextS 
OPEN  fileS  FOR  OUTPUT  AS  #253 
’PRINT  #9."  Option  chosen." 

PRINT  #9.  options 

PRINT  #9, 

PRINT  #9.  _ 

"  The  two  EXMUOET  output  files  are  filel.outS; "  and  file2.out$ 
PRINT  #9, 

GOSUB  lib.name:  ’  Try  to  detennine  which  library  was  used 
PRINT  #9.  _ 

"  The  input  (SUMMARY)  files  are  ";  filel.sumS; "  and  ";  file2.sum$ 
PRINT  #9, 


PRINT  #253. :  PRINT  #253,  "  First  EXMUCET  output  file  =  ";  filel.outS 
PRINT  #253,  "  Second  EXMUCET  output  file  *  ";  file2.out$:  PRINT  #253, 
PRINT  #253.  "  Input  (SUMMARY)  files  are  ";  filel.sumS; "  and  " 

PRINT  #253.  "  ";  file2.sum$:  ’  PRINT  #253, 

’PRINT "  Input  (OUTPUT)  files  are  filel.outS;  "  and  " 

’PRINT  "  ";  file2.outS:  PRINT  "  " 

PRINT  #9.  "  Thermo  library  used  for  first  set  was  " 

PRINT  #9.  TAB(19);  lib4tame.lS:  PRINT  #9,  "  " 

PRINT  #9,  "  Thermo  library  used  for  second  set  was  " 

PRINT  #9,  TAB(19);  libJtanv..2S 

’  Initialization  for  first  file 

OPEN  filel.sumS  FOR  INPUT  AS  #1 

OPEN  TEMPOl”  FOR  OUTPUT  AS  #2 

LINE.NO  =  0:  LK  -  0:  kase  =  0:  num. of. energies  =  0 


'  Read  a  line  of  first  input  file 
READ1:  LINE  INPUT  #1,  lyneS 

IF  (EOF(l))  THEN  GOTO  NEXT  FILE 
IF  (LEFTS(lyne$,  5)  =  "  Finis")  THEN 
PRINT  CHR$(7);  CHRS(7);  CHRS(7) 
GOTO  NEXT.FILE 
END  IF 
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’  Look  for  sentinel  (’>’) 

CONTTX:  IF  (LEFT$0yne$,  2)  o  "  >")  THEN  GOTO  READ1 

’  Sentinel  has  been  found.  Increment  case  counter, 
kase  =  lease  +  1 
GOSUB  get  energy 

’Next,  skip  3  lines 

LINE  INPUT  #1,  JUNK1J:  LINE  INPUT  #1,  JUNK2$:  LINE  INPUT  #1,  JUNK3S 
'  Read  line  of  data,  and  extract  temperature,  impetus,  and  ballergy 
LINE  INPUT  #1.  lyneS 

tempi  (kase)  *=  VAL(MID$<lyne$,  25, 9)) 
gmptusl(kase)  =  VAL(MID$Oyne$,  49,  8)) 
ballergy  1  (kase)  =  VAL(MID$Oyne$,  63.  8)) 

’  PRINT  "1  ";  kase;  num  .kase  1  (kase);  tempi  (kase);  gmptusl(kase);  _ 
ballergyl(kase) 

PRINT  #2,  num.kasel(kase);  ",";  tempi  (kase);  Y;  gmptusl(kase);  _ 

ballergy  1  (kase) 

GOTO  READ1 

’  When  end  of-of-file  has  been  found,  write  sentinel  at  end  and  close  files 
NEXT.FHJE:  num.nrgy  «  kase:  ’  PRINT  "LAST  KASE  =";  num.nrgy 
PRINT  #2.  "-9999..  -9999.,  -9999.":  PRINT  #2,  "  ":  CLOSE  #1:  CLOSE  #2 

*  Initialization  for  second  file 
OPEN  file2.sumS  FOR  INPUT  AS  #1 
OPEN  "TEMP02"  FOR  OUTPUT  AS  #2 


kase  -  0 
READ2: 

LINE  INPUT  #1,  lyne$: 

’  IF  (EOF(l)  OR  LEFT$(LYNE$,  5)  =  "  Finis")  THEN  GOTO  MAKE.TAB 

IF  (EOF(l))  GOTO  MAKE.TAB 
IF  (LEFTS(lyne$,  5)  =  "  Finis")  THEN 
PRINT  CHR$(7);  CHR$(7);  CHR$(7) 

GOTO  MAKE.TAB 
END  IF 

IF  (LEFTSOyneS,  2)  o  "  >")  THEN  GOTO  READ2 
kase  =  kase  +  1 
num.kase2(kase)  =  kase 

LINE  INPUT  #1,  JUNKS;  LINE  INPUT  #1,  JUNKS:  LINE  INPUT  #1,  JUNKS 
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LINE  INPUT  #1,  lyneS 

temp2 (lease) »  VAL(MH>$(lyne$,  25, 9)) 
gmptus2(kase)  =  VAL(MID$(lyne$,  49.  8)) 
ballergy2(kase)  -  VAUMIDSOyneS,  63.  7» 

’PRINT  "2  lease;  num.lease2(kase);  temp2(kase);  gmptus2(kase);  _ 
ballergy2  (lease) 

PRINT  #2.  num.lease2(kase);  V;  temp2(kase);  gmpcus2(lease); 

balkrgy2(kase) 

GOTO  READ2 

MAKE. TAB:  PRINT  #2.  "-9999.,  -9999.,  -9999.":  PRINT  #2,  "  " 
CLOSE  #1:  CLOSE  #2 

’  PRINT  "-9999.,  -9999..  -9999.";  nunuirgy 

HUNT  #9, :  PRINT  #9,  "  Number  of  energies  *  no.ener  *  HUNT 

’  Gel  ready  to  make  summary  table  (previously-opened) 

OPEN  "  tempo  1"  FOR  INPUT  AS  #1 
OPEN  "tempo2"  FOR  INPUT  AS  #2 
’  PRINT  "Run  option  chosen  =  imoptionS 
IF  (opt  o  3)  THEN 


head3$  =  _ 

"  T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  _ 

Gam  (Ext)" 

head33$  *  _ 

"  T(Std)  T(Ext)  Delta  T  %DelT  Ball(Std)  Ball(Ext)  DeltaB  %DelB" 
ELSE 
head3$  =  _ 

"  T(no)  T(ions)  Imp(no)  Imp(ion)  Ball(no)  BallQon)  Gam(no)_ 

Gam(ion)" 

head33$  =  _ 

"  T(no)  TOons)  Delta  T  %DdT  Ball(no)  Ball(ion)  Delta  Ball  %DdB" 
END  IF 

’  Read  a  line  of  data  from  file  1 
FOR  i  =  1  TO  999 

INPUT  #1.  num.kasel(i),  templ(i).  gmptusl(i),  ballergyl(i) 

IF  (tempi  (i)  >  0)  THEN  GOTO  NEXTT1  ELSE  GOTO  CONTT1 
NEXTT1:  NEXT 
CONTT1:  num.l  =  i  -  1 

’  Read  a  line  of  input  from  file  2 

FOR  i  =  1  TO  999:  INPUT  #2,  num.kase2(i),  temp2(i),  gmptus2(i),  ballergy2(i) 
IF  (temp2(i)  >  0)  THEN  GOTO  NEXTT2  ELSE  GOTO  CONTT2 
NEXTT2:  NEXT 
CONTT2:  NUM2  =  i  -  1 
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’  Check  that  same  number  of  lines  have  been  read 

IF  (num.l  »  NUM  .2  AND  num  kase  1  (num Jirgy)  «  n:in.Jtase2(numjurgy))  _ 

THEN  GOTO  NEXTT3 

CLS  :  PRINT  TROUBLE!!!  NUM.1  o  NUML2";  NUM.1;NUM2. :  STOP 
NEXTT3: 

*  Print  summary  table  with  statistics  after  each  group 
FOR  i.ener  *  1  TO  no.ener 


PRINT  #9. :  PRINT  #9,  TAB(30);  ejzero(i.ener);  'kJ/g  added" 
PRINT  #253. :  PRINT  #253.  TAB(25);  e^eroCi.ener);  "kJ/g  added" 
PRINT  #33. :  PRINT  #33.  TAB(30);  e.zero(i.ener);  "kJ/g  added" 
PRINT  #33.  head33$ 

PRINT  #9.  bead3$ 


'  FOR  jj  =  noxmpsns  *  i.ener  •  (no.cmpsns  - 1)  TO  no.cmpsns  *  i.ener 
FOR  jj  =  no.cmpsns  *  (i.ener  - 1)  +  1  TO  no.cmpsns  *  i.ener 
num.case  =  jj  MOD  no.cmpsns:  IF  num. case  =  0  THEN  num  .case  =  no.cmpsns 


PRINT  #9,  USING  _ 

||  um\  iyyy.il  u  aa  iiuiiuii  at 


ujyyyi  aa  aa  a^a^ata^at  at  uiwyi 

.W  Wn.nnnn 


##.####";  num.case;  tempi  (jj);  temp2(jj);  gmptusl(jj);  gmptus2(ij);  _ 
ballergylQj);  baUeigy2(jj);  1!  +  gmptusl(jjj)  /  hallergylQj);  _ 

1!  +  gmptus2(jj)  /  ballergy2(jj) 


PRINT  #33,  USING  _ 

Niyi\  uuu^^  u  uuu^^  u  uii^m  u  ujyiu 


uiyiuii  u  uu^uu  u  uuhiju  u  uu  ^<n( 
l^^nnn.n 


num.case;  templ(jj);  temp2(jj);  temp2(jj)  -  tempi  (jj);  _ 
100!  *  (temp2(ij)  -  templQj))  /  tempi  (jj);  ballergyl(jj);  _ 
ballergy2(jj);  ballergy2(jj)  -  ballergyl(jj);  _ 

100!  *  (balleigy2(u)  -  ballergyl(ij))  /  ballergyl(jj) 
NEXT  jj 
GOSUB  STAT 


NEXT  i.ener 


’  Prepare  and  print  table  of  differences 
PRINT  #33. 

PRINT  #9.:  HUNT  #9,  TAB(26);  "GRAND  TOTALS":  PRINT  #9. 

PRINT  #9.  TAB(26);  "MEAN  ABSOLUTE  DEVIATIONS" 

PRINT  #9.  TAB(13);  "Temp  (K)  Impetus  (J/g)  Ballergy  (J/g)" 
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PRINT  #33, :  PRINT  #33.  TAB(26);  “GRAND  TOTALS":  PRINT  #33, 
PRINT  #33.  TAB (26);  "MEAN  ABSOLUTE  DEVIATIONS" 

PRINT  #33.  TAB(13);  "Temp  (K)  Impetus  (J/g)  Ballergy  (J/g )" 

sum.tl  *  0!:  sum.t2  *  0!:  sum.t3  ■  01:  sum.t4  ■  0! 
sum.gl  =  0!:  sum.g2  *  0!:  sum.g3  *  0!:  sum.g4  *  0! 
sum.bl  =  0!:  sum.b2  =  0!:  sum.b3  *  0!:  sum.b4  *  0! 

FOR  i  *  1  TO  num.l 

Ddta.t0)  =  temp2C)  *  tempi  (i):  Delta. g(i)  *  gmptu«2(l)  •  gmptuslO) 

Deita.bC)  =  baUeigy20)  -  baUeigylC):  NEXT  i 

FOR  i  *  1  TO  num.l 

sum.tl  *  sum.tl  +  ABS{Delta.tC)): 
sum.gl  *  sum.gl  +  ABS(Ddta.g0)): 
sum.bl  =  sum.bl  +  ABS(Deltt.b0)): 
sum.t2  *  sum.t2  +  ABS(Delta.t0))  /  templQ) 
sum.g2  »  sum.g2  +  ABS(Delta.g(i))  /  gmptuslO) 
sum.b2  s  sum.b2  +  ABS(Delta.bO))  /  ballergy  l(i) 
sum.t3  =  sum.t3  +  Delta.  tC)  A  2 
sum.g3  =  sum.g3  +  Delta.  gO) A  2 
sum.b3  =  sum.b3  +  Delta.bC)  A  2 
sum.t4  =  sum.t4  +  (Delta.t(i)  /  tempi  0)) A  2 
sum.g4  =  sum.g4  +  (Delta.gC)  /  gmptuslO)  A  2 
sum.b4  s  sum.b4  +  (Delta.bC)  /  balkigylCO)  A  2 
NEXT  i 


PRINT  #9,  _ 

USING"  ###.#  ###.#  #####.# 

sum.tl  /(num.l  - 1);  sum.gl  /  (num.l  -  1);  sum.bl  /(num.l  - 1) 

PRINT  #9. :  PRINT  #9.  TAB(21);  "MEAN  ABSOLUTE  PERCENT  DEVIATIONS" 

PRINT  #9.  TAB(14);  "Temp  (K)  Impetus  (J/g)  Ballergy  (J/g)" 

PRINT  #9.  USING  _ 

M  UU#  M  MM  MMMM  MM 

,nnn  |  ^ 

100  *  sum.t2  /  (num.l  - 1);  100  *  sum.g2  /  (num.l  - 1);  _ 

100  *  sum.b2  /  (num.l  - 1) 

PRINT  #9. :  PRINT  #9,  TAB(23);  "ROOT  MEAN  SQUARE  DEVIATIONS" 

PRINT  #9.  TAB(13);  "Temp  (K)  Impetus  (J/g)  Ballergy  (J/g)" 

PRINT  #9.  USING  _ 

v*  ^  at  ^  h. 

SQR(sum.t3  /  (num.l  - 1));  SQR(sum.g3  /  (num.l  - 1));  _ 

SQR(sum.b3  /  (num.l  -  1)) 

PRINT  #9, :  PRINT  #9,  TAB(20);  "ROOT  MEAN  SQUARE  PERCENT  DEVIATIONS" 
PRINT  #9,  TAB(13);  "Temp  (K)  Impetus  (J/g)  Ballergy  (J/g)" 

PRINT  #9,  USING  _ 

ft  MM  MM  ^  n, 

,WW  1W#nw  WHnr*ir  f  — 

100  *  SQR(sum.t4  /  (num.l  - 1));  100  *  SQR(sum.g4  /  (num.l  - 1));  _ 

100  *  SQR(sum.b4  /  (num.l  •  1)) 
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PRINT  #33.  _ 

USING "  ###.#  «##.#  #####.«  _ 

sum.tl  /  (num.l  - 1);  sum.gl  /  (num.l  - 1);  sum.bl  /  (num.l  - 1) 

PRINT  #33. :  PRINT  #33.  TAB(21);  "MEAN  ABSOLUTE  PERCENT  DEVIATIONS" 
PRINT  #33.  TAB(14);  "Temp  (K)  Impetus  <I/g)  Baltetgy  (J/g)" 


PRINT #33,  USING. 


#.## 


####«## 


K, 

.  . 


100  *  sum.t2  /  (num.l  -  1);  100  *  sum.g2  /  (num.l  -  1);  _ 

100  *  sum.b2  /  (num.l  - 1) 

PRINT  #33. :  PRINT  #33,  TAB(23);  "ROOT  MEAN  SQUARE  DEVIATIONS" 

PRINT  #33.  TAB(13);  "Temp  (K)  Impetus  (J/g)  Ballergy  (J/g)" 


PRINT  #33.  USING  _ 

M  ^  u  na 

SQR(sum.t3  /  (num.l  - 1));  SQR(sum.g3  /  (num.l  - 1));  _ 

SQR(sum.b3  /  (num.l  - 1)) 

PRINT  #33, :  PRINT  #33.  TAB(20);  "ROOT  MEAN  SQUARE  PERCENT  DEVIATIONS" 
PRINT  #33.  TAB(13);  "Temp  (K)  Impetus  (J/g)  Balleigy  (J/g)" 

PRINT  #33,  _ 

USING"  ##.##  ##.##  ####.# 

100  *  SQR(sum.t4  /  (num.l  -  1));  100  *  SQR(sum.g4  /  (num.l  -  1));  _ 

100  *  SQR(sum.b4  /  (num.l  -  1)) 


PRINT  #33. :  PRINT  #33. :  PRINT  #33. 
PRINT  #253, .  PRINT  #253, :  PRINT  #253, 


COLOR  15. 0 

CLOSE  #1:  CLOSE  #2:  CLOSE  #33 
KILL  "TEMPOl":  KILL  "TEMP02" 
SYSTEM 


geLencrgy:  FOR  J  =  3  TO  LEN(lyne$):  SENTS  =  MID$(lyne$,  J.  1) 

IF  (SENTS  o  "  ")  THEN  GOTO  CONTT3 
NEXT:  PRINT  "TROUBLE  in  Subroutine  GETENERGY.  LYNES  =  lyneS 
STOP 

CONTT3:  lien  =  LEN(lyne$) 

FOR  11  =  lien  TO  1  STEP -1 
IF  (MIDS(lyneS,  U,  1)  =  "D  THEN  GOTO  contt4 
NEXT  11 

contt4:  IF  (MID$(lyne$,  U  - 1. 4)  s  "kJ/g")  THEN  GOTO  contt5 

CLS  :  PRINT  "RATS!  Problem  at ";  CHRS(34);  "  contt4";  CHR$(34):  STOP 
conttS:  Hj  =  11 

FOR  11  =  ll.j  TO  1  STEP  -1 

IF  (MID$(lyne$,  11. 1)  =  "+")  THEN  GOTO  contt6 

NEXT  n 

PRINT  "RATS!  Problem  at ";  CHR$(34);  "  contt5";  CHR$(34):  STOP 
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oonil6:  energy  -  VAL(MID$Oyne$,  U  +  1,  U.j  -  U  +  1)) 

PRINT  "Energy  «  energy 
IF  (energy  «  energy  jave)  THEN  GOTO  return 
energy  Jive  *  energy 
no.ener  «  noxner  +  1 
e  jeto(no  .ener) »  energy 

print  "Energy  =  energy 

return:  RETURN 

STAT.SUM: 

FOR  jj «  1  TO  no.ener  ii  *  jj  +  no.cmpens  *  (ixner  -  1) 
annul  «  sum.il  +  ABS(DelU.t(ii)) 
sum.t2  *  sunu2  +  ABS(DelU.t(U))  /  templ(li) 
sum.gl  »  sum.gl  +  ABS(Delu.gOi)) 
sum.g2  «  sum.g2  +  ABS(DeUa.gfii))  /  gmptuslQi) 
sum.bl  »  sum.bl  +  ABS(Delta.b(ii)) 
sum.b2  -  sum.b2  +  ABS(Delta.b(li))  /  baUergyl(ii) 

NEXT  jj 

*  PRINT  USING  _ 

M  J  MM  MMM  **. 

100  *  sum.t2  /  num.l;  100  *  sum.g2  /  num.l;  100  *  sum.b2  / 

RETURN 


'  A  subroutine  that  tries  to  determine  which  thermo  library  was  used 
libjuune: 

OPEN  filel.out$  FOR  INPUT  AS  #2: 
reads  1:  LINE  INPUT  #2.  lyneS 

IF  (LEFT$(lyne$.  40)  =  "  SPECIES  BEING  CONSIDERED  IN  THIS  SYSTEM")  THEN 
libjutme.l$  =  "McBride’s  standard  (CET89)  library" 

GOTO  close2 
END  IF 

FOR  jk  -  1  TO  LENOyneS) 

IF  (MID$(lyne$,  jk.  3)  o  "***")  THEN  GOTO  nexxtl 
libJiame.lS  =  RIGHT$(lyne$,  81  -  jk):  GOTO  close2 
nexxtl:  NEXT:  GOTO  readsl 
close2:  CLOSE  #2 

OPEN  file2.out$  FOR  INPUT  AS  #2 
reads3:  LINE  INPUT  #2,  lyne$ 

IF  (LEFTSOyneS.  40)  *  "SPECIES  BEING  CONSIDERED  IN  THIS  SYSTEM")  THEN 
libjuune 28  *  "McBride’s  standard  (CET89)  library" 

GOTO  dose2a 
END  IF 

FOR  jk-  1  TO  LEN(lyneS) 

IF  (MED$0yne$,  jk.  3)  o  "***")  THEN  GOTO  nexxt2 
libjuune  2S  «  RIGHTJOyneJ.  81  -  jk):  GOTO  close2a 
nexxt2:  NEXT:  GOTO  reads3 
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cloae2a:  CLOSE  #2:  RETURN 


m 


'  Check  prints 

*  GLS  :  PRINT  filel.out$,  fik2.out$ 

*  PRINT  line.coumer.1,  Une.coumer.2 

’  PRINT  libjuune.l$; "  and  Mttxname^S 

CLOSE  #2 
RETURN 

*  A  subroutine  that  computes  statistics  for  each  energy 
STAT: 

sum.tl  «  0!:  sum.t2  «  0!:  sum.gl  =  0! 
sum.g2  =  0!:  ijnax.b  =  1:  i.max.t  =  1 

sum.bl  *  0!:  sum.b2  =  0!:  (LpctbalLmax  at  0:  d.pcttempjnax  at  0 

FOR  nki  =  1  TO  no.cmpsns:  ii  *  nld  +  no.cmpsns  *  (Lener  - 1) 

Delta.g(ii)  »  gmptus2(ii)  -  gmptusl(ii) 

Delta.pct.tOi)  =  (temp2(ii)  -  templ(ii))  /  templOi) 

Deka.pctb(ii)  -  (baUagy2fii)  -  ballergyl(ii))  /  ballergyl(ii) 

IF  (ABS(Delu.pcLbOi))  >  d.pctballjnax)  THEN 
d.pcLbaHmax  *  ABS(Delta.pcLb(ii)):  ianax.b  =  nki 
d.balljnax  =  ballergy2(ii)  -  baUergylOi) 

END  IF 

IF  (ABS(Ddta.pcLtOi))  >  d.pcttempjnax)  THEN 
d.pct temp.max  =  ABS(Delta.pctt(ii)):  Lmax.t  =  nld 
temp.max  =  tempi  (ii) 

END  IF 
NEXT  nki 

PRINT  #9.  TAB(13); 

IF  (100!  *  d.pct  temp.max  >=  .1)  THEN 
PRINT  #9,  USING  _ 

"Max  abs  pet  diff  in  TEMPERATURE  =  ##.#%  for  case  no.  ##";  _ 

100!  *  d.pcttempjnax;  ijnax.t 

ELSE 

PRINT  #9,  USING  _ 

"Max  abs  pet  diff  in  TEMPERATURE  =  #.##%  for  case  no.  ##";  _ 

100!  *  d.pcttempjnax;  ijnax.t 

END  IF 

PRINT  #33.  TAB(13); 

PRINT  #33,  USING  "Max  abs  pet  diff  in  TEMPERATURE  =  ##.#%  for  case  no.  ##"; . 

100!  *  d.pct temp.max;  Lmax.t 


82 


PRINT  #253.  TAB(3); 

PRINT  #253.  USING  "Max  abc  pet  diff  in  TEMPERATURE  «  ##.#%  for  eve  no.  ##";  _ 

100!  *  (Lpet  temp.max;  Lsuuu 


PRINT  #9,  TAB(13X 

IF  (100!  *  (LpcLballmax  »  .1)  THEN 

PRINT  #9,  USING  "Max  ate  pet  diff  in  BALLERGY  -##.#%  for  case  no.  «T;  _ 

100!  *  (Lpet  ball  .max;  imax.b 

ELSE 

PRINT  #9.  USING  "Max  aba  pet  diff  in  BALLERGY  «#.##%  for  case  no.  ##";  _ 

100!  *  (Lpet  ball  .max;  imax.b 

END  IF 

PRINT  #253,  TAB(3); 

PRINT  #253,  USING  "Max  aba  pet  diff  in  BALLERGY  «##.#%  for  caw  no.  ##";  _ 

100!  *  <Lpct.baH.max;  imax.b 


PRINT  #33.  TAB(13); 

PRINT  #33,  USING  "Max  abs  pet  diff  in  BALLERGY  »##.#%  for  case  no.  ##";  _ 

100!  *  d.pctballmax;  imax.b 

*  PRINT  TAB(24);  "Mean  Absolute  Deviations" 

’  PRINT  TAB(14);  Temp  Impetus  Balletgy" 

RETURN:  *  end  of  STAT 

DATA  "d^exmuce£TSTLBSNO.STD" :  *  Summary  table,  no  ions,  standard  library 
DATA  "d^xmucefCTSTLBSNO-EXT"  :  ’  Summary  table,  no  ions,  extended  library 
DATA  "d^extnuceKTSTLBONO.STD"  :  ’  Output  table,  no  ions,  standard  library 
DATA  "d^xmuceftTSnjJONOJDCT"  :  ’  Output  table,  no  ions,  extended  library 

DATA  "dfcxmueefvTSTLBSIO.STD" :  ’  Summary  table,  w/  ions,  standard  library 
DATA  "dAexmucefcTSTLBSIOJEXT"  :  ’  Summary  table,  w/  ions,  extended  library 
DATA  "dNexmucefvTSTLBOIO.STD"  :  *  Output  table,  w/  ions,  standard  library 
DATA  "d^cxmueefcTSTLBOIOJEXT"  :  *  Output  table,  w/  ions,  extended  library 

DATA  "d:V:xmucefvTSTLBSNOJEXT"  :  ’  Summary  table,  no  ions,  extended  library 
DATA  "d.\exmuceftTSTLBSIOi£XT"  :  *  Summary  table,  w/  ions,  extended  library 
DATA  "d:^xmu«ftTSTLBONO.EXT"  :  ’  Output  table,  no  ions,  extended  library 
DATA  "d^cxmucefTSTLBOIO.EXT"  :  *  Output  table,  w/  ions,  extended  library 
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Appendix  E. 


Output  from  COMP  ARE. BAS: 

Comparing  the  Results  from  the  Standard  (CET89)  Library  with 
Those  from  McBride's  Extended  Library — Ions  Excluded 

The  two  EXMUCET  output  files  are  TSTLBONO.STD  and  TSTLBONO.EXT 

The  input  (SUMMARY)  files  are  TSTLBSNO.STD  and  TSTLBSNO.EXT 

Thermo  library  used  for  first  set  was 

***  McBride’s  Standard  (CET89)  Library  *** 

Thermo  library  used  for  second  set  was 

***  B.  McBride’s  EXTENDED-RANGE  Library  *** 

Number  of  energies  =  21 

0  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1)  30172  3010.2  1080.6  1077.9  4563.5  45832  1.2368  12352 

2)  2552.1  2549.1  944.8  943.5  3729.1  3736.5  12534  12525 

3)  3401.5  3392.2  1139.7  1136.3  5157.6  5184.0  1.2210  12192 

4)  3820.8  3811.9  1176.9  1173.5  56862  5723.1  1.2070  12050 

5)  2492.9  2490.6  929.0  927.9  3583.1  3588.6  12593  12586 

6)  3748.7  3738.6  1202.4  11982  6217 5  6278.9  1.1934  1.1908 

7)  1220.6  1216.5  474.6  472.6  2870.8  2899.8  1.1653  1.1630 

8)  3588.9  3574.6  1388.7  1382.5  65342  6590.0  1.2125  12098 

Max  abs  pet  diff  in  TEMPERATURE  =  0.4%  for  case  no.  8 
Mart  abs  pet  diff  in  B ALLERGY  =  1.0%  for  case  no.  7 

J  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1)  3320.1  3311.4  1191.1  1187.7  5074.8  5099.6  1.2347  1  2329 

2)  2881.3  2877.0  1067.5  1065.7  4274.3  4285.1  12497  12487 

3)  3674.1  3664.1  1235.8  1232.0  5609.1  5640.6  1.2203  1.2184 

4)  4022.5  4014.0  1251.2  1247.7  5996.8  6038.1  1.2086  1.2066 

5)  2828.0  2824.6  1054.8  1053.3  4129.1  4137.5  1.2555  1  2546 

6)  3903.8  3893.8  1266.6  1262.1  6479.5  6545.9  1.1955  1.1928 

7)  1306.9  1302.5  534.0  531.7  3121.9  3155.1  1.1710  1.1685 

8)  3795.4  3781.1  1475.7  14692  6948.8  7012.6  1.2124  1  2095 

Max  abs  pet  diff  in  TEMPERATURE  =  0.4%  for  case  no.  8 
Max  abs  pet  diff  in  BALLERGY  =  1.1%  for  case  no.  7 
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1  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Btll(Ext)  Gam(Std)  Gam(Ext) 

1)  36012  3591.6  1295.7  1291.9  5549.3  5578.8  1.2335  12316 

2)  3198.8  3193.1  1186.4  1184.1  48012  4816.0  1.2471  12459 

3)  3919.4  3909.3  1325.7  1321.7  6013.0  6049.3  12205  12185 

4)  4203.7  4195.5  1322.5  1318.7  62702  6315.7  1.2109  1  2088 

5)  3152.5  3147.7  1177.0  1175.0  46592  4671.3  1.2526  12515 

6)  4050.1  4040.1  1330.1  1325.3  67232  6794 S  1.1978  1.1951 

7)  1392.1  1387.5  597.5  594.8  3364.9  3400.3  1.1776  1.1749 

8)  3983.1  3969.3  15582  1551.5  7323.1  7395.0  12128  1  2098 

Max  abs  pet  diff  in  TEMPERATURE  =  0.3%  for  case  no.  8 
Max  abs  pet  diff  in  BALLERGY  =  1.1%  for  case  no.  6 

1.5  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1)  3858.0  3848.1  1394.2  1390.1  5981.7  6015.6  1.2331  12311 

2)  35002  3493.2  1300.6  1297.8  5302.5  5321.0  1.2453  12439 

3)  4139.4  4129.7  1410.5  1406.3  6371.8  6413.0  12214  12193 

4)  4371.0  4363.1  1392.1  1388.1  6517.5  6567.7  1.2136  12114 

5)  3461.6  3455.7  1294.6  12922  5165.8  5181.3  1.2506  12494 

6)  4190.1  4180.2  1393.4  1388.3  6953.6  7029.8  1.2004  1.1975 

7)  1478.4  1473.6  665.6  662 JS  3606.6  3642.4  1.1846  1.1819 

8)  4154.9  4141.9  16372  1630.4  7663.0  7743.0  1.2136  12106 

Max  abs  pet  diff  in  TEMPERATURE  =  0.3%  for  case  no.  7 
Max  abs  pet  diff  in  BALLERGY  =  1.1%  for  case  no.  6 

2  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  BaU(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1)  4090.8  4081.4  1487.1  1482.8  6372.1  6410.5  1.2334  12313 

2)  37822  3774.5  1409.7  1406.5  5773.0  5794.3  1.2442  12427 

3)  4338.5  4329.4  1491.4  1487.0  6692.3  6738.9  1.2229  12207 

4)  4528.8  4521.1  1461.1  1457.0  6746.4  68022  1.2166  12142 

5)  3752.1  3745.4  1407.1  1404.3  5643.6  5661.6  1.2493  12480 

6)  4325.3  4315.5  1456.9  1451.5  7173.7  7255.3  1.2031  12001 

7)  1567.4  1562.4  738.6  735.3  3851.5  3886.6  1.1918  1.1892 

8)  4313.8  4301.6  1713.7  1706.7  7974.7  8062.8  1.2149  12117 

Max  abs  pa  diff  in  TEMPERATURE  =  0.3%  for  case  no.  7 
Max  abs  pa  diff  in  BALLERGY  =  1.1%  for  case  no.  6 
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2.5  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  BalKExt)  Gam(Std)  Gam(Ext) 

1)  4302.8  4294.0  1575.5  1571.2  6724.7  67682  1.2343  12321 

2)  4043.6  4035.6  1513.5  1510.0  62093  62333  12437  12422 

3)  4521.3  4512.9  1569.8  15652  6982.0  7035.1  1.2248  1  2225 

4)  4680.1  4672.8  1530.3  1525.9  6961.9  7024.8  1.2198  12172 

5)  4022.4  4015.3  1514.3  15112  60892  6109.3  12487  12474 

6)  4456.9  4447.4  1520.9  15152  7385.9  7473.5  12059  12027 

7)  1660.3  1655.4  817.0  813.3  4103.8  4137.5  1.1991  1.1966 

8)  4462.4  44512  1788.4  1781.3  8263.7  8360.1  12164  12131 

Max  abs  pet  diff  in  TEMPERATURE  =  0.3%  for  case  no.  7 
Max  abs  pet  diff  in  BALLERGY  =  12%  for  case  no.  6 

3  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam  (Ext) 

1)  4497.7  4489.8  1660.6  16562  7045.6  7095.0  12357  12334 

2)  4285.0  4277.2  1612.7  1609.0  6611.9  6638.8  12439  12424 

3)  4692.1  4684.3  1646.6  1641.9  7247.9  7308.8  12272  12246 

4)  4827.1  4820.0  1599.9  1595.4  71672  7239.3  12232  12204 

5)  4272.9  4265.9  1616.8  1613.6  6502.0  6524.3  1.2487  12473 

6)  4585.8  4576.4  1585.6  1579.6  75913  7686.1  1.2089  12055 

7)  1758.8  1754.0  900.8  897.0  4368.4  4400.1  12062  12039 

8)  4603.1  4592.7  1862.1  1854.9  85343  8639.8  1.2182  12147 

Max  abs  pet  diff  in  TEMPERATURE  =  0.3%  for  case  no.  7 
Max  abs  pet  diff  in  BALLERGY  =  12%  for  case  no.  6 

33  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  BalKExt)  Gam(Std)  Gam  (Ext) 

1)  46792  4672.3  1743.4  1739.1  7340.6  7397.4  12375  12351 

2)  4508.4  4501.0  1707.9  1704.1  6982.0  7013.1  1.2446  12430 

3)  4853.8  4846.7  1722.7  1717.9  7494.9  7565.8  1.2298  12271 

4)  4971.0  4964.2  16703  1665.7  73643  7448.1  12268  12236 

5)  45052  4498.6  1715.3  1711.9  68832  6908.8  1.2492  12478 

6)  47123  4703.3  1651.0  1644.8  7791.7  7894.3  12119  1  2084 

7)  1864.4  1859.8  990.3  986.4  4650.8  4680.3  12129  12108 

8)  4737.3  4727.8  1935.2  1927.9  8790.6  8905.3  1.220?  12165 

Max  abs  pa  diff  in  TEMPERATURE  =  02%  for  case  no.  7 
Max  abs  pa  diff  in  BALLERGY  =  1.3%  for  case  no.  6 
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4  kJ/g  added 


T(Sld) 

T(Ext)  Imp(Std)  Imp(Ext) 

BaU(Std)  Ball(Ext)  Gam(Std)  Gam  (Ext) 

1) 

4850.5 

4844.6 

1824.9 

1820.6 

76143 

7681.0 

1.2397 

12370 

2) 

4716.4 

4709.6 

1800.1 

1796.3 

7322.8 

7360.3 

1.2458 

12441 

3) 

50090 

5002.4 

1798.8 

1793.8 

7726.9 

7810.4 

12328 

12297 

4) 

5113.0 

5106.4 

1742.1 

1737.1 

7554.9 

76532 

12306 

12270 

5) 

4721.9 

4715.9 

1810.6 

18072 

72352 

7266.4 

1.2502 

12487 

6) 

4837.4 

4828.6 

1717.4 

1710.9 

79872 

80992 

1.2150 

12112 

7) 

1978.7 

1974.5 

1085.2 

1081.3 

4957.0 

4984.6 

12189 

12169 

8) 

4866.5 

4857.9 

2008.0 

2000.6 

9034.6 

9159.6 

1.2223 

12184 

Max  abs  pet  diff  in  TEMPERATURE  = 

0.2%  for  case  no. 

7 

Max  abs  pet  diff  in  BALLERGY  =  1.4%  for  case  no.  6 

4.5  kJ/g  added 

T(Std) 

T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball  (Ext)  Gam(Std)  Gam  (Ext) 

1) 

5014.0 

5008.9 

1905.8 

1901.4 

7871.0 

7949.9 

1.2421 

12392 

2) 

4911.6 

4905.5 

1890.2 

1886.4 

7637.8 

7684.6 

1.2475 

12455 

3) 

5159.3 

5153.3 

1875.1 

1870.0 

7946.4 

8045.7 

12360 

12324 

4) 

5253.7 

5247.4 

1815.0 

1809.7 

7739.0 

7855.8 

1.2345 

12304 

5) 

4925.4 

4920.0 

1903.5 

1900.1 

7560.7 

7600.9 

1.2518 

12500 

6) 

4961.0 

4952.5 

1784.9 

1778.1 

8178.4 

8301.4 

12182 

12142 

7) 

2102.4 

2099.0 

1184.9 

1181.3 

5291.1 

5318.0 

1.2239 

1.2221 

8) 

4991.6 

4983.8 

2080.9 

2073.4 

9268.4 

9404.8 

1.2245 

12205 

Max  abs  pet  diff  in  TEMPERATURE  = 

0.2%  for  case  no. 

6 

Max  abs  pet  diff  in  BALLERGY  =  13%  for  case  no.  4 

5  kJ/g  added 

T(Std) 

T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam  (Ext) 

1) 

5171.6 

5167.3 

1986.4 

1982.1 

8112.8 

8207.6 

1.2448 

12415 

2) 

50965 

5091.1 

1978.9 

1975.0 

7929.8 

7990.0 

1.2496 

1.2472 

3) 

5306 2 

5300.6 

1952.1 

1946.8 

8155.0 

8274.1 

1.2394 

12353 

4) 

5393.7 

5387.6 

1889.3 

1883.7 

7916.9 

8056.8 

1.2386 

12338 

5) 

5118.3 

5113.5 

1994.9 

1991.5 

7862.4 

7916.1 

1.2537 

12516 

6) 

5083 5 

5075.3 

1853.4 

1846.3 

8365.6 

8501.4 

1.2216 

12172 

7) 

22345 

2232.7 

1287.8 

1284.9 

5651.9 

5680.7 

1.2279 

12262 

8) 

5113.5 

5106.4 

2154.1 

21463 

94933 

9642.4 

1.2269 

12226 

Max  abs  pet  diff  in  TEMPERATURE  = 

0.2%  for  case  no. 

6 

Max  abs  pet  diff  in  BALLERGY  =  1.8%  for  case  no.  4 
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5.5  kJ/g  added 


T(Std) 

T(Ext)  Imp(Std)  Imp(Ext) 

Bail(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1) 

5324.7 

5321.1 

2067.4 

2063.0 

83415 

8456.6 

1.2478 

1.2440 

2) 

5273.3 

5268.5 

2066.8 

2062.8 

82015 

8280.0 

1.2520 

12491 

3) 

54505 

5445.2 

2030.0 

2024.4 

8353.6 

8497.3 

1.2430 

12382 

4) 

55335 

5527.4 

1965.1 

1959.1 

8088.4 

8256.8 

1.2430 

12373 

5) 

5302.6 

5298.3 

2085.3 

2081.8 

8142.4 

82155 

1.2561 

12534 

6) 

5205.1 

519752 

1923.0 

1915.7 

8548.9 

8699.7 

1.2249 

12202 

7) 

2371.8 

2372.8 

1391.8 

1390.3 

6029.4 

6064.7 

1.2308 

12292 

8) 

5232.7 

5226.2 

2227.8 

2220.0 

9710.9 

9873.8 

1.2294 

12248 

Max  abs  pet  diff  In  TEMPERATURE  =  0.2%  for  case  no.  6 
Max  abs  pet  diff  in  B ALLERGY  =  2.1%  for  case  no.  4 


6  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ex!)  Ball(Std)  BaU(Ext)  Gam(Std)  Gam(Ext) 


1) 

5474.4 

5471.3 

2148.8 

2144.3 

8558.1 

8698.6 

1.2511 

12465 

2) 

5443.7 

5439.3 

2154.4 

21502 

84545 

8557.6 

1.2548 

12513 

3) 

5593.1 

5588.0 

2109.0 

2103.1 

8542.7 

8716.4 

1.2469 

12413 

4) 

5673.3 

56672 

20425 

2036.1 

8252.8 

8456.3 

12475 

12408 

5) 

54802 

5476.2 

2175.4 

2171.7 

8402.3 

8502.0 

1.2589 

12554 

6) 

5326.1 

5318.5 

1993.9 

19862 

8728.3 

8896.4 

1.2284 

12233 

7) 

2510.3 

2515.2 

1495.7 

1495.9 

6410.6 

64565 

12333 

12317 

8) 

5349.6 

5343.9 

2302.0 

2294.1 

9921.4 

10099.8 

12320 

12271 

Max  abs  pet  diff  in  TEMPERATURE  =  0.2%  for  case  no.  7 
Max  abs  pet  diff  in  B ALLERGY  =  2.5%  for  case  no.  4 


6.5  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  BaU(Ext)  Gam(Std)  Gam(Ext) 


1) 

5621.6 

5618.8 

2231.0 

2226.3 

8763.0 

8935.3 

1.2546 

1.2492 

2) 

56092 

5605.0 

2242.1 

2237.7 

8690.0 

8825.1 

1.2580 

1.2536 

3) 

5734.6 

5729.6 

21892 

2182.9 

8722.1 

8932.6 

1.2510 

12444 

4) 

5813.6 

5807.2 

2121.6 

2114.6 

8409.3 

8655.7 

1.2523 

12443 

5) 

5652.6 

5648.7 

22652 

2261.4 

8642.8 

8778.0 

1.2621 

12576 

6) 

5446.6 

5439.4 

2066.0 

2058.0 

8903.7 

9091.8 

1.2320 

12264 

7) 

2648.4 

26572 

1599.3 

1601.0 

6789.0 

6846.4 

1.2356 

12338 

8) 

5464.8 

5459.7 

2376.9 

2368.8 

101252 

10321.1 

1.2347 

12295 

Max  abs  pet  diff  in  TEMPERATURE  =  0.3%  for  case  no.  7 
Max  abs  pa  diff  in  B ALLERGY  =  2.9%  for  case  no.  4 
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7  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam  (Ext) 

1)  5767.0  5764.4  2314.1  2309.1  8955.9  9167.5  1.2584  1.2519 

2)  5771.0  5766.8  2330.3  2325 J  8908.2  9084.4  1.2616  12560 

3)  5875 5  5870.3  2270.9  2264.0  8891.1  9146.3  1.2554  12475 

4)  5954.5  5947.7  2202.5  2194.8  8556.4  88552  1.2574  1  2479 

5)  5821.0  5816.9  2355.9  2351.4  8863.8  9045.3  1.2658  1  2600 

6)  5566.8  5559.8  2139.4  2131.0  9074.8  9286.0  1.2358  1  2295 

7)  2786.6  2798.7  1703.3  1706.3  7166.0  7233.1  1.2377  1  2359 

8)  5578.6  5574.0  2452.6  2444.4  10323.6  10538.4  1.2376  1.2320 

Max  abs  pet  diff  in  TEMPERATURE  =  0.4%  for  case  no.  7 

Max  abs  pet  diff  in  B ALLERGY  =  3.5%  for  case  no.  4 

7  J  lcJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1)  59112  5908.4  2398.3  2392.8  91362  9396.1  1.2625  12547 

2)  5930.1  5925.5  2419.3  2413.8  9108.9  9336.9  1.2656  12585 

3)  60162  6010.5  2354.0  2346.4  9048.8  93582  1.2601  12507 

4)  6096.3  6088.7  22852  2276.6  8692.6  9054.8  1.2629  12514 

5)  5986.6  5981.7  2447.0  2441.7  9064.6  9305.7  1.2700  12624 

6)  5686.8  5680.0  22142  2205.4  9241.3  9479.0  1.2396  12327 

7)  29262  2940.7  1808.4  1812.4  7546.9  7620.6  1.2396  12378 

8)  5691.1  5687.1  2529.1  2520.7  10515.8  10751.9  1.2405  1.2344 

Max  abs  pet  diff  in  TEMPERATURE  =  0.5%  for  case  no.  7 

Max  abs  pet  diff  in  B ALLERGY  =  42%  for  case  no.  4 

8  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  ImpCExt)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1)  6054.5  6051.3  2483.7  2477.5  9302.9  9621.7  1.2670  12575 

2)  60872  6081.8  2509.1  2502.8  9291.0  9584.0  1.2701  12611 

3)  6157.0  6150.6  2438.7  2430.1  9193.9  9568.6  1.2653  12540 

4)  6239.3  6230.5  2369.9  2360.1  8815.9  9254.8  1.2688  12550 

5)  6150.0  6144.0  2538.9  2532.6  9243.8  95602  1.2747  12649 

6)  5806.8  5800.1  2290.3  2281.1  9402.8  9671.0  1.2436  12359 

7)  3068.3  3084.2  1914.9  1919.6  79362  8013.3  1.2413  12396 

8)  5802.8  5799.1  2606.5  2597.9  107022  10962.1  1.2435  1.2370 

Max  abs  pet  diff  in  TEMPERATURE  =  0.5%  for  case  no.  7 
Max  abs  pet  diff  in  B ALLERGY  =  5.0%  for  case  no.  4 
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8.5  kJ/g  added 

T(Std)  T(Ext)  ImjKStd)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam  (Ext) 

1)  6197.5  6193.4  2570.4  2563.3  9454.3  9844.8  1.2719  1.2604 

2)  6243.1  6236.4  2600.1  2592.6  9452.8  9826.4  1.2751  1.2638 

3)  6298.3  6290.7  2525.1  2515.4  9324.5  9777.7  1.2708  1.2573 

4)  6383.6  6373.1  2456.5  2445.3  8923.6  9455.2  1.2753  1.2586 

5)  6312.1  63042  2631.9  2624.1  9399.2  9809.8  1.2800  1.2675 

6)  5926.9  59202  2367.9  2358.1  9558.7  9861.8  1.2477  1.2391 

7)  3213.3  3229.8  2022.8  2027.8  8336.6  8414.3  1.2426  12410 

8)  5913.7  5910.4  2685.0  2676.1  10882.7  111692  12467  1.2396 

Max  abs  pci  diff  in  TEMPERATURE  =  0.5%  for  case  no.  7 
Max  abs  pet  diff  in  B ALLERGY  =  6.0%  for  case  no.  4 

9kJ/gadded 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1)  6340.4  6334.9  2658.5  26502  9588.5  10065.7  1  2773  12633 

2)  63982  6389.6  2692.3  2683.3  95922  10065.0  12807  12666 

3)  6440.4  6431.1  26132  2602.0  9438.4  9985.8  12769  12606 

4)  6529.4  6516.7  2545.3  25322  9012.9  9656.0  1.2824  12622 

5)  6473.4  6462.8  2726.0  2716.3  9527.6  10055.4  12861  12701 

6)  6047.1  6040.4  2447.0  2436.6  9708.4  10051.4  12520  12424 

7)  3360.8  3377.3  2131.6  2136.8  87482  8824.7  1.2437  12421 

8)  6024.0  6021.1  2764.4  27552  11057.1  11373.2  12500  1.2423 

Max  abs  pet  diff  in  TEMPERATURE  «  0.5%  for  case  no.  7 
Max  abs  pet  diff  in  BALLERGY  =  7.1%  for  case  no.  4 

9.5  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1)  6483.5  6476.1  2748.1  2738.1  9703.1  10284.6  12832  12662 

2)  65532  6541.8  2785.9  2774.8  9705.9  10300.2  12870  12694 

3)  6583.4  6572.0  27032  26902  9532.8  10193.0  12836  12639 

4)  6677.0  6661.3  2636.3  2620.8  90802  98572  12903  12659 

5)  6634.3  6620.1  2821.5  28092  9625.3  10297.4  12931  12728 

6)  6167.7  6160.8  2527.5  2516.5  9851.1  10239.9  12566  1  2458 

7)  3510.0  3526.1  2240.8  2246.1  9169.4  9243.4  1  2444  12430 

8)  6134.0  6131.3  2844.9  2835.3  11225.1  11574.3  12534  12450 

Max  abs  pa  diff  in  TEMPERATURE  =  0.5%  for  case  no.  7 
Max  abs  pa  diff  in  BALLERGY  «  8.6%  for  case  no.  4 
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10  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  ImpOExt)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1)  6627.1  6617.0  2839.1  2827.2  9795.0  10501.7  1.2899  1.2692 

2)  67082  6693.2  2881.0  2867.3  9790.4  10532.5  1.2943  1.2722 

3)  6727.8  6713.3  2795.2  2779.7  9604.3  10399.3  1.2910  12673 

4)  6826.7  6807.0  2729.6  2711.1  9121.1  10058.8  1  2993  12695 

5)  6795.3  67'  2918.4  2902.8  9687.7  10536.4  1.3012  1  2755 

6)  6288.7  62  2609.6  2597.7  9985.9  10427.0  1  2613  1.2491 

7)  3659.8  36/j.J  2349.8  2355.0  9597.0  9667.9  1.2448  1  2436 

8)  6243.7  6241.1  2926.5  2916.5  11386.3  11772.6  1.2570  1.2477 

Max  abs  pet  diff  in  TEMPERATURE  =  0.4%  for  case  no.  7 
Max  abs  pet  diff  in  B ALLERGY  =  10.3%  for  case  no.  4 

GRAND  TOTALS 

MEAN  ABSOLUTE  DEVIATIONS 
Temp  (K)  Impetus  (J/g)  Balleigy  (J/g) 

7.5  6.0  177.0 

MEAN  ABSOLU  i  a  PTv  LENT  DEVIATIONS 
Temp  (K)  Impetus  (J/g  Balleigy  (J/g) 

0.175  0.31  2.QI 

ROOT  MEAN  SQUARE  DEVIATIONS 
Temp  (K)  Impetus  (J/g)  Balleigy  (j/g) 

8.3  6.7  262.1 

ROOT  MEAN  SQUARE  PERCENT  DEVIATIONS 
Temp  (K)  Impetus  (J/g)  Balleigy  (J/g) 

0.20  0.33  2.8 
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Appendix  F. 


Output  from  COMPARE. BAS. 

Comparing  die  Results  from  the  Standard  (CET89)  Library  with 
Those  from  McBride's  Extended  Library — Ions  Included 

The  two  EXMUCET  output  files  are  TSTLBOIO.STD  and  TSTLBOIO.EXT 

The  input  (SUMMARY)  files  are  TSTLBSIO.STD  and  TSTLBSIO.EXT 

Thermo  library  used  for  first  set  was 

***  McBride’s  Standard  (CET89)  Library  *** 

Thermo  library  used  for  second  set  was 

**♦  B.  McBride’s  EXTENDED-RANGE  Library  *** 


Number  of  energies  =  21 


OkJ/g  added 

T(Std)  T(Ext)  Imp(Std)  lmp(Ext) 

1)  30172  3010.2  1080.6  1077.9 

2)  2552.1  2549.1  944.8  943.5 

3)  3401.5  3392.2  1139.7  1136.3 

4)  3820.8  3811.9  1176.9  1173.5 

5)  2492.9  2490.6  929.0  927.9 

6)  3748.7  3738.6  1202.4  11982 

7)  1220.6  1216.5  474.6  472.6 

8)  3588.9  3574.6  1388.7  1382.5 


BaD(Std)  BaU(Ext) 

4563.5  45832 

3729.1  3736.5 

5157.6  5184.0 
56862  5723.1 

3583.1  3588.6 
6217.5  6278.9 
2870.8  2899.8 
65342  6590.1 


Gam(Std)  Gam(Ext) 
1.2368  12352 
12534  12525 
1.2210  12192 
1.2070  12050 
12593  12586 
1.1934  1.1908 
1.1653  1.1630 
1.2125  12098 


Max  abs  pet  diff  in  TEMPERATURE  =  0.4%  for  case  no.  8 
Max  abs  pet  diff  in  B ALLERGY  =  1.0%  for  case  no.  7 


5  kJ/g  added 

T(Std)  T(Ext)  lmp(Std)  Imp(Ext)  BaU(Std)  Ball(Ext)  Gam(Std)  Gam  (Ext) 


1) 

3320.1 

3311.4 

1191.1 

1187.7 

5074.8 

5099.6 

1.2347 

12329 

2)  2881.3 

2877.0 

1067.5 

1065.7 

4274.3 

42852 

1.2497 

12487 

3) 

3674.1 

3664.1 

1235.8 

1232.0 

5609.1 

5640.6 

1.2203 

12184 

4)  4022.5 

4014.0 

1251.2 

1247.7 

5996.8 

6038.1 

1.2086 

12066 

5) 

2828.0 

2824.6 

1054.8 

1053.3 

4129.1 

4137.6 

1.2555 

12546 

6) 

3903.8 

3893.8 

1266.6 

1262.1 

6479.5 

6545.9 

1.1955 

1.1928 

7) 

1306.9 

1302.5 

534.0 

531.7 

3121.9 

3155.1 

1.1710 

1.1685 

8) 

3795.4 

3781.1 

1475.7 

14692 

6948.8 

7012.6 

1.2124 

12095 

Max  abs  pet  diff  in  TEMPERATURE  =  0.4%  for  case  no.  8 
Max  abs  pet  diff  in  B ALLERGY  =  1.1%  for  case  no.  7 
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1  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1)  3601.2  3591.6  1295.7  1291.9  5549.3  5578.8  1.2335  12316 

2)  3198.8  3193.1  1186.4  1184.1  4801.2  4816.0  1.2471  1.2459 

3)  3919.4  3909.3  1325.7  1321.7  6013.0  6049.3  1.2205  1.2185 

4)  4203.7  4195.5  1322.4  1318.7  62702  6315.8  1.2109  1.2088 

5)  31525  3147.7  1177.0  1175.0  4659.2  4671.3  1.2526  12515 

6)  4050.1  4040.1  1330.1  1325.3  67232  6794.5  1.1978  1.1951 

7)  1392.1  1387.5  597.5  594.8  3364.9  3400.3  1.1776  1.1749 

8)  3983.1  3969.3  15582  1551-5  7323.1  7395.0  12128  12098 

Max  ahs  pet  diff  In  TEMPERATURE  =  0.3%  for  case  no.  8 
Max  abs  pet  diff  in  B ALLERGY  «  1.1%  for  case  no.  6 

15  kJ/g  added 

T(Std)  T(Ext)  ImpCStd)  Imp(Ext)  BaB(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1)  3858.0  3848.1  13942  1390.1  5981.7  6015.7  1  2331  12311 

2)  35002  3493.2  1300.6  1297.8  53025  5321.0  1.2453  12439 

3)  4139.4  4129.7  14105  1406.3  6371.8  6413.0  12214  12193 

4)  4371.0  4363.1  1392.1  1388.1  65175  6567.7  1.2136  12114 

5)  3461.6  3455.7  1294.6  12922  5165.8  5181.3  1.2506  12494 

6)  4190.1  41802  1393.4  1388.3  69535  7029.8  1.2004  1.1975 

7)  1478.4  1473.6  665.6  6625  3606.6  3642.4  1.1846  1.1819 

8)  4154.9  4141.9  16372  1630.4  7663.0  7743.0  1.2136  12106 

Max  abs  pet  diff  in  TEMPERATURE  =  0.3%  for  case  no.  7 
Max  abs  pet  diff  in  BALLERGY  =  1.1%  for  case  no.  6 

2  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam  (Ext) 

1)  4090.8  4081.4  1487.1  1482.8  6372.1  6410.6  1  2334  12313 

2)  37822  3774.5  1409.7  14065  5773.0  5794.4  1.2442  12427 

3)  43385  4329.4  1491.4  1487.0  66922  6738.9  1.2229  1  2207 

4)  4528.8  4521.1  1461.1  1457.0  6746.4  6802.3  1.2166  12142 

5)  3752.1  3745.4  1407.1  1404.3  5643.6  5661.6  1.2493  1  2480 

6)  4325.3  43155  1456.9  14515  7173.7  7255.3  1.2031  12001 

7)  1567.4  1562.4  738.6  735.3  38515  3886.6  1.1918  1.1892 

8)  4313.8  4301.6  1713.7  1706.7  7974.7  8062.8  12149  12117 

Max  abs  pet  diff  in  TEMPERATURE  =  0.3%  for  case  no.  7 
Max  abs  pet  diff  in  BALLERGY  =  1.1%  for  case  no.  6 
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2JkJ/fe  added 

T(Std)  T(Ext)  Imp(Std)  ImpCExt)  Ball(Std)  Ball(Ext)  Gam(Sld)  Gam(Ext) 

1)  4302.8  4294.0  1575.5  1570  6724.7  6768.2  1.2343  0321 

2)  4043.6  4035.6  1513.5  1510.0  6209.5  6233.5  0437  0422 

3)  4521.3  4512.9  1569.8  1565.2  6982.0  7035.1  1.2248  1.2225 

4)  4680.1  4672.8  1530.3  1525.9  6961.9  7024.8  1.2198  0172 

5)  4022.4  4015.3  1514.3  1510  6089.2  6109.3  1.2487  1.2474 

6)  4456.9  4447.4  1520.9  1515.2  7385.9  7473.5  1.2059  12027 

7)  1660.3  1655.4  817.0  813.3  4103.8  4137.5  1.1991  1.1966 

8)  4462.4  4450  1788.4  1781.3  8263.7  8360.1  0164  0131 

Max  abs  pet  diff  in  TEMPERATURE  *  0.3%  for  case  no.  7 
Max  abs  pet  diff  in  B ALLERGY  =  0%  for  case  no.  6 

3kJ/gadded 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1)  4497.7  4489.8  1660.6  16562  7045.6  7095.0  12357  0334 

2)  4285.0  42772  1612.7  1609.0  6611.8  6638.9  12439  12424 

3)  4692.0  4684.3  1646.6  1641.9  7247.8  7308.8  12272  12246 

4)  4827.0  4820.0  1599.9  1595.4  71672  7239.3  12232  12204 

5)  4272.9  4265.9  1616.8  1613.6  6502.0  6524.3  1.2487  0473 

6)  4585.8  4576.4  1585.6  1579.6  7590  7686.1  1.2089  0055 

7)  1758.8  1754.0  900.8  897.0  4368.4  4400.1  12062  12039 

8)  4603.0  4592.7  1862.1  1854.9  8534.5  8639.8  0182  0147 

Max  abs  pet  diff  in  TEMPERATURE  *=  0.3%  for  case  no.  7 
Max  abs  pet  diff  in  BALLERGY  *  12%  for  case  no.  6 

3J  kj/g  added 

T(Std)  T(Ext)  Imp(Std)  ImpQExt)  BaD(Std)  BalKExt)  Gam(Std)  Gam(Ext) 

1)  46792  4672.3  1743.4  1739.1  7340.5  7397.4  12375  0351 

2)  4508.4  4501.0  1707.9  1704.1  6981.9  7013.1  1.2446  12430 

3)  4853.8  4846.7  1722.7  1717.9  7494.9  7565.8  1.2298  12271 

4)  4971.0  49642  1670.5  1665.7  7364.4  7448.1  1.2268  12236 

5)  45052  4498.6  1715.3  1711.9  68832  6908.8  12492  12478 

6)  4712.4  4703.3  1651.0  1644.8  7791.7  7894.3  0119  12084 

7)  1864.4  1859.8  990.3  986.4  4650.8  4680.3  0129  0108 

8)  4737 3  4727.8  19352  1927.9  8790.6  8905.3  12201  0165 

Max  abs  pet  diff  in  TEMPERATURE  *  02%  for  case  no.  7 
Max  abs  pet  diff  in  BALLERGY  =  1.3%  for  case  no.  6 
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4  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Exi)  Cam(Std)  Gam  (Ext) 

1)  4850.5  4844.6  1824.9  1820.6  7614.5  7680.9  1.2397  12370 

2)  4716.3  4709.6  1800.1  1796.3  7322.8  7360.3  1.2458  12441 

3)  5008.9  5002.4  1798.7  1793.8  7726.8  7810.4  1.2328  12297 

4)  5112.9  5106.4  1742.1  1737.1  7554.8  7653.1  1.2306  12270 

5)  4721.9  4715.9  1810.5  18072  7235.1  7266.4  1.2502  1  2487 

6)  4837.4  4828.6  1717.4  1710.9  79872  80992  1.2150  12112 

7)  1978.7  1974.5  10852  1081.3  4957.0  4984.6  12189  12169 

8)  4866.5  4857.8  2008.0  2000.6  9034.6  9159.6  1.2223  12184  . 

Max  abs  pet  diff  in  TEMPERATURE  =  02%  for  case  no.  7 
Max  abs  pet  diff  in  B ALLERGY  =  1.4%  for  case  no.  6 

4.5  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1)  5014.0  5008.9  1905.7  1901.4  7871.0  7949.9  1.2421  12392 

2)  4911.5  4905.5  1890.2  1886.4  7637.7  7684.6  1.2475  12455 

3)  5159.3  5153.3  1875.1  1870.0  7946.3  8045.7  12360  12324 

4)  5253.6  5247.3  1814.9  1809.7  7738.9  7855.7  1.2345  1  2304 

5)  4925.4  4920.0  1903.5  1900.1  7560.6  7600.8  1.2518  12500 

6)  4960.9  4952.4  1784.8  1778.1  8178.3  8301.4  12182  12142 

7)  2102.4  2099.0  1184.9  1181.3  5291.1  5318.0  1.2239  12221 

8)  4991.6  4983.7  2080.9  2073.4  9268.4  9404.7  1.2245  12205 

Max  abs  pet  diff  in  TEMPERATURE  =  02%  for  case  no.  6 
Max  abs  pet  diff  in  BALLERGY  =  1.5%  for  case  no.  4 

5  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  BaD(Std)  Ball(Ext)  Gam(Std)  Gam  (Ext) 

1)  5171.6  5167.3  1986.4  1982.1  8112.7  8207.6  1.2449  12415 

2)  5096.4  5091.1  1978.9  1975.0  7929.8  7990.0  1.24%  12472 

3)  5306.1  5300.5  1952.1  1946.7  8154.9  8274.0  1.2394  1  2353 

4)  5393.6  5387.5  18892  1883.6  7916.8  8056.6  1.2386  1  2338 

5)  51182  5113.5  1994.9  1991.4  7862.3  7916.1  1.2537  12516 

6)  5083.4  507 52  1853.3  1846.3  8365.5  8501.4  12215  12172 

7)  2234.5  2232.7  1287.8  1284.9  5651.9  5680.7  1.2279  12262 

8)  5113.4  5106.3  2154.1  2146.5  9493.4  9642.4  1.2269  12226 

Max  abs  pet  diff  in  TEMPERATURE  =  0.2%  for  case  no.  6 
Max  abs  pet  diff  in  BALLERGY  =  1.8%  for  case  no.  4 


% 


5.5  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1)  5324.6  5321.0  2067.3  2062.9  8341.4  84565  1.2478  12439 

2)  52732  5268.4  2066.7  2062.8  8201.4  8280.0  1.2520  12491 

3)  5450.4  5445.2  2030.0  2024.4  8353.5  84972  1.2430  12382 

4)  5533.3  5527.3  1965.0  1959.1  80882  8256.6  1.2429  12373 

5)  5302.5  5298.3  2085.3  2081.8  8142.3  8215.5  1.2561  12534 

6)  5205.0  51972  1923.0  1915.6  8548.8  8699.6  12249  12202 

7)  2371.8  2372.8  1391.8  1390.3  6029.4  6064.7  1  2308  12292 

8)  5232.6  52262  2227.7  2220.0  9710.8  9873.7  12294  12248 

Max  abs  pet  diff  in  TEMPERATURE  =  0.1%  for  case  no.  6 
Max  abs  pet  diff  in  B  ALLERGY  x  2.1%  for  case  no.  4 

6kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1)  5474.3  54712  2148.8  21442  8558.0  8698.5  12511  12465 

2)  5443.6  5439.2  2154.3  21502  8454.4  85575  1.2548  12513 

3)  5593.0  5587.9  2108.9  2103.0  85425  8716.3  1.2469  12413 

4)  56732  5667.1  2042.4  2036.0  8252.7  84562  12475  12408 

5)  5480.1  5476.2  2175.3  2171.7  84022  8501.9  1.2589  12554 

6)  5326.0  5318.5  1993.8  19862  8728.1  8896.3  1.2284  12233 

7)  2510.3  2515.2  1495.7  1495.9  6410.6  64565  1.2333  12317 

8)  5349.6  5343.8  2302.0  2294.1  9921.3  10099.6  12320  12271 

Max  abs  pet  diff  in  TEMPERATURE  =  0.2%  for  case  no.  7 
Max  abs  pet  diff  in  B ALLERGY  =  25%  for  case  no.  4 

65  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam  (Ext) 

1)  56215  5618.7  2230.9  22262  8762.8  8935.1  12546  12492 

2)  5609.1  5605  0  2242.1  2237.7  8689.8  8825.0  1.2580  12536 

3)  57345  5729.4  2189.2  2182.8  8721.9  8932.4  1.2510  12444 

4)  5813.4  5807.1  21215  21145  8409.1  86555  1.2523  12443 

5)  5652.4  5648.6  2265.4  2261.4  8642.6  8777.8  1.2621  12576 

6)  54465  5439.3  2065.9  2057.9  89035  9091.6  1.2320  12264 

7)  2648.4  26572  1599.3  1601.0  6789.0  6846.4  1.2356  12338 

8)  5464.7  5459.6  2376.8  2368.8  10125.4  10320.9  12347  1.2295 

Max  abs  pa  diff  in  TEMPERATURE  =  0.3%  for  case  no.  7 
Max  abs  pa  diff  in  B ALLERGY  =  2.9%  for  case  no.  4 
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7  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  BaU(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1)  5766.9  5764.2  2314.0  2309.0  8955.7  9167.3  1.2584  1.2519 

2)  5770.8  5766.7  2330.2  2325.4  8908.0  9084.2  1.2616  1.2560 

3)  5875.3  5870.1  2270.7  2263.9  8890.8  9146.1  1.2554  1.2475 

4)  5954.3  5947.5  2202.3  2194.7  8556.2  8854.9  1.2574  1.2479 

5)  5820.8  5816.8  2355.8  2351.3  8863.6  90452  1.2658  12600 

6)  5566.7  5559.7  2139.3  2131.0  9074.6  9285.8  1.2357  12295 

7)  2786.6  2798.7  1703.3  1706.3  7166.0  72332  1.2377  12359 

8)  5578.5  5573.9  2452.5  2444.3  10323.4  10538.1  12376  12319 

Max  abs  pet  diff  in  TEMPERATURE  =  0.4%  for  case  no.  7 

Max  abs  pet  diff  in  B ALLERGY  =  35%  for  case  no.  4 

7.5  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1)  5911.0  5908.2  23982  2392.7  9136.0  9395.8  1.2625  12547 

2)  5929.8  5925.3  2419.1  2413.7  9108.6  9336.7  1.2656  12585 

3)  6015.9  6010.3  2353.8  2346.3  9048.6  9357.9  1.2601  12507 

4)  6096.0  6088.5  2285.0  22765  8692.3  90545  1.2629  12514 

5)  5986.3  5981.6  2446.8  2441.6  9064.4  9305.4  1.2699  1  2624 

6)  5686.7  5679.9  2214.1  2205.3  9241.1  9478.8  1.2396  12327 

7)  29262  2940.7  1808.4  1812.4  7546.9  7620.6  1.2396  12378 

8)  5691.0  5686.9  2529.0  2520.6  10515.6  10751.7  1.2405  1.2344 

Max  abs  pa  diff  in  TEMPERATURE  =  0.5%  for  case  no.  7 

Max  abs  pa  diff  in  B ALLERGY  &  42%  for  case  no.  4 

8  kJ/g  added 

T(Std)  T(Ext)  Imp(Std)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1)  6054.3  6051.1  2483.6  2477.4  9302.6  9621.4  1.2670  12575 

2)  6086.9  6081.6  2509.0  2502.7  9290.7  9583.7  1.2701  12611 

3)  6156.7  6150.3  2438.5  2430.0  9193.6  95682  1.2652  12540 

4)  6238.9  6230.2  2369.6  2359.9  8815.5  9254.4  1.2688  1  2550 

5)  6149.7  6143.8  2538.7  25323  92433  9559.9  1.2746  12649 

6)  5806.6  5799.9  2290.2  2281.0  94023  9670.6  1.2436  1  2359 

7)  3068.3  3084.2  1914.9  1919.6  79362  8013.3  1.2413  12396 

8)  5802.6  5799.0  2606.4  2597.8  10701.9  10961.8  1.2435  1.2370 

Max  abs  pet  diff  in  TEMPERATURE  =  0.5%  for  case  no.  7 

Max  abs  pet  diff  in  B ALLERGY  =  5.0%  for  case  no.  4 
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8.5  kJ/g  added 

T(Std)  T(Ext)  ImpCStd)  Imp(Ext)  BaU(Std)  Ball(Ext)  Gam(Std)  Gam  (Ext) 

1)  61972  6193.1  2570.2  25632  9454.0  9844.4  1.2719  12604 

2)  6242.7  6236.1  2599.9  2592.5  9452.5  9826.1  1.2750  12638 

3)  6297.9  6290.4  2524.8  25152  93242  97772  1.2708  1  2573 

4)  6383.1  6372.7  2456.2  2445.1  8923.3  9454.6  1.2753  12586 

5)  6311.7  6304.0  2631.6  2624.0  9398.8  9809.5  1.2800  12675 

6)  5926.6  5920.0  2367.8  2358.0  9558.4  9861.4  1.2477  12391 

7)  3213.3  3229.8  2022.8  2027.8  8336.6  8414.4  1.2426  12410 

8)  5913.5  59102  2684.8  2675.9  10882.4  11168.8  12467  12396 

Max  abs  pet  diff  In  TEMPERATURE  =  0.59  for  case  no.  7 
Max  abs  pet  diff  in  B ALLERGY  =  6.0%  for  case  no.  4 

9  kJ/g  added 

T(Std)  T(Ext)  ImpCStd)  ImpCExt)  BaU(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1)  6340.0  6334.6  2658.3  2650.0  95882  1006S.1  12772  12633 

2)  6397.8  6389.3  2692.0  2683.1  9591.8  10064.5  12807  12666 

3)  6439.9  6430.8  2612.9  2601.8  9438.0  99852  1.2768  12606 

4)  6528.8  6516.2  2545.0  2532.0  9012.6  9655.3  1.2824  12622 

5)  6472.9  6462.5  2725.7  27162  9527.3  10054.9  12861  12701 

6)  6046.8  6040.1  2446.8  2436.4  9708.1  10051.0  12520  12424 

7)  3360.8  3377.3  2131.6  2136.8  87482  8824.8  1.2437  12421 

8)  6023.8  6020.9  2764.2  2755.0  11056.8  11372.8  1.2500  1.2422 

Max  abs  pet  diff  in  TEMPERATURE  =  0.5%  for  case  no.  7 
Max  abs  pet  diff  in  B ALLERGY  =  7.1%  for  case  no.  4 

9.5  kJ/g  added 

T(Std)  T(Ext)  ImpCStd)  Imp(Ext)  Ball(Std)  Ball(Ext)  Gam(Std)  Gam  (Ext) 

1)  6483.0  6475.7  2747.8  2737.9  9702.7  10283.9  12832  12662 

2)  6552.6  6541.5  2785.6  2774.6  9705.6  10299.6  12870  12694 

3)  6582.8  6571.5  2702.9  2689.9  9532.4  101923  12835  12639 

4)  6676.3  6660.7  2635.9  26205  9079.9  9856.3  1.2903  12659 

5)  6633.7  6619.7  2821.1  2809.0  9625.1  10296.9  1  2931  12728 

6)  6167.4  6160.4  2527.3  2516.3  9850.7  10239.3  1  2566  1  2457 

7)  3510.0  3526.1  2240.8  2246.1  9169.4  9243.4  1  2444  1  2430 

8)  6133.7  6131.0  2844.7  2835.1  11224.7  11573.8  1.2534  1.2450 

Max  abs  pa  diff  in  TEMPERATURE  =  0.5%  for  case  no.  7 
Max  abs  pa  diff  in  B ALLERGY  =  8.6%  for  case  no.  4 
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10  kJ/g  added 


T(Std) 

T(Ext)  Imp(Std)  Imp(Ext) 

Ball(Std)  Ball(Ext)  Gam(Std)  Gam(Ext) 

1) 

6626.4 

6616.5 

2838.8 

2826.9 

9794.6 

10500.9 

1.2898 

12692 

2) 

6707.6 

6692.8 

2880.6 

2867.0 

9790.1 

10531.8 

1.2942 

12722 

3) 

6727.0 

6712.7 

2794.8 

2779.4 

9604.0 

10398.5 

12910 

12673 

4) 

6825.8 

6806.2 

2729.0 

2710.7 

9120.9 

10057.8 

12992 

1.2695 

5) 

6794.6 

6775.8 

2917.9 

2902.5 

9687.6 

10535.7 

1.3012 

12755 

6) 

6288.3 

6281.0 

2609.4 

2597.5 

9985.4 

10426.3 

1.2613 

12491 

7) 

3659.8 

3675.3 

2349.8 

2355.0 

9597.0 

9667.9 

1.2448 

12436 

8) 

6243.4 

6240.8 

2926.3 

2916.3 

11385.9 

11772.0 

1.2570 

1.2477 

Max  abs  pet  diff  in  TEMPERATURE  =  0.4%  for  case  no.  7 
Max  abs  pet  diff  in  B ALLERGY  =  10.3%  for  case  no.  4 


GRAND  TOTALS 

MEAN  ABSOLUTE  DEVIATIONS 
Temp  (K)  Impetus  (J/g)  BaUergy  (J/g) 

7.5  6.0  176.9 

MEAN  ABSOLUTE  PERCENT  DEVIATIONS 
Temp  (K)  Impetus  (J/g)  Ballergy  (J/g) 

0.174  0.31  2.01 

ROOT  MEAN  SQUARE  DEVIATIONS 
Temp  (K)  Impetus  (J/g)  Balleigy  (J/g) 

8.3  6.7  262.0 

ROOT  MEAN  SQUARE  PERCENT  DEVIATIONS 
Temp  (K)  Impetus  (J/g)  Balleigy  (J/g) 

0.20  0.32  2.8 
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Appendix  G. 


Output  from  COMPARE.BAS:  The  Effects  of 
Excluding/Including  Ions  with  McBride’s  Extended  Library 

The  two  EXMUCET  output  files  are  TSTLBONO.EXT  and  TSTLBOIO.EXT 

The  input  (SUMMARY)  files  are  TSTLBSNO.EXT  and  TSTLBSIO.EXT 

Thermo  library  used  for  first  set  was 

***  B.  McBride’s  EXTENDED-RANGE  Library  *** 

Thermo  library  used  for  second  set  was 

***  B.  McBride’s  EXTENDED-RANGE  Library  *•* 

Number  of  energies  =  21 

0  kJ/g  added 

T(no)  T(ions)  Imp(no)  Imp(ion)  Ball(no)  Ball(ion)  Gam(no)  Gam(ion) 

1)  30102  3010.2  1077.9  1077.9  45832  45832  1.2352  12352 

2)  2549.1  2549.1  943.5  9435  37365  37365  12525  12525 

3)  33922  3392.2  1136.3  1136.3  5184.0  5184.0  1.2192  12192 

4)  3811.9  3811.9  11735  11735  5723.1  5723.1  1.2050  12050 

5)  2490.6  2490.6  927.9  927.9  3588.6  3588.6  12586  12586 

6)  3738.6  3738.6  11982  11982  6278.9  6278.9  1.1908  1.1908 

7)  12165  1216.5  472.6  472.6  2899.8  2899.8  1.1630  1.1630 

8)  3574.6  3574.6  13825  13825  6590.0  6590.1  1.2098  12098 

Max  abs  pet  diff  in  TEMPERATURE  =  0.00%  for  case  no.  1 
Max  abs  pa  diff  in  B ALLERGY  =  0.00%  for  case  no.  8 

2  kJ/g  added 

T(no)  T(ions)  Imp(no)  Imp(ion)  Ball(no)  Ball(ion)  Gam(no)  Gam(ion) 

1)  4081.4  4081.4  1482.8  1482.8  64105  6410.6  1.2313  12313 

2)  37745  3774.5  14065  14065  5794.3  5794.4  1.2427  12427 

3)  4329.4  4329.4  1487.0  1487.0  6738.9  6738.9  1.2207  12207 

4)  4521.1  4521.1  1457.0  1457.0  68022  6802.3  1.2142  12142 

5)  3745.4  3745.4  1404.3  1404.3  5661.6  5661.6  1.2480  12480 

6)  43155  43155  1451.5  14515  7255.3  7255.3  1.2001  12001 

7)  1562.4  1562.4  735.3  735.3  3886.6  3886.6  1.1892  1.1892 

8)  4301.6  4301.6  1706.7  1706.7  8062.8  8062.8  1.2117  12117 

Max  abs  pa  diff  in  TEMPERATURE  =  0.00%  for  case  no.  1 
Max  abs  pa  diff  in  B ALLERGY  -  0.00%  for  case  no.  2 
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4  kj/g  added 

T(no)  T(ions)  Imp(no)  Imp(ion)  Ball(no)  Ball(ion)  Gam(no)  Gam(ion) 

1)  4844.6  4844.6  1820.6  1820.6  7681.0  7680.9  1.2370  1.2370 

2)  4709.6  4709.6  1796.3  1796.3  7360.3  7360.3  1.2441  1.2441 

3)  5002.4  5002.4  1793.8  1793.8  7810.4  7810.4  1.2297  1.2297 

4)  5106.4  5106.4  1737.1  1737.1  7653.2  7653.1  1.2270  1.2270 

5)  4715.9  4715.9  1807.2  1807.2  7266.4  7266.4  1.2487  1.2487 

6)  4828.6  4828.6  1710.9  1710.9  80992  80992  12112  12112 

7)  1974.5  1974.5  1081.3  1081.3  4984.6  4984.6  12169  12169 

8)  4857.9  4857.8  2000.6  2000.6  9159.6  9159.6  12184  12184 

Max  abs  pet  diff  in  TEMPERATURE  =  0.00%  for  case  no.  8 
Max  abs  pet  diff  in  B ALLERGY  &  0.00%  for  case  no.  4 

6  kJ/g  added 

T(no)  T(ions)  Imp(no)  Imp(ion)  Ball(no)  Ball(ion)  Gam(no)  Gam(ion) 

1)  5471.3  5471.2  2144.3  21442  8698.6  8698.5  1.2465  12465 

2)  5439.3  5439.2  21502  21502  8557.6  8557.5  1.2513  12513 

3)  5588.0  5587.9  2103.1  2103.0  8716.4  8716.3  12413  12413 

4)  56672  5667.1  2036.1  2036.0  8456.3  84562  1.2408  12408 

5)  54762  54762  2171.7  2171.7  8502.0  8501.9  1.2554  12554 

6)  5318.5  5318.5  19862  19862  8896.4  8896.3  1.2233  12233 

7)  25152  2515.2  1495.9  1495.9  6456.5  6456.5  12317  12317 

8)  5343.9  5343.8  2294.1  2294.1  10099.8  10099.6  1.2271  1.2271 

Max  abs  pet  diff  in  TEMPERATURE  =  0.00%  for  case  no.  8 
Max  abs  pet  diff  in  B ALLERGY  =  0.00%  for  case  no.  8 

8  kJ/g  added 

T(no)  T(ions)  Imp(no)  Imp(ion)  Ball(no)  Ball(ion)  Gam(no)  Gam  (ion) 

1)  6051.3  6051.1  2477.5  2477.4  9621.7  9621.4  1.2575  12575 

2)  6081.8  6081.6  2502.8  2502.7  9584.0  9583.7  1  2611  12611 

3)  6150.6  6150.3  2430.1  2430.0  9568.6  95682  1.2540  12540 

4)  6230.5  6230.2  2360.1  2359.9  9254.8  9254.4  1.2550  12550 

5)  6144.0  6143.8  2532.6  2532.5  95602  9559.9  1.2649  1  2649 

6)  5800.1  5799.9  2281.1  2281.0  9671.0  9670.6  1.2359  12359 

7)  30842  3084.2  1919.6  1919.6  8013.3  8013.3  1  2396  1  2396 

8)  5799.1  5799.0  2597.9  2597.8  10962.1  10961.8  1.2370  1.2370 

Max  abs  pet  diff  in  TEMPERATURE  =  0.00%  for  case  no.  3 
Max  abs  pet  diff  in  B ALLERGY  =  0.00%  for  case  no.  4 
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10  kJ/g  added 

T(no)  T(ions)  Imp(no)  Imp(k»)  Ball(oo)  BalKkm)  Gam(no)  Gam(ion) 

1)  6617.0  66163  2827.2  2826.9  10S01.7  10S00.9  1.2692  1.2692 

2)  6693.2  6692.8  2867.3  2867.0  10532.5  10531.8  1.2722  1.2722 

3)  6713.3  6712.7  2779.7  2779.4  10399.3  10398.5  1.2673  1.2673 

4)  6807.0  6806.2  2711.1  2710.7  10058.8  10057.8  1.2695  1.2695 

5)  6776.3  6775.8  2902.8  2902.5  10536.4  10535.7  1.2755  1.2755 

6)  6281.4  6281.0  2597.7  2597.5  10427.0  10426.3  1.2491  12491 

7)  3675.3  3675.3  2355.0  2355.0  9667.9  9667.9  1.2436  1.2436 

8)  6241.1  6240.8  2916.5  2916.3  11772.6  11772.0  12477  1.2477 

Max  abs  pet  diff  in  TEMPERATURE  =  0.01%  for  case  no.  4 
Max  abs  pa  diff  in  B  ALLERGY  *  0.01%  for  case  no.  4 

12  kJ/g  added 

T(no)  T(ions)  Imp(no)  Imp(ion)  Ball(no)  Ball(ion)  Gam(no)  Gam(ion) 

1)  71793  7178.4  3193.1  31923  11354.1  11352.6  1.2812  12812 

2)  7292.8  7291.9  3244.5  3244.0  11437.1  11435.8  1.2837  12837 

3)  7284.7  7283.6  3151.8  3151.1  11216.7  11215.1  1.2810  1.2810 

4)  7400.1  7398.6  3087.9  3087.0  10869.3  10867.3  12841  12841 

5)  73903  7389.3  3281.4  3280.8  11466.6  11465.0  12862  1.2862 

6)  6768.4  6767.6  2937.6  2937.1  11159.7  11158.6  1.2632  12632 

7)  42572  42572  2780.0  2780.0  11363.7  11363.7  12446  12446 

8)  6679.3  6678.7  32522  3251.7  12536.1  12535.1  12594  12594 

Max  abs  pa  diff  in  TEMPERATURE  =  0.02%  for  case  no.  4 
Max  abs  pa  diff  in  BALLERGY  =  0.02%  for  case  no.  4 

14  kJ/g  added 

T(no)  T(ions)  Imp(no)  ImpQon)  Ball(no)  Ball(ion)  Gam(no)  GamQon) 

1)  7738.3  7736.3  3570.7  35693  12181.8  12179.1  12931  12931 

2)  7879.6  7877.6  3627.9  3626.7  12308.3  12305.7  12948  12947 

3)  7863.8  7861.6  3542.1  3540.7  12020.0  120172  12947  12946 

4)  8006.1  8003.2  34852  34833  11685.0  11681.4  1.2983  1.2982 

5)  7981.9  7979.8  3660.0  3658.7  123623  12359.6  1.2961  12960 

6)  7265.4  7264.1  33012  33002  118592  11857.3  12784  12783 

7)  4792.3  4792.3  3182.0  3182.0  12958.8  12958.8  12455  12455 

8)  7119.1  7118.1  3605.8  3605.1  13246.9  13245.4  12722  1.2722 

Max  abs  pa  diff  in  TEMPERATURE^  0.04%  for  case  no.  4 
Max  abs  pa  diff  in  BALLERGY  »  0.03%  for  case  no.  4 
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16  kJ/g  added 

T(no)  T(ioos)  lmp(no)  Imp(ion)  Ball(no)  Ball(ion)  Gam(no)  Gam(ion) 

1)  8285.8  8282.5  3952.6  3950.5  12984.1  12979.8  1.3044  1.3044 

2)  84425  8439.0  40075  4005.3  13149.9  13145.4  1.3048  1.3047 

3)  8441.0  8437.0  3942.0  39395  12806.4  12801.9  1.3078  1.3077 

4)  86122  86072  3893.6  38905  125042  12498.3  1.3114  1.3113 

5)  8540.0  8536.2  4029.3  4027.0  13231.6  132265  1.3045  1.3045 

6)  7774.9  77725  3688.1  3686.4  12511.0  12508.4  12948  1.2947 

7)  52785  52785  3566.1  3566.1  14408.9  14408.9  12475  12475 

8)  7563.7  7562.0  3977.4  3976.1  13892.6  138905  12863  1.2862 

Max  abs  pet  diff  in  TEMPERATURE  =  0.06%  for  case  no.  4 
Max  abs  pet  diff  in  B ALLERGY  =  0.05%  for  case  no.  4 

18  kJ/g  added 

T(no)  T(kms)  Imp(no)  Imp(ion)  Ball(no)  Ball(ion)  Gam(no)  Gam  (ion) 

1)  8812.9  88075  43315  4328.1  137622  13755.8  1.3147  1.3146 

2)  8972.0  89662  4376.6  4373.1  13970.1  13963.0  1.3133  1.3132 

3)  9001.8  8995.3  43415  4337.3  13576.0  13569.1  1.3198  1.3196 

4)  9201.1  9192.9  43015  4296.3  133272  13318.1  1.3228  1.3226 

5)  9059.1  9053.0  4385.7  4382.0  14087.0  14078.6  1.3113  1.3113 

6)  8296.3  8292.3  4096.4  40935  13098.1  13094.8  1.3127  1.3126 

7)  5726.1  5726.1  3940.2  39402  15716.9  15716.8  12507  1.2507 

8)  8014.0  8011.3  4366.0  4363.8  144562  14453.7  1.3020  1.3019 

Max  abs  pet  diff  in  TEMPERATURE  =  0.09%  for  case  no.  4 
Max  abs  pet  diff  in  B  ALLERG  Y  =  0.07%  for  case  no.  4 

20  kJ/g  added 

T(no)  T(ions)  Imp(no)  Imp(ion)  Ball(no)  BallQon)  Gam(no)  Gam  (ion) 

1)  9313.7  9305.6  4703.6  46985  145225  14513.4  1.3239  15237 

2)  9465.9  94572  47345  4728.8  14785.6  14774.6  1.3202  15201 

3)  9534.9  9525.0  4733.2  4726.7  14337.1  143275  1.3301  1.3299 

4)  9759.1  9746.7  4700.6  4692.7  14162.7  14149.0  1.3319  1.3317 

5)  9541.3  9532.2  4730.8  47252  14948.9  14935.8  1.3165  15164 

6)  88265  8819.9  4522.3  45175  136025  13599.4  1.3325  1.3322 

7)  6145.3  61452  4310.3  43102  16892.1  16892.0  12552  1.2552 

8)  8469.3  8465.1  4769.7  4766.3  149185  14916.3  1.3197  1.3195 

Max  abs  pet  diff  in  TEMPERATURE  *  0.1%  for  case  no.  4 
Max  abs  pet  diff  in  B  ALLERGY  =  0.10%  for  case  no.  4 
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22  kJ/g  added 

T(no)  T(ions)  Imp(no)  ImpCion)  Ball(no)  Ball(ion)  Gtm(no)  Gam(ion) 

1)  9787.7  9776.1  5068.6  5060.8  15277.8  15265.9  1.3318  1.3315 

2)  9928.0  9915.5  5083.3  50752  156202  15603.9  1.3254  1.3253 

3) 10036.4  10021.9  51152  5105.5  151103  15096.7  1.3385  1.3382 

4)  10281.5  10263.8  5088.6  5077.0  15034.1  15013.8  1.3385  1.3382 

5)  9992.5  9979.6  50685  5060.3  15843.4  158235  1.3199  1.3198 

6)  9358.8  9348.9  4960.9  4953.4  140105  14009.4  1.3541  1.3536 

7)  6543.9  6543.7  4680.1  4680.0  17939.4  17939.0  12609  12609 

8)  8927.6  8921.3  51862  5180.9  152605  15260.3  1.3398  1.3395 

Max  abs  pet  diff  In  TEMPERATURE  =  02%  for  case  no.  4 
Max  abs  pet  diffin  BALLERGY  =  0.1%  for  case  no.  4 

24  kJ/g  added 

T(no)  T(ions)  Imp(no)  ImpCion)  Ball(no)  Ball(ion)  Gam(no)  Gam(ion) 

1)  10238.1  10222.1  5428.9  5417.8  16050.8  16035.6  1.3382  1.3379 

2)  10364 2  10347.0  54272  5415.7  16505.6  164820  15288  1.3286 

3)  10508.4  10488.3  54892  5475.4  159312  159120  1.3446  1.3441 

4) 10770.9  10746.7  5466.9  5450.6  15984.8  15954.9  1.3420  1.3416 

5)  10419.1  10401.6  5402.4  5391.0  168032  16773.9  1.3215  1.3214 

6)  98872  9872.3  54072  5395.7  14319.8  14324.6  1.3776  1.3767 

7)  6927.3  6927.0  5052.1  5051.8  18857.4  18856.7  12679  1.2679 

8)  93862  9376.9  56132  56052  154672  15472.1  1.3629  1.3623 

Max  abs  pet  diff  in  TEMPERATURE  =  02%  for  case  no.  4 
Max  abs  pet  diff  in  BALLERGY  =  02%  for  case  no.  4 

26  kJ/g  added 

T(no)  T(ions)  Imp(no)  ImpCion)  BalKno)  BaU(ion)  Gam(no)  Gam  (ion) 

1)  10669.7  106483  5787.1  5771.9  16876.9  16858.0  1.3429  1.3424 

2) 10779.7  10756.9  5769.0  5753.4  17484.4  17451.0  1.3300  13297 

3)  10955.6  10928.8  5857.7  5838.8  16852.6  16825.9  1.3476  1.3470 

4) 112323  112003  5838.0  5815.9  17083.8  170393  13417  1.3413 

5)  10826.4  10803.3  5735.4  5719.8  17871.7  178293  1.3209  1.3208 

6)  10406.0  10384.3  5857.6  5840.4  145502  145672  1.4026  1.4009 

7)  7299.3  7298.7  5427.4  5426.8  19639.7  19638.7  12763  12763 

8)  9843.1  9829.6  6049.0  60373  15532.3  15547.9  1.3894  1.3883 

Max  abs  pet  diff  in  TEMPERATURE  =  0.3%  for  case  no.  4 
Max  abs  pet  diff  in  BALLERGY  =  0.3%  for  case  no.  4 
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28  kJ/g  added 

T(no)  T(ioos)  Imp(no)  Imp(ion)  Ball(no)  Ball(ion)  Gam(no)  Gam(ion) 

1)  11086.8  110S8.6  6145 3  6124.9  17814.7  17790.6  1.3450  1.3443 

2)  11178.7  11149.1  6110.8  6089.9  18616.0  18569.7  1.3283  1.3279 

3)  11382.5  11347.5  6222.8  6197.4  17954.3  17917.2  1.3466  1.3459 

4)  11670.7  11629.5  6202.7  6173.7  18435.7  18369.9  1.3365  1.3361 

5)  11218.1  11188.3  6068.8  6048.1  19109.2  19049.7  1.3176  1.3175 

6)  10912.2  10881.7  6309.6  6284.8  14757.7  14796.0  1.4275  1.4248 

7)  7662.4  7661.4  5807.0  5806.0  20276.6  20275.3  1.2864  1.2864 

8)  10296.9  10277.8  6493.1  6476.1  15465.9  15501.6  1.4198  1.4178 

Max  ate  pet  diff  in  TEMPERATURE  *  0.4%  for  case  no.  4 
Max  ate  pet  diff  in  B ALLERGY  *=  0.4%  for  case  no.  4 

30  kJ/g  added 

T(no)  T(ions)  Imp(no)  Imp(ion)  Ball(no)  Ball(ion)  Gam(no)  Gam  (ion) 

1)  11492.5  11455.9  6504.6  6477.3  18964.9  18930.9  1.3430  1.3422 

2)  11563.6  11526.0  6452.8  6425.3  19985.2  19921.7  1.3229  1.3225 

3)  11791.3  11746.7  6584.0  6550.9  19355.2  193014  1.3402  1.3394 

4)  12087.0  12035.2  6559.1  65210  201953  200962  1.3248  1.3245 

5)  11596.3  11558.6  6402.7  6375.7  20600.2  20517.8  1.3108  1.3107 

6)  11404.3  11362.4  67612  6726.3  15059.8  15129.8  1.4490  1.4446 

7)  8018.5  8016.8  6191.4  6189.8  20756.0  20754.7  1.2983  1.2982 

8)  10747.5  10720.8  6945.3  6920.9  15309.8  15378.8  1.4537  1.4500 

Max  ate  pet  diff  in  TEMPERATURE  =  0.4%  for  case  no.  4 
Max  ate  pet  diff  in  BALLERGY  *  0.5%  for  case  no.  4 

32  kf/g  added 

T(no)  T(ions)  Imp(no)  ImpOon)  Ball(no)  Ball(ion)  Gam(no)  Gam  (ion) 

1)  11886.9  11840.3  6861.7  6826.4  205052  20444.7  1.3346  1.3339 

2)  11934.9  11887.8  67932  6757.9  21712.3  21625.5  1.3129  1.3125 

3)  12181.6  12125.9  6937.9  6895.6  21233.4  21152.4  1.3267  1.3260 

4)  12478.8  124153  6900.7  6854.9  22588.7  22430.4  1.3055  1.3056 

5)  11961.3  11914.4  6735.3  6700.8  22463.1  22349.6  12998  12998 

6)  11879.0  11822.7  7206.7  7159.1  15687.3  15791.8  1.4594  1.4533 

7)  8368.9  8366.3  6581.3  6578.7  210643  21064.0  1.3124  1.3123 

8)  11194.6  11157.7  7404.7  7370.3  15171.7  152893  1.4881  1.4820 

Max  ate  pa  diff  in  TEMPERATURE  =  0.5%  for  case  no.  4 
Max  ate  pa  diff  in  BALLERGY  =  0.8%  for  case  no.  8 
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34  kJ/g  added 

T(no)  T(ions)  Imp(no)  lmp(kn)  Ball(no)  BaU(k»)  Gam(no)  GamQon) 

1)  12265.7  12208.2  7208.6  7165.0  22752.7  22614.3  1.3168  1.3168 

2)  12290.8  12232.8  7127.8  7083.4  23969.7  23844.9  12974  12971 

3)  12548.7  12481.0  7276.0  72242  238542  23708.9  1.3050  1.3047 

4)  12838.6  12764.0  7215.5  7162.4  259362  25655.9  12782  12792 

5)  123112  12253.9  7062.1  7019.0  24863.9  24703.1  12840  12841 

6)  12326.3  122542  76312  7570.3  17084.0  171805  1.4467  1.4406 

7)  8714.8  8710.8  6977.3  6973.4  21187.9  21190.0  1.3293  1.3291 

8)  11635.2  11585.4  7866.0  7818.7  15307.3  154752  12139  12052 

Max  abs  pet  diff  in  TEMPERATURE  *  0.6%  for  case  no.  6 
Max  abs  pa  diff  in  B  ALLERGY  *  1.1%  for  case  no.  8 

36  kJ/g  added 

T(no)  T(ions)  Imp(no)  Imp(ion)  Ball(no)  BaIl(ion)  Gam(no)  Gam(ion) 

1)  12617.3  12550.1  7528.6  7479.1  26237.6  258862  12869  12889 

2)  12626.3  12556.9  7448.3  7394.4  27002.8  26799.9  12758  12759 

3)  12883.3  12804.8  7584.3  7525.1  27599.0  27291.7  12748  12757 

4)  13154.4  13072.1  7486.7  7431.0  30658.7  30111.5  12442  12468 

5)  12641.4  12573.0  7375.4  7323.4  28034.7  277862  12631  12636 

6)  12724.9  12640.3  8006.7  7936.9  20011.7  19921.3  1.4001  1.3984 

7)  90572  9051.4  7380.7  7374.8  21111-5  21119.6  1.3496  1.3492 

8)  12058.3  11994.1  8312.3  8251.5  16285.1  16430.7  15104  15022 

Max  abs  pa  diff  in  TEMPERATURE  ■  0.7%  for  case  no.  6 
Max  abs  pa  diff  in  BALLERGY  *  1.8%  fa  case  no.  4 

38  kJ/g  added 

T(no)  T(ions)  Imp(no)  Imp(ion)  BaU(no)  Ball(ion)  Gam(no)  GamOon) 

1)  12923.0  12851.6  77972  7749.0  316375  308255  12465  12514 

2) 12933.1  12853.4  7742.1  7680.7  31149.4  30766.7  12485  12496 

3) 13172.7  13088.1  78445  7784.3  32926.6  32258.1  12382  12413 

4) 13415.3  13332.0  7699.0  7648.9  37147.8  36105.1  12073  12119 

5)  12944.0  12865.6  7663.4  76042  32290.6  31858.8  12373  12387 

6)  13049.0  129615  8298.0  8232.1  25333.3  24710.3  15276  15331 

7)  93975  9389.0  7792.8  77842  20821.0  20840.6  1.3743  1.3735 

8) 124385  12363.9  8707.8  8639.7  19125.9  189942  1.4553  1.4549 

Max  abs  pa  diff  in  TEMPERATURE  *  0.7%  for  case  no.  6 
Max  abs  pa  diff  in  BALLERGY  =  2.8%  fa  case  no.  4 
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T(no) 

40  kJ/g  added 
T(ions)  Imp(no)  Imp(ion) 

Ball(no)  Ball(ion)  Gam(no)  Gam  (ion) 

1) 

13169.0 

13100.7 

7996.9 

7958.8 

39349.4 

378326 

12032 

12104 

2) 

13200.9 

13114.9 

7993.5 

7930.4 

36797.3 

36040.3 

12172 

12200 

3) 

13408.3 

13324.6 

8044.0 

7991.5 

40147.0 

38868.1 

12004 

12056 

4) 

13619.9 

13542.0 

7849.6 

7812.3 

45484.4 

43748.9 

1.1726 

1.1786 

5) 

13209.4 

13124.8 

79112 

7850.3 

37994.7 

37196.7 

12082 

12110 

6) 

132912 

13210.8 

8495.5 

8444.9 

33274.1 

31838.4 

12553 

12652 

7) 

9736.1 

9724.6 

82152 

8202.8 

20302.5 

20342.8 

1.4046 

1.4032 

8) 

12746.4 

126712 

9011.9 

8949.9 

24833.7 

24038.9 

1.3629 

1.3723 

Max  abspct 

diff  in  TEMPERATURE  = 

0.7%  for  case  no. 

2 

Max  abs  pet 

diff  in  BALLERGY  «  4.3%  for  case  no.  6 

GRAND  TOTALS 


MEAN  ABSOLUTE  DEVIATIONS 
Temp  (K)  Impetus  (J/g)  B  allergy  (J/g) 

20.7  15.3  1052 

MEAN  ABSOLUTE  PERCENT  DEVIATIONS 
Temp  (K)  Impetus  (J/g)  B  allergy  (J/g) 

0.174  0.22  0.36 

ROOT  MEAN  SQUARE  DEVIATIONS 
Temp  (K)  Impetus  (J/g)  Ballergy  (J/g) 

34.6  25.7  297.4 

ROOT  MEAN  SQUARE  PERCENT  DEVIATIONS 
Temp  (K)  Impetus  (J/g)  Ballergy  (J/g) 

0.28  0.35  0.9 
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Appendix  H. 

Comparison  of  the  Species 
with  Mole  Fractions  Greater  than  1.0  x  10' 
Produced  by  Two  Different  Thermo  Libraries 


NOTE:  These  products  are  for  Composition  No.  1  with  5000  J!{  of  added  Energy. 


OLD  (STD)  LIB 

NEW  (EXTENDED)  LIB 

C2H2 

3.936E-07 

C2H2 

3.391E-07 

C 

2.005E-06 

C 

2.031E-06 

C2HRAD 

2.412E-07 

C2H 

1.251E-07 

COO  RAD 

2.128E-06 

CCORAD 

1.685E-06 

CH 

2.979E-06 

CH 

2.642E-06 

CH2 

4.722E-06 

CH2 

4.812E-06 

CH3 

5.728E-06 

CH3 

5.250E-06 

CH4 

1.553E-06 

CH4 

1.830E-06 

CN 

4.147E-0S 

CN 

3.729E-Q5 

CO 

2.813E-01 

CO 

2.823E-01 

C02 

3.857E-02 

C02 

3.774E-02 

CH20 

8.209E-05 

CH20 

6.732E-05 

HCOOH 

8.562E-05 

HCOOH 

8.320E-05 

H 

4.489E-02 

H 

4.477E-02 

H2 

1.335E-01 

H2 

1.316E-01 

H2N2 

1.U3E-06 

— 

.... 

H20 

1.680E-01 

H20 

1.708E-01 

H202 

7.128E-OS 

H202 

8.200E-05 

HCN 

3.197E-04 

HCN 

3.241E-04 

HCORAD 

1.053E-03 

HCO 

1.035E-03 

HNCO 

8.314E-05 

HNCO 

8.182E-05 

HNO 

2.890E-O4 

HNO 

2.999E-04 

HN02 

2349E-05 

HN02 

2.410E-05 

H02 

2.487E-04 

H02 

2.332E-04 

CH30 

9.268E-07 

— 

-  -  -  - 

KETENE 

3.037E-07 

KETENE 

2.735E-07 

i 

a 

OLD  (STD)  UB 

NEW  (EXTENDED)  LIB 

H 

H 

CH30H 

1369E-07 

CH30H 

1.509E-07 

H 

N 

3.874E-04 

N 

3.873E-04 

1 

N2 

2.488E-01 

N2 

2.494E-01 

H 

N2H2 

1.847E-06 

N2H2 

1.789E-06 

n 

N20 

4.212E-05 

N20 

4.106E-05 

u 

N3 

2.328E-07 

N3 

2.303E-07 

H 

N3H 

2.425E-07 

N3H 

2.372E-07 

1 

NCNRAD 

2.704E-07 

NCN 

1.931E-05 

1 

NCO 

2.678E-05 

NCO 

1.752E-05 

1 

NH 

8.050E-04 

NH 

5215E-04 

1 

NH2 

4.908E-04 

NH2 

4.789E-04 

n 

NH20H 

3.052E-06 

NH20H 

2.88QE-06 

R 

NH3 

3.109E-04 

NH3 

4.012E-04 

B 

NO 

1.746E-02 

NO 

1.712E-Q2 

It 

N02 

4.068E-05 

NQ2 

3.892E-05 

R 

0 

8.990E-03 

0 

8.842E-03 

I 

02 

4.309E-03 

02 

4.179E-03 

R 

03 

1.912E-07 

03 

6346E-07 

H 

OH 

4.982E-02 

OH 

4.877E-Q2 

R 

— 

-  -  -  - 

CH20H 

1J84E-06 

R 

— 

.... 

CNN 

1.957E-06 

R 

— 

-  -  -  - 

CHCO 

2.409E-07 

R 

— 

.  -  -  - 

COOH 

2.663E-04 

R 

... 

.... 

HNC 

1.368E-05 

R 

... 

.... 

CH30 

1.734E-07 

H 
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Intentionally  left  blank. 
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Appendix  1. 

Listing  of  Thermodynamic  Tables  Supplied  by  NIST 


ARL  Requested  NIST  Sent  Tables  ARL  Requested  NIST  Sent  Tables 

COEFS  for  for  COEFS  for  for 


c 

Carbon,  monatomic 

c+ 

Carbon,  unipositive 

c- 

Carbon,  uninegative 

C2 

CCH  (g)  Radical 

C2- 

CH 

CH  (g) 

CH+ 

CH+  (g) 

CH2 

CH2  (g) 

CN 

CN  (g) 

CN+ 

CN-f  (g) 

CN- 

CN-  (g) 

CO 

C02- 

C02-  (g) 

— 

electron  gas 

— 

Hydrogen,  monatomic 

H+ 

Hydrogen,  unipositive 

H- 

Hydrogen,  uninegative 

H2+ 

H2- 

N 

Nitrogen,  monatomic 

N- 

Nitrogen,  uninegative 

N2 

N+ 

Nitrogen,  unipositive 

N2+ 

N2- 

N20f 

NO 

NO  (g) 

NOf  (g) 

N02  (g) 

N02-  (g) 

0 

Oxygen,  monatomic 

Of 

Oxygen,  unipositive 

0- 

Oxygen,  uninegative 

02 

02- 

OH 

OH  (g) 

0H+ 

0H+  Hydroxyl  positive 

OH- 

OH-  Hydroxide 

ONO- 

Ti 

Titanium,  monatomic 

Ti+ 

Titanium,  unipositive 

TiO 

TiO  (g) 

— 

Titanium,  uninegative 

Intentionally  left  blank. 
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Appendix  J. 

Comparison  of  Thermodynamic  Tables 
Produced  by  Back-Calculation  from  the  LRC  Coefficients 
with  Those  Furnished  by  NIST 

Units  are  Joules,  moles,  and  Kelvins 
NAME/ 

HEAT  CAPACITY  ENTROPY  ENTHALPY  Number 


MnAvDev  Pet  RMS  MnAvDev 

Pet  RMS 

MnAvDev 

Pet  RMS 

of  Pts 

c+ 

Carbon,  unipositive  ion 

0.002  0.020  0.000 

0.001 

1.2 

0.00 

102 

electron  gas 

0.004  0.019  -0.000 

0.005 

0.0 

0.01 

102 

H+ 

Hydrogen,  unipositive  ion 

0.004  0.019  -0.000 

0.002 

0.0 

0.01 

102 

H- 

Hydrogen,  uninegative  ion 

0.004  0.019  -0.000 

0.002 

0.0 

0.01 

102 

O 

Oxygen,  monatomic 

-0.001  0.073  -0.002 

0.002 

-8.3 

0.01 

102 

o+ 

Oxygen,  unipositive  ion 

0.001  0.029  0.001 

0.002 

1.1 

0.01 

102 

N+ 

Nitrogen,  unipositive  ion 

0.000  0.016  0.001 

0.001 

0.3 

0.01 

102 

N- 

Nitrogen,  uninegative  ion 

0.001  0.013  0.001 

0.001 

0.2 

0.01 

102 

N 

Nitrogen,  monatomic 

-0.003  0.097  -0.000 

0.002 

-18.1 

0.02 

102 

H 

Hydrogen,  monatomic 

-0.001  0.167  -0.003 

0.003 

-17.0 

0.03 

102 

0- 

Oxygen,  uninegative  ion 

•0.001  0.026  0.002 

0.002 

-11.1 

0.03 

102 
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C  Carbon,  monatomic 


0.007 

0.268 

-0.004 

OH-  Hydroxyl  negative  km 

0.133 

0.409 

0.0S2 

N02-  Nitrogen  dioxide  negative  km 

0.331 

1.090 

-0.656 

CN-  Cyanide 

0.133 

0.461 

1.131 

CH2 

Methylene 

-0.070 

0.907 

1.103 

NO+  Nitric  oxide  positive  ion 


0.969  5.001 

0.186 

OH  Hydroxyl 

-1.654  9.444 

0.187 

NO  Nitric  oxide 

1.668  7.135 

0.409 

OH+  Hydroxyl  positive  ion 
0.858  4.395 

0.925 

CN  Cyanogen 

1.794  5.198 

1.042 

CH  Methylidyne 
-3.148  12.080 

-0.999 

CCH  Ethynyl  Radical 
-1.749  5.545 

9.548 

N02  Nitrogen  dioxide 

4.763  10.896 

3.027 

C-  Carbon,  uninegative  ion 
-3.654  17.626 

-1.720 

CN+  Cyanogen,  ion 

1.978  8.993 

-5.291 

CH+  Methylidyne,  km 

3.852  23.368 

8.505 

0.004 

-26.3 

0.04 

102 

0.032 

407.3 

0.20 

102 

0.207 

221.0 

0.52 

62 

0.405 

795.9 

038 

102 

0.393 

419.9 

0.72 

62 

0.103 

1221.3 

0.75 

102 

0.148 

3.7 

1.07 

102 

0.213 

2798.2 

138 

102 

0.423 

461 6 J5 

2.24 

102 

0.430 

5074.3 

2.29 

102 

0.480 

-5007.8 

2.69 

102 

2.966 

-3400.1 

4.20 

102 

1.007 

7630.9 

4.43 

62 

1.188 

-5994.0 

8.22 

62 

2.514 

13159.2 

10.87 

102 

3.711 

30723.0 

16.33 

102 
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Appendix  K. 

Listing  of  Program  TK3FIT5 


PROGRAM  TIOFTTS 


This  version:  26  January  1993 

Specifically  adapted  to  creating  coefficients  for  EXBLAKE.  This 
program  is  die  filth  of  a  aeries  of  programs  named  either  TIGFm 
or  TIGFmX  0  *  1,6).  This  program  is  specifically  designed  to 
get  its  input  data  from  the  tables  given  to  ARL/WTD/APCB  by  Ms. 

Bonnie  McBride  (NASA  Lewis  Research  Center).  Her  tables  have  9 
columns  per  temperature.  For  use  with  this  program  they  have  been 
copied  into  TABLES  JN;  the  sentinel  line, 

9999.,  999.,  99„  9„  0.9, 0.09  , 99.,  99.,  99. 
has  been  added,  and  the  title  line  (which  contains  the  name  of  the 
species)  has  been  changed.  A  *$’  has  been  placed  in  ad  1  as  a 
sentinel,  and  the  name  of  die  species  staffs  in  ad  8. 

See  TIGF1T6  for  a  version  designed  for  input  from  the  JANAF  Tables 
tape.  It  would  not  be  difficult  to  modify  the  present  program  to 
accept  one  or  the  other  of  these  inputs,  but  it  is  easier  to  have  a 
separate  program. 

This  one  has  an  of  die  features  of  TIGFTT4  (use  of  the  *$’ 
sentinel:  computation  of  measures  of  goodness  of-fit);  A  N  D  it 
includes  the  cut-off  at  15  000  K. 

C  See  the  Comments  at  the  beginning  of  T1GFIT4  for  further  details. 

C  In  this  version,  die  output  files  have  been  changed  from  TIGFIT.OUT 
C  and  TIGFIT.AUX  to  TK3FTTS.OUT  and  TKHTTS AUX. 

IMPLICIT  REAL*8  (A-H.O-Z) 

C 

DIMENSION  A(7,7),  B(9),  KX(7),  JUNK(80) 

DIMENSION  TK(500),  CPP(500),  SS(500),  HH(500) 

CHARACTER  XNAME*  10PHASE*6,  JUNK*1,  SENT*1 
PARAMETER  (N6-0) 

LOGICAL  HITEMP 

COMMON  /COEFS/  B.  NUM 

COMMON  /FUNCS/  TK,  CPP,  SS,  HH,  THX,  TSTART,  ENTRP0,  HOFJ, 
$  XNAME,  PHASE 
COMMON  /IOSTUF/  U  LO 

COMMON  /LSTSQ/  T,  T2,  T3,  T4,  T5,  T6,  TM2,  TM3,  TM4,  TM5,  TM6. 
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$  Y,  YT2,  YT3,  TM,  SN,  YTM2,  YT,  YTM,  YTM3 

DATA  PHASE/ ’GAS  ’/ 

C 

LO-6 

U-5 

OPEN  (LIJFILE*  TEMPS  JN,^TATUS=,OLD’) 

OPEN  CSJTL^TIGFTrS.OUT’^TATUS-’UNKNOWN’) 

OPEN  (LO,FILE»’TIGFIT5.AUX’,STATUS» ’UNKNOWN’) 

OPEN  (4,FILE**CHEKK5’,STAT,US« ’UNKNOWN’) 

OPEN  OFILE-‘CALORIC‘ACCESS«‘APPEND‘) 

WRITE  <2,  81) 

81  FORMAT  (’  The  data  cards  for  the  TIGER  code  are  (units  =  calories’ 

$  ’.  mols,  A  K):  ’) 

C 

C  CALL  LOGO 
C 

REWIND  U 

READ  (LI.  *  END=33)  TFIX.  TSTART 
CLOSE  (LD 

OPEN  (UJFILE=TABLESJN’^TATUS=’OLD’) 

REWIND  U 

10  READ(LI,1020J5ND=33)  SENT.  XNAME 
1020  FORMAT(Al  ,6X,  A6  ) 

IF  (SENT  .NE.  ’$’)  GOTO  10 

WRITE  (N6.2008)  XNAME 
2008  FORMAT(lX,’Processing  — >  ’.A10) 

C  Initialize  die  dements  that  will  make  up  the  matrix 
CALL  FMATR  (0.  XI.  X2) 

C 

C  Skip  the  header  line  in  the  table 
READ  (LL  499.  END*33)  JUNK 
499  FORMAT  (80A1) 

C  Read  the  first  line  of  data. 

IS  READOL*.  END-33)  TEM.  CP,  XI,  HMH298,  ENTRPO,  X2,  X3,  HOFJ,  X4 
C  Is  temperature  the  fixed  temperature?  If  not,  read  another  line. 

IF  (TEM.  LT.  TFDC-0.05  .OR.  TEM  .GT.  TFIX+0.05  )  GOTO  15 
C  Yes,  it  is  die  fixed  T.  Is  T  Fix  *  T  Start? 

C  If  it  is,  then  start  saving  data  and  continue  to  read  the  table 
20  IF  (TEM  .GE  TST ART-0.05  .AND.  TEM  LE.  TSTART+0.05  )  THEN 
KT«1 

TK(1)»  TSTART 
CPP(1)*CP 
SS(1)  =  ENTRPO 
HH(1)  *  HMH298 
HTTEMP  *  .FALSE. 

GOTO  35 

END  IF 
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C  No,  T  Fix  want  T  Start.  Keep  looting  for  T  Start 
25  READ  (U  *,  END-33,  ERR-999)  TEM,  CP.  XI.  HMH298,  ENTROP.  X2. 
$  X3.  ENTHALP.  X4 

IF  (TEM  IT.  TSTART-O.Q5  .OR .  TEM  .GT.  TSTART+0.05  )  GOTO  25 
C  At  last  T  Suit  has  been  found;  save  values  and  keep  reading  table. 

29KT-1 

TK(1)  «  TSTART 
CPP(1)«CP 
SS(1)  -  ENTROP 
HH(1)  «  HMH298 
HTTEMP  *  .FALSE. 

CALL  FMATR  (1.  TK(1).  CP) 

35  KT  *  KT  +  1 

REAEKLI  *.  END-33)  TEM,  CP.  XI,  HMH298.  ENTROP,  X2, 

$  X3,  ENTHALP.  X4 
TK(KT)  *  TEM 
CPP(KT)  *  CP 
SS(KT)  *  ENTROP 
HH(KT)  *  HMH298 

C  WRITE  (LO .♦)  KT.  TEM,  CP,  ENTROP 

IF  (IEM  JJL  14999 .DO)  GOTO  49 
IF  (IEM  .GE.  15000. DO)  THEN 
HTTEMP  =  TRUE. 

GOTO  60 
END  IF 

49  IF  (CP  £Q.  999.0D0)  GOTO  60 
CALL  FMATR  (1.  TK(KD,  CP) 

IF  (NOT.  HTTEMP)  GOTO  35 
60  NUM=KT-1 


IF  (HTTEMP)  NUM=KT 

B(l)  *  Y 
B(2)  *  YT 
B(3)  =  YT2 
B(4)  *  YT3 
B(5)  *  YTM 
B(6)  «  YTM2 
B(7)«  YTM3 
A(l,l)  -  SN 
A(2J5)  «  SN 
A(3,6)  -  SN 
A(4,7)  -  SN 
A(5,2)  >  SN 
A(63)  -  SN 
A(7A)  «  SN 
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A(2,l)  -  T 
A(33)«T 
A(4,6)  ■  T 
A(U)-T 
A(53)  -  T 
A(M)«T 
A(3,l)  «  T2 
A02)-T2 
A(4jS)  ■  T2 
A(13)  -  T2 
A(5,4)  «  T2 
A(4,l)  *  T3 
AC23)«T3 
A(l,4)*T3 
A(3^)-T3 
A(5,l)  *  TM 
A06)«TM 
A(3,7)  *  TM 
A(13)  «  TM 
A(63)  *  TM 
A(73)  «  TM 
A(6,l) «  TM2 
A(2,7)  =  TM2 
A(53)*TM2 
A(l,6)  «  TM2 
A(73)  =  TM2 
Aa.l)  «  TM3 
A(53)  -  TM3 
A(l,7)  =  TM3 
A(63)  *  TM3 
A(2,4)  =  T4 
A(33)*T4 
A(43)  «  T4 
A(3,4)  *  T5 
A(43)»T5 
A(4,4)  =  T6 
A(5,7)  *  TM4 
A(63)  *  TM4 
A(73)  *  TM4 
A(6,7)  «  TM5 
Af73)-TM5 
A(7.7)  «  TM6 

CALL  SOLVEX(7,  A,  B,  KX,  0,  KERR) 

IF  (KERR  IE.  0)  GOTO  66 
64  STOP 
66  X-  TFDC/1000. 

B(8)  *  HOFJ  -  TFIX  *  831451  *  ( (  (  B(4)  *  X  /4.0  +  B(3V 
A3X>)*X  +  B(2) / 2.0 ) * X  ♦  B(l)  -  ( ( B(7)  /  (X  *  2.0) 

A  +  B(6))  /  X  -  B(5)  *  DLOG(X)  )  /  X ) 
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BC9)  *  ENTRPO  -  8.31451  *  (  (  (  B(4)  *  X  /  3.0  +  B(3)  /  2D  ) 

A*X  +  B<2))*X  ♦  B(l)  *  DLOCKX)  -  (  (  B(7)  /  (3D  •  X)  ♦ 

A  B(6)  /  2.0 )  /  X  +  B(5) )  /  X  ) 

WRITE  (LO.  70)  XNAME,  PHASE 
70  FORMATC/5X,  A10,  A  5X,  13HThe  phase  ii ,  A6,  / ) 

WRITE  (JLQJ2)  TFDC,  TSTART,  ENTRPO,  HOFJ 
22  FORMAT  (’  Rt  fixed  at  T  -’.FlU/ 

$  ’  Hi  sums  mT  *\F1U/ 

S  ’  Entropy  at  T  Fix  *\  F11.2/ 

$  '  Enthalpy  M  T  Fix  *\  F10.I  /  ) 

C  ***************************************** 

C  WRITE  (4. 71)  XNAME,  PHASE 
C  71  FORMATC/5X.  A10,  Phase  -  \  A6/) 

C  DO  9877  UK*1,NUM,4 
C  IIK1  *  UK 
C  UK2  -  UK1+3 
C  IF  (DK2  .GT.  NUM)  UK2  -  NUM 
C9877  WRITE  (4. 9878)  (U,  TK(U),  U=UK1JJK2) 

C9878  FORMAT  (4(I5,y,F8.1.3X)) 

C  *************************************************** 

WRITE  (LO.  80) 

80  FORMAT  (*  The  data  cards  for  die  TIGER  code  are  (units  =  Joules,  ’ 
$  ’  mols,  &  K):  7) 

1*1 

WRITE  (LO.90)  XNAME,  PHASE,  I,  B(l),  B(2),  B(3) 

WRITE  (3, 95)  XNAME,  PHASE,  I,  B(l),  B(2),  B(3) 

WRITE  a  95)  XNAME,  PHASE,  I,  B(l),  B(2),  B(3) 

1*2 

WRITE  (LO.90)  XNAME,  PHASE,  I,  B(4),  B<5),  B(6) 

WRITE  0,  95)  XNAME,  PHASE,  l  B<4),  B(5),  B(6) 

WRITE  (2, 95)  XNAME,  PHASE,  I,  B(4),  B(5),  B(6) 

1*3 

WRITE  (LO.90)  XNAME,  PHASE,  I,  B(7),  B(8),  B(9) 

WRITE  (3. 95)  XNAME,  PHASE,  I,  B(7),  B{8),  B(9) 

B8  *  B(8)/4.184D0 
B9  *  BC9V4.184D0 

WRITE  (2, 95)  XNAME,  PHASE,  I,  B(7),  B8,  B9 
90  FORMATT5X,  5HSTR, ,  A6 , 7,  A6,  7  42, 3(2H,  .1PD14.7)) 

95  FORMAT(  5HSTR, ,  A6  ,  7.  A6,  7, 12,  3(2H,  ,1PD14.7)) 

C  WRITE  (4. 110) 

110  FORMAT^  T  CP  S 
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$  ’H-H298  MU/RT') 

120  F0RMAT(1XP8AF15  J  ^5.3^4.3^14.4) 


T  «  298.15 

CALL  CALC  CTJiOFJ.CP^ JD 
C  WRITE  (4. 120)  T,  CP.  S,  H 
T«  300.0 
DO  130  I  *  1. 18 
CALL  CALC  O' JI0FJ,CP3JD 
C  WRITE  (4. 120)  T.  CP,  S.  H,  UU 
130T-T+  100.0 

C  WRITE  (LO,  70)  XNAME,  PHASE 
C  WRITE  (LO.  110) 

DO  140  I  *  1, 30 

CALL  CALC  (T.HOFLCP.S.H) 

C  WRITE  (4, 120)  T,  CP,  S,  H,  UU 
140  T  «  T  ♦  100.0 
C  GOTO  10 
T»  6000.0 
DO  150  I »  1.9 
CALL  CALC  (T.HOFJ.CP.S.H) 

C  WRITE  (4, 120)  T,  CP,  S,  H,  UU 
150  T  *  T  ♦  1000.0 

CALL  GOODF5 

GOTO  10 

999  WRITE  (LO.  *)  TEM.  CP.  XI.  HMH298.  ENTROP,  X2 
WRITE  (4.  *)  TEM,  CP.  XI,  HMH298,  ENTROP.  X2 
33  STOP 

END 

SUBROUTINE  SOLVEX(N.  A,  B,  KX,  K.  KERR) 

C 

C 

IMPLICIT  REAL*8  (A-H.O-Z) 

DIMENSION  A(7,7),  B(9),  KX(7),  T(7) 

EPS  «  OX) 

KERR  >  0 
OF  (K)  10,  iai60 
10  DO  20  J  *  1,  N 
20  KX(J)  *  1 
DO  ISO  I  *  1,  N 
IM1 - 1  -  1 
0F0M1)  50,50,30 
30  DO  40  l.IMl 
KJ  *  KX(J) 

Q  «  A(L  KJ) 
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JP1  -  J  ♦  1 
DO  40  L  ■  JP1,  N 
KL«  KX(L) 

40  AO.  KL)  -  AO.  KL)  -  Q  *  A(J,  KL) 
SOKI«KX(D 
Q  -  DABS(Aa  Kl)) 

IF  (N  -  D  100,  100,  60 
gol-i 
ipi  -i+ 1 
DO  80  J  *  IP1,  N 
KJ  «  KX(J) 

AKJ  -  DABS(Aa  KJ» 

IF  (AKJ  -  Q)  80,  8a  70 
70  Q-  AKJ 
L  *  J 

80  CONTINUE 
IF  (L  -  D  90.100,90 
90  J  «  KX(L) 

KX(L)«KX(I) 

KX(D- J 

100  IF  (Q-  EPS)  liana  120 
110  KERR  -  1 
RETURN 
120  Kl «  KX(I) 

Q  —  Ad,  Kl) 

IF(N-D  150,150.130 
130 DO  140  J-IP1.N 
KJ  —  KX(J) 

140  Aa  KJ)  -  Aa  KJ)  /  Q 
150  CONTINUE 
160  DO  200  I-l.N 
Q«B(D 
IM1 -  I  -  1 

IF0M1)  190,190,170 
170  DO  180  J-1.IM1 
kj  sb  Kxrn 

180  Q  -  Q  -  AG.  KJ)  *  B(J) 

190  Kl «  KX(I) 

200  B(D  *  Q  /  Aa  Kl) 

DO 230  IC-l.N 
I  -  N  -  IC  +  1 
Q-BO) 

IF  (N  -  U  230,230,210 
210 IP1  -  I  +  1 

DO  220  J-IPl.N 
KJ-KXCD 

220Q-Q-T(D*  AaKD 
230  TO) -Q 

DO  240  J-l.N 


KJ  ■  KX(J) 

240  B(KJ)  -  T( J) 

RETURN 

END 

SUBROUTINE  CALC  (TJIOFJ.CPAH) 

IMPUCTT  REAL*8  (A-H.O-Z) 

COMMON  /COEFS  /B(9) 

TH  -  T  / 1000.0 
TH2  *  TH  *  TH 
TH3 -TH2  •  TH 
RTH«  1.0/TH 
RTH2-RTH  *RTH 
RTH3  ■  RTH2  *  RTH 

CP  KB<1)  +  B(2)*TH  +  B(3)*TH2  +  B(4)*TH3  +  B(5)*RTH  +  B(6)*RTH2 
$  +  B(7)*RTH3)*  831451 
X-DLOGCTH) 

S  *  8.31451  *  (  B(1)*X  +  Bam*  +  B(3)*TH2/2.0  +  B(4)*TH3/3.0 
$  -  B(5)*RTH  -  B(6)*RTH2/2.0  -  B(7)*RTH3/3.0  )  +  B(9) 

H*(T*8ji45i  *(bo)+ Barra/2.0  +  B(3rra2/3.o+ 

$  B(4)*TH3/4.0  +  B(5)*RTH*X  -  B(6fRra  -  B(7)*RTH3/2.0) 

$  +  B(8)  *  HOFJ) 

C  UU«((H*  1000.0 +  HOFD/T-S)  /  831451 
RETURN 
END 

SUBROUTINE  FMATR  (IOPT.  TEM.  CP) 

IMPUCIT  REAL*8  (A-H.O-Z) 

COMMON /LSTSQ/T,  T2,  T3,  T4,  T5,  TO,  TM2,  TM3,  TM4,  TM5,  TM6, 
$  Y,  YT2,  YT3,  TM,  SN,  YTM2,  YT,  YTM,  YTM3 

IF  (IOPT  JNE.  0)  GOTO  10 

Y  *  0.0 

YT  *  0.0 

YT2  *  0.0 

YT3  *  0.0 

YTM  *  OX) 

YTM2  *  0.0 
YTM3  *  0.0 
SN  *  0.0 
T*0X) 

T2-0.0 
T3-0.0 
T4-0.0 
T5-0.0 
TO*  OX) 
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TM-OO 

TM2-O0 

TM3-O0 

TM4-O0 

TMS-OO 

TM6-0.0 

GOTO  20 

10  X  *  CP  /  8.31451 
C  CHECK  print 

C  WRITE  (4. 9999)  TOM.  CP.  T,  X.  Y 
C9999  FORMAT  (’  TOM  ft  Q»  -  \  2F10.2,  3F15.4) 


TH  ■  TOM  /  1000.0 
TH2  «  TH  •  TH 
TH3  «  TH  *  TH2 
RTH  «  1.0/TH 
RTH2  «  RTH  *  RTH 
RTH3  -  RTH  *  RTH2 
Y  *  Y  +  X 
YT  *  YT  +  X  *  TH 
YT2  *  YT2  +  X  *  TH2 
YT3  «  YT3  +  X  *  TH3 
YTM  «  YTM  +  X  *  RTH 
YTM2  «  YTM2  +  X  *  RTH2 
YTM3  »  YTM3  +  X  *  RTH3 
SN  *  SN  +  1.0 
T-T  +  TH 
T2«T2  +  TH2 
T3  »T3  +  TH3 
T4  «  T4  +  TH2  *  TH2 
T5-T5+  TH2  *TH3 
T6  *  T6  +  TH3  *  TH3 
TM  «  TM  +  RTH 
TM2  «TM2  +  RTH2 
TM3  *  TM3  +  RTH3 
TM4  «TM4  +  RTH2  *  RTH2 
TM5  *TM5  +  RTH2  *  RTH3 
TM6  «  TM6  +  RTH3  *  RTH3 
20  RETURN 
END 
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Intentionally  left  blank. 
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PROGRAM  TK3F1T5X 


Appendix  L. 

Listing  of  Program  TK3FTT5X 


This  version:  26  January  1993 

Specifically  adapted  to  creating  coefficients  for  EXBLAKE.  This 
program  is  die  fourth  of  a  aeries  of  programs  named  either  Thjrfii 
or  nGFlTlX  (i «  1,4).  This  program  is  specifically  designed  to 
get  its  input  from  tables  created  by  expanding  the  extended-range 
fitting  coefficients  given  to  ARL/WTD/APCB  by  Ms.  Bonnie  McBride 
(NASA  Lewis  Research  Center).  The  expansion  was  done  by  Program 
EXCALCTH4  (q.v.). 

This  program  has  all  of  the  features  of  its  predecessors  (use  of  the 
*$’  sentinel:  computation  of  measures  of  goodness-of-fit);  but  it 
does  not  include  die  cut-off  at  IS  000  K. 

In  this  version,  die  output  files  have  been  dunged  from  TIGHT-OUT 
and  HGFIT.AUX  to  TIGFITSX.OUT  and  HGFrrSJtAUX. 

C  This  one  has  all  of  die  features  of  its  predecessors  (use  of  the  *$* 

C  sentinel;  computation  of  measures  of  goodnessof-fit);  but  it  does 
C  not  indude  the  cutoff  at  IS  000  K. 

C  For  further  important  details,  see  the  Comment  lines  at  the  beginning 
C  of  TIGFIT4X. 

IMPLICIT  REAL*8  (A-H.O-Z) 

C 

DIMENSION  A(7.7).  B(9),  KX(7) 

DIMENSION  TK(500),  CPP(500).  SS(500),  HH(500) 

CHARACTER  XNAME*10,  PHASE*d,  SENT*1 
PARAMETER  (N6=0) 

LOGICAL  HITEMP 

COMMON  /COEFS/  B.  NUM,  NNUM 

COMMON  /FUNCS/  TK.  CPP,  SS,  HH,  TFIX,  TSTART,  ENTRP0,  HOFJ, 
$  XNAME,  PHASE 
COMMON  /IOSTUF/  LI,  LO 

COMMON  /LSTSQ/  T.  T2,  T3,  T4.  TS.  T6.  TM2,  TM3,  TM4,  TMS,  TM6, 
$  Y,  YT2,  YT3,  TM,  SN.  YTM2,  YT,  YTM,  YTM3 

DATA  PHASE/' GAS  ’/ 

C 
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LO*  6 
U-S 

OPEN  (UJFILE«TEMPSJN,fSTATUS«,OLD,) 

OPEN  (3,FILE*,TIGFIT5XOUT,,STATUS«,UNKNOWN') 

OPEN  (LOJTLE-TKHTTSX-AUX’^TATUS-’UNKNOWN’) 

OPEN  (4mE*,CHEKK5X’^TATUS»,UNKNOWN’) 

OPEN  (2JTLE-*CALORIC’.ACCESS«’APPEND,,STATUS=’lJNKNOWN’) 
WRITE  a  81) 

81  FORMAT  (’  The  data  cards  for  the  TIGER  code  are  (units  *  calories’ 

$  \  mols,  &  K):  ’) 

C 

C  CALL  LOGO 
C 

REWIND  LI 

READ  (LI.  *.  END-33)  TFIX,  TSTART 
CLOSE  (LI) 

OPEN  (Ujn^’EXO^LCTH.OUT’^TATUS-’OLD’) 

REWIND  U 

10  READOL1020.END-33)  SENT.  XNAME 
1020  FORMAT(Al.  A10  ) 

IF  (SENT  J®.  ’$’)  GOTO  10 

WRITE  (N6.2008)  XNAME 
2008  FORMAT(lX.Trocessing  **>  ’.A10) 

C  Initialize  die  elements  that  will  make  up  foe  matrix 
CALL  FMATR  (0,  XI.  X2) 

C  Read  foe  first  line  of  data. 

13  READ(U.*,  END-33)  TEM.  CP.  ENTRP0,  HMH298.  HOFJ 
C  WRITE  (LO.*)  "Input  1".  TEM.  CP.  ENTROP.  HMH298.  HOFJ 
C  Is  temperature  foe  pinning  temperature?  If  not,  read  another  line. 

IF  OEM.  LT.  TFIX-0.05  .OR.  TEM  .GT.  TFIX+O.Q5  )  GOTO  15 
C  Yes.  it  is  the  pinning  T.  Is  T  Fix  «  T  Start? 

C  If  it  is,  then  start  saving  data  and  continue  to  read  foe  table 
20  IF  (TEM  .GE.  TSTART-0.Q5  .AND.  TEM  XE.  TSTART+0.05  )  THEN 
KT«  1 

TK(1)  -  TSTART 
CPP(1)  =  CP 
SS(1)  =  ENTRPO 
HH(1)  -  HMH298 
GOTO  35 

END  IF 

C  No.  T  Fix  wasn’t  T  Start  Keep  looking  for  T  Start 
25  READ  (LI,*,  END-33,  ERR-999)  TEM,  CP,  ENTROP.  HMH298,  ENTHALP 
C  WRITE  (LO,*)  "Input2\  TEM,  CP,  ENTROP,  HMH298,  ENTHALP 
IF  (TEM  XT.  TSTART-0.05  .OR  .  TEM  .GT.  TSTART+0.05  )  GOTO  25 
C  At  last  T  Start  has  been  found:  save  values  and  keep  reading  table. 

29KT-  1 
TK(1)  «  TSTART 
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CPP(1)-CP 
SS(1)  -  ENTROP 
HH(1)  -  HMH298 
HTIEMP*  .FALSE. 

CALL  FMATR  (1,  TK(1),  CP) 

35  KT  *  KT  ♦  1 

READ(LL*.  END«33)  TEM.  CP.ENTROP,  HMH298.  ENTHALP 
C  WRITE  (LO.  *)  "Input3",  TEM,  CPfNTROP,  HMH298,  ENTHALP 
TK(KT)  *  TEM 
CPP(KT)  *  CP 
SS(KT)  -  ENTROP 
HH(KT)  *  HMH298 

C  WRITE  (LO  ,*)  TEM,  CP.  ENTROP 

IF  (TEM  UE.  14999.)  GOTO  49 
DP  OHM  £Q.  15000.)  THEN 
HTIEMP*  .TRUE. 

GOTO  60 
END  IF 

49  OF  (CP  EQ.  999.0)  GOTO  60 
CALL  FMATR  (1,  TK(KT>,  CP) 

IF  (.NOT.  HTIEMP)  GOTO  35 

C  The  end  of  the  table  hu  been  reached.  Is  the  list  temperature 
C  equal  to  or  greater  than  10  000  K? 

60  NUM-KT-1 


IF  (HTIEMP)  NUM=KT 

IF  (TK(NUM)  GE.  1.D4)  THEN 
NNUM=0 
GOTO  61 
END  IF 

61  CALL  EXTEND  (TK,  CPP.  NUM,  NNUM) 
B(1)«Y 
B(2)»YT 
B(3)  *  YT2 
B(4)«  YT3 
B(5)  ■  YTM 
B(6)  *  YTM2 
B(7)  *  YTM3 
A(l,l)  *  SN 
AQJS)  *  SN 
A(3,6)  =  SN 
A(4,7)  *  SN 
M52)  *  SN 
A(63)  *  SN 
ACM)  *  SN 
A(2,l)  *  T 
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AOS)-T 
A(4jS)  rn  T 
A(U)  «  T 
A(53)  *  T 
A(6,4)  «  T 
A(3,l)  «  T2 
A(23)-T2 
A(43)-T2 
A(13)-T2 
A(5,4)  »  T2 
A(4,1)«T3 
A(23)  *  T3 
A(l,4)  «  T3 
A(33)  «  T3 
A(5,1)«TM 
A(2,6)  »  TM 
A(3,7)  »  TM 
A(13)  «  TM 
A(63>  *  TM 
A(73)  *  TM 
A(6,l)  «  TM2 
A(2,7)  «  TM2 
A(53)  *  TM2 
A(M)  -  TM2 
A(73)  =  TM2 
A(7,l)  =  TM3 
A(5,6)  *  TM3 
A(l,7)  -  TM3 
A(63)  *  TM3 
A(2,4)  *  T4 
A(33)  *  T4 
A(43)  *  T4 
A@A)  *  T5 
A(43)«T5 
A(4,4)  *  T6 
A(5,7)  *  TM4 
A(6,6)  -  TM4 
A(73)  *  TM4 
A(6,7)  -  TM5 
A(7,6)  -  TM5 
A(7,7)  *  TM6 

CALL  SOLVEX(7,  A,  B,  KX,  0,  KERR) 

IF  (KERR  JLE.  0)  GOTO  66 
64  STOP 
66  X*  TFIX/1000. 

B(8)  «  HOF)  -  TFIX  •  831451  *  (  (  (  B(4)  *  X  /  4.0  +  B(3)/ 
A  3.0 )  * X  +  B<2)/2.0)*X  +  B(l) - ( (  B(7) / (X  *  2.0) 

A  ♦  B(6))  /  X  -  B(5)  *  DLOG(X)  )  /  X  ) 

B(9)  ■  ENTRPO -  8.31451  *(((B(4)  *X/3.0+B(3)/ 2.0) 
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A  *  X  +  B(2)  )  *  X  ♦  B(l)  *  DLOG(X)  -  (  (  B(7)  /  (3.0  *  X)  ♦ 
AB(6)/2.0)/X  +  B(5)  )  /  X  ) 

WRITE  (LO,  70)  XNAME.  PHASE 
70  FORMAT(5X.  A10.  /,  5X,  13HThc  phase  is .  A6,  /  ) 

WRITE  (L0.22)  TFIX.  TSTART,  ENTRPO,  HOFJ 
22  FORMAT  (’ Fit  fixed  at  T  «’,F11.2/ 

$  ’Fit  starts  at  T  »’,F11.2/ 

$  ’  Entropy  at  T  Fix  -’.F11J/ 

$  ’  Enthalpy  at  T  Fix  «\  F10.1 /) 


C  ••M*****»****MM«M*M*MMMM«*M**M 

C  WRITE  (4. 71)  XNAME,  PHASE 
C  71  F0RMATC/5X,  A10,  Thasc  =  ’,  A6/) 

C  DO  9877  UK-1.NUM.4 
C  UK1  =  UK 
C  UK2  =  UK  1+3 
C  IF  (UK2  .GT.  NUM)  UK2  «  NUM 
C9877  WRITE  (4.  9878)  (U.  TK(U),  U=UK1JJK2) 
C9878  FORMAT  (4(I5,’)*JF8.1.3X)) 

c  **M****M********«*«M*M*M«**Hmm* 


WRITE  (LO.  80) 

80  FORMAT  (’  The  data  cards  for  die  TIGER  code  are  (units  =  Joules,  * 
$  ’  mols,  &  K):  ’/) 

1=1 

WRITE  (LO.90)  XNAME,  PHASE.  I.  B(l),  B(2),  B(3) 

WRITE  (3,  95)  XNAME,  PHASE,  I,  B(l),  B(2),  B(3) 

WRITE  (2,95)  XNAME.  PHASE,  I,  B(l),  B(2),  B(3) 

1  =  2 

WRITE  (LO.90)  XNAME,  PHASE,  I,  B(4),  B(5),  B(6) 

WRITE  (3.  95)  XNAME,  PHASE,  I.  B(4),  B(5),  B(6) 

WRITE  (2.  95)  XNAME,  PHASE,  L  B(4),  B(5),  B(6) 

1  =  3 

WRITE  (LO,90)  XNAME,  PHASE,  I,  B(7),  B(8),  B(9) 

WRITE  (3,  95)  XNAME,  PHASE,  I,  B(7),  B(8),  B(9) 

B8  =  B(8)/4.184D0 
B9  =  B(9)/4.184D0 

WRITE  a  95)  XNAME,  PHASE,  I,  BC7),  B8,  B9 
90  FORMAT(2X,  5HSTR, ,  A6,  A6,  J13(2H.  .1PD15.7)) 

95  FORMAT(  5HSTR, ,  A6,  A6,  I1,3(2H,  .1PD15.7)) 

C  WRITE  (4, 110) 

110  FORMAT^  T  CP  S’ 

$  ’H-H298  MU/RT’) 

120  FORMAT(lXJ8ZF15.34f15.3JJ14.34J14.4) 
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T  -  298.15 

CALL  CALC  (T.HOFJ.CP.S.H) 

C  WRTIE  (4, 120)  T,  CP,  S,  H 
T-  300.0 
DO  130  I  -  1.  18 
CALL  CALC  (T*HOFJ,CPtSJH) 

C  WRITE  (4, 120)  T,  CP,  S,  H,  UU 
130  T*  T  +  100.0 

C  WRTTE  (LO,  70)  XNAME,  phase 
C  WRITE  (LO.  110) 

DO  140  I  «  1.  30 

CALL  CALC  CT»HOFJ,CP,S,H) 

C  WRTTE  (4. 120)  T,  CP,  S,  H,  UU 
140  T  «  T  +  100.0 
C  GOTO  10 
T«  6000.0 
DO  150  I«  1,9 
CALL  CALC  (TJlOFJ.CP^Ji) 

C  WRTTE  (4. 120)  T,  CP.  S.  H.  UU 
150  T  *  T  +  1000.0 
CALL  GOODF5X 
GOTO  10 

999  WRTTE  (LO.  *)  TEM,  CP.  XI.  HMH298,  ENTROP,  X2 
WRTIE  (4.  *)  TEM,  CP.  XI.  HMH298,  ENTROP.  X2 
33  STOP 
END 

SUBROUTINE  CALC  fl\HOFJ,CP,S,H) 

IMPLICIT  REAL*8  (A-H.O-Z) 

COMMON  /COEFS  /B(9) 

TH  *=  T  /  1000.0 
TH2«TH*TH 
TH3  =  TH2  *  TH 
RTH  =  1.0/TH 
RTH2  *  RTH  *  RTH 
RTH3  «  RTH2  *  RTH 

CP  «(B(1)  +  B(2)*TH  +  B(3)*TH2  +  B(4)*TH3  +  B(5)*RTH  B(6)*RTH2 
$  +  B(7)*RTH3)*  8.31451 
X«DLOG(TH) 

S  *  8.31451  *  (  B(1)*X  +  B(2)*TH  +  B(3)*TH2/2.0  +  B(4)*TH3/3.0 
$  -  B(5)*RTH  -  B(6)*RTH2/2.0  -  B(7)*RTH3/3.0  )  +  B(9) 

H  «  (  T  *  8.31451  *  (  B(l)  +  B(2)*TH/2.0  +  B(3)*TH2/3.0  + 

$  B(4)*TH3/4.0  +  B(5)*RTH*X  -  B(6)*R1H2  -  B(7)*RTH3/2.0) 

$  +  B(8)  -  HOFJ  ) 

RETURN 

END 
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SUBROUTINE  FMATR  (IOPT,  TEM,  CP) 
IMPUCIT  REAL*8  (A-H.O-Z) 


COMMON  /LSTTSQ/  T,  T2,  T3,  T4,  T5.  T6,  TM2, 1X3,  TM4,  TM5,  TM6, 
$  Y,  YT2,  YT3,  TM,  SN,  YTM2,  YT,  YTM,  YTM3 

IF  OOPT  .NE.  0)  GOTO  10 

Y  =  0.0 

YT  *  0.0 

YT2  *  0.0 

YT3  *  0.0 

YTM  *  0.0 

YTM2  *  0.0 

YTM3  =  0.0 

SN  =  0.0 

T  *  0.0 

T2  m  0.0 

T3  *=  0.0 

T4  =  0.0 

TS  *  0.0 

T6  =  0.0 

TM  =  0.0 

TM2  =  0.0 

TM3  *  0.0 

TM4  *  0.0 

TM5  *  0.0 

TM6  *  0.0 

GOTO  20 

10  X  =  CP  /  8.31431 


C  CHECK  print 

C  WRITE  (4, 9999)  TEM,  CP,  T,  X,  Y 
C9999  FORMAT  (*  TEM  &Cp  =  ',  2F102,  3F15.4) 


TH  *  TEM  /  1000.0 
TH2  =  TH  *  TH 
TH3  *  TH  *  TH2 
RTH  =  1.0/TH 
RTH2  =  RTH  ♦  RTH 
RTH3  =  RTH  *  RTH2 
YsY  +  X 
YT  =  YT  +  X  *TH 
YT2  =  YT2  +  X  *  TH2 
YT3  -  YT3  +  X  *  TH3 
YTM  *  YTM  +  X  *  RTH 
YTM2  «  YTM2  +  X  *  RTH2 
YTM3  =  YTM3  +  X  ♦  RTH3 
SN  *  SN  +  1.0 
T  =  T  +  TH 
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T2-T2  +  TH2 
T3-T3  +  TH3 
T4  «  T4  ♦  TH2  *  TH2 
TS-T5  +  TH2*TH3 
T6-T6  +  TH3  * TH3 
TM-TM  +  RTH 
TM2  *TM2  +  RTH2 
TM3-TM3  +  RTH3 
TM4  »  TM4  +  RTH2  *  RTH2 
TM5  «TM5  +  RTH2  *  RTH3 
TM6  «  TM6  +  RTH3  *  RTH3 
20  RETURN 
END 

SUBROUTINE  SOLVEX(N,  A,  B,  KX,  K,  KERR) 
C 

IMPUCTT  REAL*8  (A-H.O-Z) 

DIMENSION  A(7.7).  B(9).  KX(7),  T(7) 

EPS  «  0.0 
KERR  *  0 
IF  (K)  10,10,160 
10  DO  20  J  *  1,  N 
20  KX(J)  »  J 
DO  150  I  *  1,  N 
IM1  *1-1 
IF  (IM1)  50.50.30 
30  DO  40  J*l,  IM1 
KJ  *  KX(J) 

Q“  A(I,  KJ) 

JP1  -  J  +  1 
DO  40  L  «  JP1,  N 
KL  =  KX(L) 

40  A(I,  KL)  =  AG.  KL)  -  Q  *  A(J,  KL) 

50KI  =  KX(I) 

Q  *  DABS(A(L  KI)) 

IF  (N  -  D  100.100.60 
60  L  =  I 
IP1  «  I  +  1 
DO  80  J  *  IP1,  N 
KJ«KX(J) 

AKJ  *  DABS(A(L  KJ)) 

IF  (AKJ  -  Q)  80.80.70 
70  Q  *  AKJ 
L»J 

80  CONTINUE 
OF  (L  - 1)  90.100.90 
90  J  *  KX(L) 

KX(L)  *  KX(I) 

KX(I)>  J 

100  IF  (Q- EPS)  110,  110,120 
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110  KERR  -  1 
RETURN 
120  K1  *  KX(Q 
Q  -  AO,  KI) 

IF  (N  -  D  150,150.130 
130 DO  140  J-IP1.N 

ifi « ifYn^ 

140  AG,  KJ)  »  A(I,  KJ)  /  Q 
150  CONTINUE 
160  DO  200  I-l.N 
Q*B(I) 

IM1 « I  -  1 

IF0M1)  190.190.170 
170  DO  180  J-l.IMl 
1CT  w  icxm 

180  Q  -  Q  -  AO.  KJ)  *  B(J) 
190  KI  -  KX(I) 

200  B(D  -  Q  /  A(I.  KD 
DO  230  IC-l.N 
I «  N  -  IC  +  1 
Q«B(I) 

BF(N  -D  23a  230.  210 
210 IP1  =  I  +  1 

DO  220  J  *  IP1,  N 

KJ  s  KJCCJf^ 

220  Q  »  Q  -  T(J)  *  A(I,  KJ) 
230  TO) -Q 

DO  240  J=  l.N 
KJ  *  KX(J) 

240  B(KJ)  *  T(J) 

RETURN 

END 
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Code  G301.D.  Wilson 
Code  G32,  Guns  Systems  Division 
Code  E23,  Technical  Library 
Dahlgren.  VA  22448-5000 

2  Commander 

Naval  Surface  Warfare  Center 
Crand  Division 
ATTN:  Code  4052. 

Steve  Backer 
D.  Holmes 

Bldg.  108 

Crane.  IN  47522-5000 

2  Commander 

Naval  Underwater  Systems  Center 
Energy  Conversion  Dept 
ATIN:  Code  5B331,  R.  S.  Lazw 
Technical  Library 
Newport,  RI 02840 

1  Commander 

Naval  Surface  Weapon  Center 
Indian  Head  Division 
ATTN:  Code  270P1,  Mr.  Ed  Chan 
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ATTN:  Jeff  Brown 

312  Mechanical  Engineering  Bldg. 
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Alexandria.  VA  22312 

2  General  Dynamics  Land  Systems 
ATTN:  Mr.  F.  Lunsford 

Dr.  M.  Weidner 
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1  University  of  Illinois  at  Chicago 
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Dr.  O.  Cook 
Fart  Halstead 
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2  ftaunhota-Imtitnt  fder  Chemische  Technologic 
ATTN:  Dr.  F.  Volk 

Mr.  H.  Batbdt 
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